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ABSTRACT 

An  experimental  study  of  laminar  flow  heat  transfer  of  an  ethylene  glycol/water 
mixture  in  an  electrically  heated  horizontal  tube  using  wire  mesh  (HEATEX)  and 
twisted  tape  inserts  was  investigated.  Twelve  thermocouples,  inserted  in  the  tube  wall 
at  four  longitudinal  locations,  enabled  a  mean  inside  experimental  heat-transfer 
coefficient  to  be  accurately  measured.  A  constant  wall  heat  flux  boundary  condition 
was  placed  on  the  wall  by  wrapping  six  200  W  flexible  heater  tapes  tightly  around  the 
tube.  The  ethylene  glycol/water  mixture  provided  a  coolant  Reynolds  number 
between  200-5000  and  a  Prandtl  number  between  30-140.  Two  smooth  inside 
diameters  and  a  roped  tube  profile  were  tested  with  and  without  the  inserts.  Heat- 
transfer  correlations  for  tubes  without  inserts  were  developed  and  compared  with 
theory  for  both  thermally  and  hydrodynamically  developing  flow.  Correlations  were 
also  developed  for  the  two  types  of  inserts.  Nusselt  numbers  for  fully  developed  flow 
were  found  to  be  a  function  of  Reynolds  and  Prandtl  numbers  for  the  wire  mesh 
insert  and  a  function  of  tape  twist  ratio,  Reynolds  and  Prandtl  numbers  for  the 
twisted  tape  insert.  Heat  transfer  enhancements  of  over  7  for  the  wire  mesh  insert 
and  over  4  for  the  twisted  tape  insert  at  high  Reynolds  numbers  were  obtained  over 
the  empty  tube.  By  using  these  correlations  in  conjunction  with  earlier  obtained 
refrigerant  condensation  data  (using  the  same  tubes,  inserts,  and  coolant),  more 
accurate  values  of  the  outside  condensation  heat-transfer  coefficients  were  obtained. 
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I.   INTRODUCTION 

A.   BACKGROUND 

With  the  increasing  complexity  of  shipboard  weapons  and  combat  systems  used 
aboard  United  States  Naval  vessels,  the  need  for  more  energy  efficient,  lightweight, 
and  high  capacity  refrigeration,  air  conditioning,  and  chilled  water  systems  has  grown 
tremendously.  One  way  to  meet  this  need  is  to  redesign  the  refrigeration  condenser. 
In  redesigning  these  condensers,  changes  in  condenser  tube  material  can  be 
combined  with  various  heat  transfer  enhancement  methods  to  achieve  a  lighter 
weight,  more  efficient  condenser  tube. 

Emphasis  on  refrigerant  condensation  research  reached  a  much  higher  level  of 
importance  due  to  the  pressing  environmental  issues  of  ozone  depletion  and  global 
warming.  In  September  1987,  an  international  conference  was  held  in  Montreal, 
Canada,  in  which  24  nations  representing  the  United  Nation's  Environment  Program 
(UNEP),  signed  the  Montreal  Protocol  on  Substances  that  Deplete  the  Ozone  Layer. 
This  agreement  called  for  a  near  term  freeze  on  the  production  and  consumption  of 
chlorofluorocarbons  (CFCs),  the  major  contributor  to  ozone  depletion 
[Ref.  1].  In  August  1988,  the  Environmental  Protection  Agency  (EPA) 
adopted  the  regulations  of  the  Montreal  Protocol  under  the  Clean  Air  Act  [Ref.  2]. 
These  regulations  required  a  freeze  in  both  production  and  consumption  of  CFC-11, 


12,  113,  114  and  115  at  1986  levels.  The  U.S.  Navy  currently  uses  CFC-12  and  CFC- 
114  in  its  shipboard  air  conditioning  and  refrigeration  systems.  In  June  1990  at  a 
progress  meeting  held  in  London,  England,  UNEP  delegates  agreed  to  even  larger 
reductions,  calling  for  a  50%  reduction  by  1995,  85%  reduction  by  1997,  and  for 
CFCs  to  be  totally  phased  out  by  the  year  2000  [Ref.  3].  In  the  spring  of  1992, 
President  Bush  moved  up  the  CFC  phase  out  deadline  to  1995. 

Because  of  the  Navy's  commitment  to  the  CFC  phase  out  policy  and  the 
increasing  need  for  lighter  weight,  more  efficient  refrigeration  systems,  refrigerant 
condensation  research  has  surfaced  to  the  forefront  of  Naval  Research. 

B.      REFRIGERANT  CONDENSATION  RESEARCH  AT  NPS 

In  order  to  examine  heat  transfer  characteristics  of  various  condenser  tube 
designs  and  alternative  refrigerants,  a  condenser/evaporator  test  platform  was 
constructed  at  the  Naval  Postgraduate  School.  Construction  was  begun  by  Zebrowski 
[Ref.  4]  and  completed  by  Mabrey  [Ref.  5]  in  1988.  The  test  platform,  shown  in 
Figure  1.1,  consisted  of  an  evaporator  and  a  condenser  section  with  an  associated 
ethylene  glycol/water  mixture  coolant  system. 

The  evaporator  section  was  designed  so  that  in  addition  to  supplying  the 
refrigerant  vapor  to  the  condenser,  it  could  also  be  utilized  to  conduct  bundle  boiling 
experiments.  The  condenser  section  consisted  mainly  of  four  instrumented  condenser 
tubes  and  a  series  of  auxiliary  condenser  coils  used  to  control  pressure  in  the 
apparatus.  Details  of  design,  construction  and  operation  of  the  condenser/evaporator 


test  platform  are  provided  by  Zebrowski  [Ref.  4],  Mabrey  [Ref.  5]  and  Mazzone 
[Ref.  6]. 

Unlike  a  typical  vapor  compression  refrigeration  system  used  aboard  ships,  the 
test  platform  at  NPS  lacked  a  compressor  between  the  evaporator  and  the  condenser. 
This  compressor  is  normally  required  to  increase  the  pressure  (and  hence 
temperature)  of  the  refrigerant  vapor  sufficiently  enough  so  that  it  can  give  up  its 
latent  heat  to  the  ambient  sea  water,  the  primary  coolant  Because  there  was  no 
compressor  in  the  NPS  test  platform,  the  condenser  must  operate  at  the  saturation 
conditions  of  the  evaporator.  Therefore,  a  much  colder  coolant,  a  refrigerated 
ethylene  glycol/water  mixture,  was  needed  to  condense  the  vapor,  leading  to  laminar 
flow  conditions  withia  the  condenser  tubes.  This  in  turn  leads  to  poor  heat  transfer 
due  to  a  large  inside  thermal  resistance. 

The  primary  objective  of  the  condensation  experiments  of  Mabrey  [Ref.  5]  and 
Mazzone  [Ref.  6]  was  to  accurately  determine  outside  heat  transfer  coefficients. 
Briefly,  this  was  accomplished  by  first  determining  the  heat  transfer  rate  to  the 
coolant  by: 

q  =  m  cp  LT      '  (1.1) 

where  m  is  the  coolant  mass  flow  rate  and  AT  the  coolant  temperature  rise. 
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Figure  1.1    NPS  Evaporator/Condenser  Refrigeration  Apparatus. 


The  overall  heat  transfer  coefficient,  U0,  could  then  be  determined  by: 


Ua  =  2 

°      A_  LMTD 


(1.2) 


where  A,,  is  the  outside  surface  area  of  the  tube  and  LMTD  is  the  log  mean 
temperature  difference,  defined  as. 


LMTD  = 
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The  overall  heat  transfer  coefficient  can  also  be  thought  of  in  terms  of  a  sum  of 
thermal  resistances.  The  overall  thermal  resistance  is  given  by: 
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The  overall  thermal  resistance  can  be  measured  experimentally;  the  wall  resistance, 
R^,  is  known,  and  the  inside  resistance  can  be  determined  from  a  chosen  correlation. 
Therefore,  the  outside  heat  transfer  coefficient,  h0,  can  be  calculated  by: 


h    = 
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(1.5) 


There  is  one  major  drawback  to  this  method.  As  mentioned  earlier,  when 
laminar  flow  heat  transfer  conditions  exist  within  the  condenser  tube,  it  leads  to  poor 
heat  transfer  performance  and  therefore  represents  the  dominant  thermal  resistance 
within  the  condenser.  This  inside  resistance  then  becomes  the  controlling  factor  in 


the  determination  of  the  outside  heat  transfer  coefficient  Indeed,  if  the  inside 
resistance  is  a  large  proportion  of  the  overall  resistance,  then  the  accuracy  of  the 
calculated  outside  heat  transfer  coefficient  becomes  uncertain.  Therefore,  the 
accuracy  of  the  outside  heat  transfer  coefficient  relies  heavily  on  an  accurate 
correlation  to  determine  hj. 

An  alternative  is  to  reduce  the  inside  controlling  resistance  by  augmenting  the 
inside  heat-transfer  coefficient  Augmentation  techniques  for  lowering  the  inside 
thermal  resistance  primarily  focus  on  repeatedly  bringing  fresh  fluid  from  the  coolant 
bulk  stream  to  the  heated  inside  surface.  This  results  in  breaking  up  the  thermal 
boundary  layer  and  thereby  reducing  the  thermal  resistance.  The  technique  utilized 
in-the^NPS-condensation  experiments  to  break  up-this^ thermal  boundary  layer  made 
use  of  two  types  of  inserts,  a  twisted  tape  insert  and  a  wire  mesh  (HEATEX)  insert. 
These  inserts  and  their  use  will  be  discussed  in  more  detail  in  later  chapters. 

The  results  of  condensation  experiments  conducted  by  Mazzone  [Ref.  6]  show 
signs  of  inaccuracies  related  to  the  choice  and  accuracy  of  the  inside  heat  transfer 
coefficient  correlation.  For  example,  in  attempting  to  determine  the  outside  heat 
transfer  coefficient  for  a  smooth  tube  with  no  insert,  Mazzone  found  that  the 
correlation  used  to  determine  the  inside  resistance  yielded  values  that  were  greater 
than  the  actual  measured  overall  resistance!  Obviously,  no  sensible  values  of  h0  could 
be  determined.  Similar  problems  existed  (but  to  a  lesser  extent  due  to  the  greatly 
reduced  inside  thermal  resistance)  when  using  various  types  of  inserts.  Mazzone 
concluded  that  "...  laminar  flow  correlations  for  twisted  tape,  HEATEX,  and  no 


insert  conditions  must  be  developed  which  have  better  accuracy  than  those  used  ... 
allowing  for  more  accurate  determination  of  outside  heat  transfer  coefficients." 
Therefore,  the  determination  of  accurate  inside  heat  transfer  coefficients  is 
paramount  in  the  final  outcome  of  the  refrigeration  condensation  project 

C.      OBJECTIVES 

The  main  focus  of  this  thesis  was  therefore  to  accurately  determine  inside  heat 
transfer  coefficients  in  support  of  refrigerant  condensation  experiments  being 
conducted  at  NPS. 

The  specific  objectives  of  this  study  were  to: 


1.  Design  and  build  a  test  apparatus  to  accurately  determine  inside  heat  transfer 
coefficients  for  the  same  tubes  and  inserts  used  in  the  condensation 
experiments. 

2.  Develop  correlations  for  inside  heat  transfer  coefficients  under  laminar  flow 
conditions  both  with  and  without  inserts. 

3.  Using  the  correlations  mentioned  above,  reprocess  existing  condensation  data 
to  determine  more  accurate  values  of  the  outside  heat  transfer  coefficients. 


II.  LITERATURE  SURVEY 

A.      INTRODUCTION 

Laminar  fully  developed  and  developing  forced  convection  flows  have  been 
analyzed  in  great  detail  for  various  boundary  conditions.  The  study  of  heat  transfer 
in  laminar  flow  through  a  closed  duct  was  first  made  by  Graetz  [Ref.  7]  in  1883  and 
later  by  Nusselt  [Ref.  8]  in  1910.  An  excellent  and  thorough  study  of  the  theoretical 
development  (with  both  analytical  and  experimental  solutions)  for  various  laminar 
forced  convection  internal  flow  boundary  conditions  is  given  by  Shah  and  London 
[Ref.  9].  

This  survey  will  be  concerned  with  laminar  flow  heat  transfer  in  fully  developed 
and  developing  flows  subjected  to  both  a  constant  wall  temperature  and  constant 
heat  flux  boundary  condition.  It  will  also  summarize  some  of  the  augmentation 
techniques  used  to  enhance  tubeside  laminar  flow  heat  transfer  and  discuss  the 
details  involved  when  using  twisted  tape  and  wire  mesh  inserts  to  enhance  inside  heat 
transfer. 


B.      LAMINAR  FLOW  IN  CIRCULAR  DUCTS 
1.      Boundary  Layers  and  Entrance  Lengths 

a.      Hydrodynamic  Boundary  Layer 

For  laminar  flow  in  a  circular  tube,  fluid  enters  the  tube  with  a 
uniform  velocity  as  shown  in  Figure  2.1;  but  when  the  fluid  particles  make  contact 
with  the  tube  wall,  viscous  effects  cause  them  to  assume  a  zero  velocity.  These 
viscous  effects,  associated  with  shear  stresses  acting  in  planes  parallel  to  the  tube 
wall,  tend  to  retard  the  motion  of  particles  in  adjoining  fluid  layers,  causing  a 
momentum  boundary  layer  to  develop.  The  momentum  boundary  layer  thickness,  5, 
increases  until  its  growth  is  stopped  by  symmetry  at  the  tube  centerline.  The  viscous 
effects  then  extend  over  the  entire  cross  section  of  the  tube.  The  flow  is  then  said  to 
be  hydrodynamically  fully  developed.  At  this  point,  the  fluid  velocity  distribution  at 
a  given  cross  section  is  independent  of  the  axial  distance,  x  (i.e.,  du/dx=0)  [Refs.  9, 
10,  11,  12]. 

The  hydrodynamic  entry  length  for  laminar  flow  (ReD  <  2300)  is 
defined  as  that  region  where  the  hydrodynamic  boundary  layer  is  developing  and  may 
be  obtained  from  the  approximation  [Ref.  10]: 


<x   ^ 


(2.1) 

D 

lam 


«  0.05ReD 


This  expression  shows  that  for  ReD=2300,  the  entrance  length  is  115  times  the  tube 
diameter,  a  considerable  length  for  most  experimental  facilities. 
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Figure  2. 1    Laminar  Hydrodynamic  Boundary  Layer  Development  in  a  Circular  Tube 
(Courtesy  [Ref.  10]) 


b.     Thermal  Boundary  Layer 

When  a  fluid  with  a  uniform  temperature  distribution  enters  a  circular 
duct  with  a  wall  temperature  that  is  greater  than  the  uniform  fluid  temperature, 
temperature  gradients  in  the  fluid  develop  and  convective  heat  transfer  occurs,  as  shown 
in  Figure  2.2.  The  region  of  fluid  where  these  temperature  gradients  exist  is  called  the 
thermal  boundary  layer.  The  thermal  boundary  layer  thickness,  5t,  grows  until  a 
dimensionless  fluid  temperature  distribution  is  independent  of  x.  The  thermal  boundary 
layer    is    then    considered    fully    developed.    This    point    is    usually    expressed    in 


10 


Sur!»ce  condition 
r-Tt  >T(r,  0)  flj-i 


•7  M  r    ~  r 


L.< 


Tfr.O)        T,     |         ro.o;         r. 


rrr.  o;      rw 


Thermal  entrance  region  |         Fully  developed  region 

I 


0 


Figure  2.2    Thermal  Boundary  Layer  Development  in  a  Heated  Circular  Tube 
(Courtesy  of  [Ref.  10]) 


terms  of  a  non-dimensional  temperature  profile  given  by  [Ref.  11]: 


J_ 
ax 


T,{x)-T{r*) 

m-Tm(x) 


=  0 


(2.2) 


where  T,  is  the  tube  surface  temperature,  T  is  the  local  fluid  temperature,  and  TB 
is  the  mean  temperature  of  the  fluid  over  the  cross  section  of  the  tube.  The  exact 
shape  of  the  fully  developed  temperature  profile  is  dependent  upon  the  boundary 
conditions  imposed,  either  a  uniform  surface  temperature  or  constant  heat  flux. 

For  laminar  flow,  the  thermal  entry  length  may  be  expressed  as  the 
following  [Ref.  10]: 


^f<u 


Q.Q5ReDPr 


(2.3) 


It  should  be  noted  that  the  rate  of  development  of  the  velocity  and 
temperature  profiles  in  the  entry  region  depend  upon  the  fluid  Prandtl  number.  For 
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Pr  =  1,  both  velocity  and  temperature  profiles  develop  at  the  same  rate.  However, 
for  Pr  >  1,  the  hydrodynamic  boundary  layer  develops  more  rapidly  than  the  thermal 
boundary  layer  (Xfdh  <  Xf4t).  For  very  large  Prandtl  number  (Pr  >  100),  Xf4h  is  so 
much  smaller  than  Xf4t  that  it  is  reasonable  to  assume  that  a  fully  developed  velocity 
profile  exists  throughout  the  thermal  entry  region.  [Refs.  9,  10,  11,  12] 

2.      Fully  Developed  Flow 

Fully  developed  laminar  flow  of  an  incompressible  fluid  through  a  circular 
duct  exists  when  both  the  thermal  and  velocity  profiles  are  fully  developed.  Two 
boundary  conditions  must  now  be  addressed,  that  of  a  uniform  wall  temperature  and 
constant  wall  heat  flux. 

For  these  conditions,  the  general  form  of  the  energy  equation  in  cylindrical 
coordinates  for  an  incompressible  fluid  having  constant  properties  (with  no  heat 
generation  and  neglecting  body  forces,  pressure  gradients  and  viscous  dissipation), 
is  given  by  [Ref.  11]: 


dx        dr      r  dr 


frill  (2.4) 


It  can  be  shown  that  for  these  two  cases,  the  energy  equation  reduces  to 
the  following  differential  equations, 


Constant  Wall  Heat  Flux:    J_  _<3_  (rf[\       _" 

r  dr       dr  a 
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Uniform  Wall  Temperature:    }_^__  (r—\  =!L 

r  dr       dr        a 
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(2.6) 


For  the  constant  wall  heat  flux  boundary  condition,  Kays  [Ref.  11]  solves 
the  differential  equation  (equation  2.5)  for  a  fully  developed  velocity  and 
temperature  profile  and  expresses  the  Nusselt  number  as: 


Nu  =  —  =  4.364 
k 


(2.7) 


Note  that  for  this  specified  case,  the  heat  transfer  coefficient,  h,  depends 
only  on  the  thermal  conductivity,  k,  and  the  tube  diameter,  D. 

A  solution  to  the  uniform  wall  temperature  case  is  more  complex,  and 
involves  evaluating  the  infinite  series  solution  to  the  differential  equation  (equation 
2.6)  for  the  hydrodynamically  developed,  thermally  developing  flow  situation.  In 
evaluating  the  series  solution,  the  resulting  asymptotic  Nusselt  number  is  given  by 
[Ref.  11]: 


Nu  =  3.657 


(2.8) 


This  result  is  16%  less  than  the  solution  for  the  constant  wall  heat  flux  boundary 
condition.  The  asymptotic  values  in  equations  2.7  and  2.8  are  shown  in  Figure  2.3. 
Note  that  the  inverse  of  the  Graetz  number  is  equivalent  to  the  dimensionless  axial 
distance  x\ 
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3.      Entry  Length  Problem 

The  solution  to  the  energy  equation  for  the  entry  region  is  far  more 
difficult  to  obtain  because  velocity  and  temperature  vary  both  axially  and  radially. 
There  are  two  entry  length  problems  to  be  considered,  a  combined  entry  length, 
where  both  temperature  and  velocity  profiles  develop  simultaneously,  and  the 
thermal  entry  length  problem,  where  the  thermal  conditions  develop  in  the  presence 
of  a  fully  developed  velocity  profile.  This  latter  case  would  exist  if  the  place  at  which 
heat  transfer  begins  is  preceded  by  an  unheated  starting  length,  a  situation  that  exists 
in  most  experimental  facilities.  It  should  be  noted  that  for  Pr  <<  1  (i.e.,  liquid 
metals),  the  thermal  boundary  layer  would  develop  faster  than  the  hydrodynamic 
boundary  layer.  However,  this  case  has  little  practical  application  and  is  therefore 
rarely  considered  in  the  literature. 

Solutions  have  been  obtained  for  both  entry  length  conditions  mentioned 
above  [Refs.  9,  10,  11,  12]  and  these  results  are  also  presented  in  Figure  2.3.  The 
Nusselt  number  is  plotted  against  a  dimensionless  axial  coordinate  x"  (the  inverse  of 
the  Graetz  number)  given  by: 

x'  =  -f^L  (2.9) 

Refr 

Fully  developed  conditions  are  reached  at  x*  «  0.05  where  the  solutions  for  the 
thermal  entry  length  converge  with  those  for  the  combined  entry  length.  For  x"  > 
0.05,  the  two  types  of  entry  length  asymptotically  reach  the  fully  developed  solutions 
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of  4.36  and  3.66  for  the  constant  wall  heat  flux  and  uniform  wall  temperature 
boundary  conditions  respectively. 
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Figure  2.3     Local  Nusselt  Number  Obtained  from  Entry  Length  Solutions  for 
Laminar  Flow  in  a  Circular  Tube  (Courtesy  [Ref.  10]) 

a.      Thermal  Entry  Region  -  Uniform  Wall  Temperature 

The  thermal  entry  length  problem,  first  studied  by  Graetz  in  1883  and 
later  by  Nusselt  in  1910,  is  more  commonly  known  as  the  Graetz  Problem.  They 
considered  an  incompressible  fluid  with  constant  properties  flowing  through  a 
circular  tube  having  a  fully  developed  laminar  velocity  profile  and  a  developing 
laminar  temperature  profile  subjected  to  a  uniform  wall  temperature.  Shah  and 
London  [Ref.  9]  present  a  numerically  derived  infinite  series  solution  to  the  Graetz 
Problem. 
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A  more  useful  correlation  presented  by  Incropera  and  Dewitt  [Ret 
10]  for  the  uniform  wall  temperature  thermal  entry  length  problem  was  developed 
by  Hausen  and  is  of  the  form: 

0.0668  (DIL)ReDPr 

Nu     =  3.66  + -J 1— £ (2. 10) 

U0.Q4[(D/L)ReDPr]2/3 

This  correlation  was  used  by  Mazzone  [Ref.  6]  in  his  refrigerant  condensation 
experiments  to  predict  the  average  inside  heat  transfer  coefficient  for  a  smooth  tube 
without  any  type  of  insert 

Shah  and  London  [Ref.  9]  present  a  more  complete  set  of  correlations 
derived  from  experimental  data  for  heat  exchanger  tubes  in  laminar  flow  for  local 
(NuJ  and  mean  (Num)  Nusselt  numbers  based  on  various  x"  ranges.  The  correlations 
for  local  Nusselt  number  are  given  by: 


Nux  =  l.Qllix'Y^-Q.l  forx'  <  0.01 

Nux  =  3.657  +6.874(103A:-)-0-488e-37-2r'  forx'  >  0.01 


(2.11) 


and  for  mean  Nusselt  numbers  by: 


Num  =  l.6l5(x'yu3-0.1  forx'  <Z  0.005 

M/m  =  1.615(jc')-1/3~0.2  for  0.005<*'  <  0.03 

0  0499 

Nu    =  3.657  +U-U4*y  for  x'>  0.03 


(2.12) 


The  reader  is  reminded  that  both  the  Hausen  correlation  and  the 
Shah  and  London  correlations  presented  above  are  for  a  thermally  developing  flow 
with  a  fully  developed  velocity  profile  subjected  to  a  uniform  wall  temperature. 
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b.  Thermal  Entry  Region  -  Constant  Wall  Heat  Flux 

For  the  thermal  entry  length  problem  under  a  constant  wall  heat  flux 
boundary  condition,  Shah  and  London  [Ref.  9]  present  several  numerical  solutions 
but  recommend  the  following  approximate  equations  for  the  local  Nusselt  number: 

Nux  =  1.302(^-)-^-l  forx*  <;  5xl0"5 

Nux  =  1302(x')-y3-05  for  SxW5  <;  x*  <;  1.5 xlO"3  (2.13) 

Nux  =  4.364+8.68(10V)-°J06<?-4U'  for  x*  >  1.5xl0-3 

Shah  and  London  also  recommend  the  following  approximate  equations  for  the  mean 
Nusselt  number: 

Num  =  1.953(r')-w  for  x'  <  0.03 

Nu     =  4.364  + for  x    >  0.03 

x 

The  above  two  sets  of  equations  are  for  a  thermally  developing, 

hydrodynamically  developed  flow  subjected  to  a  constant  wall  heat  flux. 

c.  Combined  Entry  Length 

As  discussed  earlier,  the  rate  at  which  the  temperature  and  velocity 
profiles  develop  is  dependent  upon  the  Prandtl  number.  At  medium  to  high  Prandtl 
numbers  (Pr  >  100),  the  velocity  profile  develops  much  more  rapidly  than  the 
temperature  profile  and  a  fully  developed  velocity  profile  assumption  is  valid;  for 
such  conditions  the  equations  given  above  should  be  used. 

For  1  <  Pr  <  100,  however,  a  fully  developed  velocity  profile  is  no 
longer  a  valid  assumption  and  a  combined  entry  length  solution  must  be  developed. 
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Again,  Shah  and  London  [Ref.  9]  present  a  comprehensive  study  of  various  solutions 
to  the  combined  entry  length  problem,  including  two  solutions  for  local  Nusselt 
numbers  developed  by  Churchill  and  Ozoe  [Ref  13],  the  first,  for  a  uniform  wall 
temperature: 


Nux+1.7 


5.357U+[(388/<n-)x,]-*9P/8 

tt/(284c*) 

k  (1  +(Pr/0.0468)2/3}1/2{l  +[(388/tt)*  '^P'4  j 


1  3/8 


4/3 


(2.15) 


and  the  second,  for  a  constant  wall  heat  flux: 


Nux+l 


1  + 


5.364(1 +[(220/7T)r-]-10/9P/l0 

7r/(115.2Q 

{ (1  +(Pr/0.0207)2/3}1/2{l  +[(220/tt>:  -]-l^P5  J 


1  3/10 


x5/3 


(2.16) 


From  the  preceding  paragraphs,  it  can  be  seen  that  the  subject  of 
laminar  flow  forced  convection  in  circular  ducts  has  been  studied  in  great  detail 
yielding  both  analytical  and  experimental  results.  The  reader  is  referred  to  Shah  and 
London  [Ref.  9]  and  Kays  and  Crawford  [Ref.  11]  for  a  more  complete  and  thorough 
review  of  this  subject. 
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C.      AUGMENTATION    TECHNIQUES    TO    ENHANCE    TUBESIDE    HEAT 
TRANSFER 

1.      General  Introduction 

As  discussed  earlier,  when  laminar  flow  heat  transfer  occurs  within  tubes, 
it  usually  represents  the  dominant  thermal  resistance  in  a  heat  exchanger.  These 
large  resistances  result  in  low  inside  heat  transfer  coefficients.  In  recent  years,  the 
requirement  for  more  efficient  heat  exchangers  has  stimulated  considerable  interest 
in  techniques  to  augment  or  enhance  tubeside  heat  transfer. 

The  various  methods  to  enhance  tubeside  heat  transfer  can  be  classified 
into  two  categories,  passive  and  active  techniques.  Passive  techniques  involve  no 
external  energy  (except  increased  pump  work)  to  produce  the  necessary 
enhancement  Passive  techniques  are  divided  into  the  following  methods: 


1.  Internally  finned  tubes  (longitudinal  or  spiral  fins)  designed  to  increase  the 
inside  heat  transfer  surface  area. 

2.  Surface  roughness  techniques  aimed  at  agitating  flow  rather  than  increasing 
heat  transfer  surface  area. 

3.  Swirl  flow  techniques  which  involve  the  use  of  inlet  vortex  generators, 
periodically  spread  propellers  or  twisted  tape  inserts  designed  to  cause 
enhancement  by  creating  a  rotating  and/or  secondary  flow. 

4.  Displaced  promoters  which  alter  the  fluid  flow  near  the  surface  by  the  use  of 
inserts  such  as  wire  mesh,  static  mixer  elements,  rings  or  disks. 

5.  Coolant  additives  designed  to  increase  the  conductivity  of  the  coolant. 

6.  Compound  techniques  that  involve  more  than  one  of  the  above  mentioned 
methods. 
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Active  techniques  utilize  an  external  energy  source  to  promote  the 
required  augmentation.  These  techniques  include:  (1)  mechanical  aids,  (2)  heated 
surface  vibration,  (3)  fluid  pulsation,  (4)  electrostatic  fields,  and  (5)  suction  or 
injection. 

All  of  the  above  techniques  are  aimed  at  creating  turbulence  in  the 
relatively  slow  moving  or  stagnant  boundary  layer  fluid  by  displacing  this  fluid  from 
the  heated  surface  and  mixing  it  with  the  bulk  fluid.  Essential  in  all  of  these 
techniques  is  a  balance  between  increased  heat  transfer  rate  and  increased  pumping 
power  requirements/costs.  Bergles  and  Joshi  [Ref.  14]  present  an  excellent  review 
of  augmentation  techniques  and  compare  heat  transfer  data  taken  at  low  Reynolds 
number  flow  subjected  to  both  unifornrwalf  temperature  and  constant  wall  heat  flux. 

Of  particular  interest  to  this  thesis  are  two  passive  techniques,  the  swirl 
flow  technique  using  twisted  tape  inserts  and  the  displaced  promoter  technique  using 
wire  mesh  inserts.  These  two  methods  of  enhancement  of  tubeside  heat  transfer  are 
discussed  in  more  detail  in  the  following  sections. 

2.      Heat  Transfer  Enhancement  with  Twisted  Tape  Elements 

Among  the  most  common  swirl  flow  augmentation  techniques  is  the  use 
of  twisted  tape  inserts.  The  tapes  usually  consist  of  a  thin  stainless  steel,  brass  or 
copper  strip  which  is  the  width  of  the  inside  diameter  of  the  tube.  The  tapes  are 
twisted  by  clamping  one  end  to  an  overhead,  attaching  weights  to  the  other  end  and 
then  twisting  to  the  desired  specifications  [Ref.  15].  Enhancement  occurs  not  only 
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due  to  the  increased  path  length  of  flow,  but  also  due  to  secondary  flow  effects  and 
fin  effects  caused  by  tape  contact  with  the  inside  tube  wall. 

Analytical  studies  of  twisted  tape  inserts  by  Date  and  Singham  [Ref.  16] 
suggest  that  swirl  flows  improve  heat  transfer  by  as  much  as  a  factor  of  70;  however, 
no  experimental  data  have  confirmed  their  unrealistically  large  predictions.  The  first 
experimental  work  of  swirl  flow  using  twisted  tapes  was  conducted  by  Hong  and 
Bergles  [Ref.  15].  Their  primary  objective  was  to  develop  experimentally  based 
correlations  for  predicting  heat  transfer  coefficients  for  laminar  flow  of  water  and 
ethylene  glycol  in  an  electrically  heated  tube  (heated  length  of  1.22  m,  inside 
diameter  of  10.2  mm)  with  two  sizes  of  twisted  tape  inserts.  They  showed  that  the 
Nusselt  number  was  independent  of  axial  location,  .suggesting  that  fully  developed 
conditions  existed.  Therefore,  the  mean  Nusselt  number  should  take  the  form: 

Num  =/(*e„iV,v).  (2.17) 

Re,  is  the  Reynolds  number  for  the  swirl  flow  given  by: 

Re    =  . 1^ (2.18) 

iru.(D(-4S) 

and  8  is  the  tape  thickness.  The  twist  ratio,  y,  is  defined  as: 

y  -.  !L  (2.19) 
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where  H  is  the  pitch  for  a  180°  twist  of  the  tape.  Hong  and  Bergles  then  combined 
Res  and  y  as  a  single  parameter  and  came  up  with  the  following  correlation  for 
values  of  Re,/y  >  10: 

Num  =  0.383Pr0J5(^^)a622  (2.20) 

For  values  of  Rejy  <  10,  the  following  correlation  was  obtained: 

Num  =  5.172(1  +5.484  x  lO^Pr^ReJy)125)^  (2.21) 

It  should  be  noted  that  all  fluid  properties  in  equations  2.20  and  2.21 
should  be  evaluated  at  the  bulk  fluid  temperature.  One  additional  important 
observation  was  made  by  Hong  and  Bergles  regarding  tube  wall  temperature.  They 
found  a  significant  variation  of  temperature  at  one  axial  location  with  two  different 
tape  orientations  at  that  position  of  90°  and  120°,  suggesting  that  the  temperature 
profile  was  related  to  the  tape  orientation.  Hong  and  Bergles  concluded  that  for  a 
fully  developed  flow,  inside  heat  transfer  coefficients  depend  on  twist  ratio  and  could 
be  improved  by  a  factor  of  two  or  three  over  an  empty  tube.  These  conclusions  were 
confirmed  by  later  studies  of  Saha  et.  al.  [Ref.  17]  and  Marner  and  Bergles  [Ref. 
18]. 

3.      Heat  Transfer  Enhancement  with  Wire  Mesh  (HEATEX)  Elements 

The  use  of  wire  mesh  inserts  to  promote  inside  heat  transfer  falls  into  the 
displaced  promoter  technique  mentioned  earlier.  One  such  wire  mesh  insert  was 
commercially  developed  by  CAL-GAVIN  Limited  (UK)  and  called  "HEATEX." 
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This  HEATEX  insert  consists  of  a  central  wire  core  on  which  a  series  of 
wire  loops  or  petals  are  attached.  The  petals  are  attached  to  the  core  at  an  inclined 
angle  that  face  the  oncoming  fluid  flow.  Ideally,  each  petal  makes  contact  with  the 
tube  wall  in  the  form  of  an  arc,  thereby  maximizing  disturbance  of  fluid  along  the 
tube  wall.  In  addition  to  the  boundary  layer  disturbance,  radial  mixing  occurs  caused 
by  repeated  redistribution  and  mixing  of  fluid  through  the  mesh  of  loops  near  the 
tube  axis  [Ref.  19].  Various  wire  loop  densities  and  nominal  diameters  are 
available,  depending  on  the  application. 

Initial  studies  by  Gough  et  al.  [Ref.  19]  with  HEATEX  inserts  suggest 
enhancements  of  a  factor  greater  than  ten  times  that  of  a  similar  tube  with  no  insert 
Additional  experiments  were  conducted  by  Oliver  and  Aldington  [Ref.  20]  in  which 
a  correlation  was  determined: 

Nu  =  0.232Re0MPr0A6  (2-22) 

It  must  be  remembered,  however,  that  this  correlation  is  for  a  specific  loop  density 
and  nominal  diameter  HEATEX  element. 

The  two  types  of  HEATEX  mixing  elements  supplied  by  CAL-GAVIN 
(nominal  diameters  of  13.26  and  10.16  mm)  did  not  match  the  loop  densities  of 
previous  studies;  therefore  new  correlations  needed  to  be  developed.  Mazzone  [Ref. 
6]  developed  correlations  for  these  inserts  from  limited  data  supplied  by  CAL- 
GAVIN.  For  the  13.26  mm  and  10.16  mm  elements,  these  correlations  were 
respectively: 
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Nu  =  0.226te°-65Pr0-46     (13.26mm)  (2-23) 

and 

Nu  =  0.063fo°-76/V0-46     (10.16mm)  (2-24) 

Since  these  correlations  were  developed  from  a  very  limited  amount  of 
data  taken  under  conditions  that  did  not  exactly  match  that  of  Mazzone's 
condensation  experiments,  they  were  always  subject  to  inaccuracy.  Therefore,  the 
need  for  more  accurate  correlations  under  the  exact  conditions  used  in  the 
condensation  experiments  are  required. 
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III.  EXPERIMENTAL  APPARATUS 

A.      SYSTEM  OVERVIEW 

The  test  apparatus  utilized  is  shown  in  Figures  3.1  and  3.2.  The  apparatus  is  a 
closed  loop,  low  pressure  system  utilizing  an  ethylene  glycol/water  mixture  as  the 
coolant  The  coolant  is  circulated  by  a  centrifugal  pump  from  a  supply  tank,  through 
a  reducing  manifold  and  flowmeter  into  the  test  section.  Coolant  is  then  returned  to 
the  supply  tank  via  a  return  header. 

The  test  section  consisted  of  an  inlet  and  outlet  thermocouple  well/mixing 
chamber  to  measure  inlet  and  outlet  coolant  temperatures  and  the  test  tube  itself. 
The  tube  was  instrumented  with  twelve  thermocouples  (embedded  in  the  tube  wall) 
to  determine  inside  wall  temperatures.  A  constant  heat  flux  boundary  condition  was 
achieved  using  a  series  of  flexible  electrical  heating  elements  wrapped  around  the 
tube.  Monitoring  the  system  as  well  as  data  acquisition  and  data  processing  were 
provided  by  a  Hewlett-Packard  microcomputer  and  Data  Acquisition  System. 
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Figure  3.1   Primary  and  Secondary  Coolant  Systems. 
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Figure  3.2   Data  Acquisition  and  Control  Systems. 
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B.   TUBES  AND  INSERTS  TESTED 

1.      Tubes  Tested 

Three  types  of  tubes  were  utilized  in  the  experiments,  a  copper  smooth 
tube,  a  copper/nickel  externally  finned  tube,  and  a  copper/nickel  corrugated  or  roped 
tube.  These  tubes  correspond  to  the  tubes  utilized  in  the  refrigerant  condensation 
experiments  of  Mabrey  [Ref.  5]  and  Mazzone  [Ref.  6]. 

The  smooth  copper  tube  provided  a  baseline  set  of  experimental  data 
which  were  then  compared  with  previous  experimental  and  theoretical  work  for  the 
inside  heat  transfer  coefficient  conducted  by  Shah  and  London  [Ref.  9]  for  no  insert, 
Oliver  and  Aldington  [Ref.  20]  for  the  HEATEX  insert,  and  Hong  and  Bergles  [Ref. 
15]  for  the  twisted  tape  insert 

The  Cu/Ni  finned  tube  provided  a  smaller  inside  diameter  to  be  used  in 
these  comparisons.  Note  that  the  fins  play  little  part  in  the  experiments.  The  wall 
thermocouples  were  buried  in  the  tube  wall  and  were  used  to  measure  the  inside  wall 
temperature,  which  had  a  smooth  bore.  Consequently,  for  the  purposes  of  these 
experiments,  the  tube  behaved  like  a  smooth  copper  tube  with  a  different  internal 
diameter. 

The  final  tube  tested  was  a  Cu/Ni  corrugated  (or  roped)  tube, 
commercially  referred  to  as  a  "Korodense"  tube.  The  tube  was  manufactured  and 
provided  by  Wolverine  Tube  Company  of  Decatur,  Al.  All  tube  dimensions  and 
specifications  were  taken  from  the  Wolverine  Technical  Bulletin  No.  4020  [Ref. 
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21].  Figure  3.3  is  a  schematic  of  the  Korodense  tube.  Tube  specifications  for  all 
three  tubes  are  listed  in  Table  E.1  of  Appendix  E. 
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Figure  3.3   Wolverine  Korodense  Tube  (Courtesy  [Ref.  21]) 

2.      Insert  Elements  Tested 

Two  types  of  insert  mixing  elements  were  utilized  in  this  investigation,  a 
twisted  tape  insert  and  a  wire  mesh  insert  commercially  referred  to  as  a  "HEATEX" 
tube  insert.  In  both  cases,  two  sizes  were  required  due  to  the  two  different  tube 
inside  diameters  used. 

The  twisted  tape  inserts  were  made  of  a  strip  of  brass  0.559  mm  thick  and 
either  10.16  mm  or  13.26  mm  wide.  The  tapes  were  manufactured  in  the  NPS 
Mechanical  Engineering  machine  shop  by  clamping  them  to  the  overhead,  attaching 
a  weight  to  the  opposite  end,  and  twisting  them  to  the  desired  twist  ratio.  A  twist 
ratio  of  3  and  4  corresponding  to  the  two  respective  tube  inside  diameters  were 
manufactured.  Once  installed  in  the  tube,  the  edges  of  the  tape  should  ideally  have 
been  in  contact  with  the  tube  wall. 
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The  second  type  of  insert  used  was  a  commercially  available  product 
referred  to  as  HEATEX  radial  mixing  elements.  These  elements  were  supplied  by 
Cal  Gavin  of  Birmingham,  U.K.  They  are  manufactured  from  stainless  steel  and 
consist  of  a  central  wire  core  onto  which  a  series  of  wire  loops  or  petals  are  attached. 
Each  petal  was  inclined  at  an  angle  facing  the  oncoming  flow.  The  outside  diameter 
of  the  elements  were  slightly  oversized  for  the  tube,  ensuring  that  close  contact  with 
the  tube  wall  in  the  form  of  an  arc,  rather  than  a  point,  was  achieved.  As  with  the 
twisted  tape  elements,  different  diameter  HEATEX  elements  were  needed,  13.26 
mm  used  in  the  smooth  and  Korodense  tubes  and  the  other  10.15  mm  for  the  Cu/Ni 
finned  tube.  Figure  3.4  shows  both  the  twisted  tape  and  HEATEX  inserts. 


Figure  3.4   Twisted  Tape  and  HEATEX  Inserts 
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C      COOLANT  SYSTEM 

1.  Primary  Coolant  System 

An  ethylene  glycol/water  mixture  (54/46%  by  weight)  was  stored  in  a  1.5 
m3  insulated  supply  tank.  The  mixture  was  circulated  through  the  primary  system  by 
a  l¥t  hp,  constant  speed,  centrifugal  pump.  The  coolant  was  discharged  from  the 
pump  via  a  76  mm  diameter  PVC  pipe  to  a  manifold  where  it  was  reduced  to  a  15.9 
mm  diameter  Tygon  flexible  tube.  A  rotameter  type  flowmeter  with  a  ball  valve  at 
its  entrance  was  used  to  control  the  flow  rate  through  the  system.  Flowmeter 
calibration  results  are  discussed  in  Appendix  A. 

Flow  entered  the  test  section  through  a  90°  bend  which  housed  the 
thermocouple  well  and  exited  through  a  mixing  chambex/the_rmpcouple  well 
combination  (see  Figure  3.1).  Prior  to  entering  the  exit  mixing  chamber,  the  coolant 
passed  through  a  10  cm  length  of  flexible  tubing  containing  a  high  density  wire  mesh. 
This  wire  mesh,  in  conjunction  with  the  mixing  chamber,  served  to  minimize  radial 
temperature  variations  and  provide  a  uniform  outlet  coolant  temperature.  The  flow 
then  exited  the  test  section  through  more  flexible  tubing  to  a  manifold  and  return 
PVC  header,  suspended  above  the  apparatus,  and  back  to  the  coolant  supply  tank. 

2.  Secondary  Coolant  System 

In  order  to  achieve  a  wide  range  of  Reynolds  number,  coolant  inlet 
temperatures  were  varied  from  -10°C  to  +20°C.  This  was  accomplished  by  a 
secondary  coolant  system  comprised  of  a  28  kW  (8  ton)  external  refrigeration  system. 
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A  30  gpm,  .75  hp  centrifugal  pump  continuously  circulated  coolant  between  the 
supply  tank  and  a  chiller  barrel  where  it  was  cooled  to  the  desired  temperature. 
Temperature  was  maintained  by  a  thermostatically  controlled  solenoid  valve  at  the 
discharge  side  of  the  chiller  barrel.  This  system  allowed  for  a  range  of  coolant 
temperatures  between  -20°C  and  ambient  to  be  selected  while  maintaining  the 
coolant  inlet  temperature  approximately  constant  (±2°C)  during  each  run. 

D.      INSTRUMENTATION  AND  TEST  SECTION 

1.      Instrumentation 

All  temperature  measurements  were  made  using  copper/constantan  type-T 
teflon  coated  thermocouples  and  were  previously  calibrated  by  Mabrey  [Ref.  5].  The 
inlet  coolant  temperature  was  measured  by  a  single  thermocouple  while  the  outlet 
coolant  temperature  was  measured  by  two  thermocouples  located  immediately  after 
the  mixing  chamber. 

Each  test  tube  was  fitted  with  twelve  wall  thermocouples  to  determine  the 
average  inside  wall  temperature.  Thermocouples  were  placed  at  four  different 
longitudinal  positions  along  the  axis  of  the  tube.  Three  thermocouples  at  each 
location  were  spread  evenly  around  the  tube  at  120°C  intervals.  Figure  3.5  shows  a 
typical  tube  with  the  exact  longitudinal  and  nodal  thermocouple  locations.  Actual 
longitudinal  locations  varied  slightly  for  each  tube  and  are  given  in  Table  E.l  as  x/D; 
where  x  is  the  length  from  the  start  of  the  heated  length  to  the  thermocouple 
location  and  D;  in  the  inside  diameter  of  the  tube.  The  two  extreme  longitudinal 
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locations  of  the  thermocouples  corresponded  to  the  beginning  and  end  of  the 
"effective"  condensing  length  of  the  tubes  being  used  in  the  condensation  experiments 
of  Mabrey  [Ref.  5]  and  Mazzone  [Ref.  6]. 


INSTRUMENTED  TUBE 
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Figure  3.5    Instrumented  Tube 
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2.      Thermocouple  Attachment  Method 

In  order  to  attach  the  thermocouples  to  the  wall  of  the  tube,  short 
longitudinal  grooves  were  cut  at  each  location.  These  longitudinal  grooves  were  cut 
to  a  depth  of  xh  the  tube  wall  thickness.  It  was  decided  that  a  thermocouple 
embedded  within  the  tube  wall  vice  being  attached  to  the  surface  would  lead  to  a 
more  accurate  tube  wall  temperature  and  minimize  the  effect  of  direct  contact 
heating  from  the  heating  elements.  Thermocouples  were  then  placed  in  tnese  grooves 
and  either  soldered  in  place  or  tack  welded  and  covered  with  a  liquid  metal  epoxy. 
Different  methods  for  attaching  and  covering  the  thermocouples  were  used  for  each 
tube  based  on  the  tube  material  and  configuration  (discussed  below). 

In  the  case  of  the  smooth  copper  tube,  the  thermocouples  were  silver 
soldered  in  place  and  the  remainder  of  the  groove  then  covered  with  silver  solder. 
Since  the  tubes  were  copper,  this  was  relatively  easy  to  accomplish.  Care  was  taken 
to  ensure  that  the  thermocouple  remained  in  contact  with  the  tube  wall  and  that  no 
air  pockets  were  created  (which  would  cause  an  additional  large  resistance  to  heat 
flow).  Also,  great  care  was  taken  to  prevent  any  melting  of  the  teflon  insulation  past 
the  thermocouple  joint.  This  condition  could  have  caused  an  electrical  contact 
between  the  two  wires  of  the  thermocouple  prior  to  the  tube  wall,  leading  to  an 
unrealistically  high  observed  wall  temperature. 

A  different  method  was  required  to  attach  the  thermocouples  to  the 
copper/nickel  finned  tube  due  not  only  to  the  external  fins  but  also  the  10%  nickel 
content.  Initially,  a  small  diameter  hole  was  machined  into  the  tube  to  minimize  fin 
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damage  (rather  than  cutting  a  longitudinal  groove).  Since  the  thermocouple  was  to 
be  embedded  half  way  into  the  tube  wall,  this  hole  proved  to  be  relatively  deep. 
Also,  because  of  the  nickel  content  of  these  tubes,  the  preheat  temperature  required 
to  ensure  that  the  silver  solder  would  adhere  to  the  tube  wall  was  much  higher. 
These  higher  required  temperatures  and  the  hole  depth  caused  a  significant  amount 
of  the  thermocouple  teflon  insulation  to  be  melted,  allowing  for  potential  contact 
between  the  thermocouple  wires  at  some  point  other  than  the  tube  wall.  This 
technique  was  therefore  abandoned  in  favor  of  the  longitudinal  groove  method. 
However,  thermocouple  insulation  was  again  being  destroyed  due  to  the  higher 
required  preheat  temperatures.  An  alternative  solder,  soft  solder,  with  lower  preheat 
and  melting  temperatures,  was  therefore  used.  The  final  combination  of  longitudinal 
grooves  and  soft  solder  was  successful. 

For  the  copper/nickel  Korodense  tube,  the  groove  and  soft  solder  method 
was  again  tried.  Initial  tests  in  the  apparatus  showed  wall  temperatures  at  one 
longitudinal  position  varying  by  as  much  as  15°C.  On  closer  examination  of  the 
thermocouple  connections  (using  a  magnifying  glass),  it  was  revealed  that  insulation 
damage  was  evident  on  over  half  of  the  thermocouples.  Several  more  attempts  at  the 
groove  and  soft  solder  method  were  tried,  all  leading  to  the  same  results. 

A  third  method  utilizing  a  tack  welding  machine  (acquired  and  adapted 
for  our  use)  was  tried.  Thermocouples  were  tack  welded  into  the  bottom  of  the 
groove  resulting  in  no  damage  to  the  insulation  and  a  liquid  steel  epoxy  was  used  to 
fill  the  remainder  of  the  groove,  providing  additional  strength  and  protection  to  the 
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connection.  This  final  method  proved  to  be  the  most  successful  and  is  recommended 
for  use  in  any  follow-up  experiments. 

3.      Test  Section 

In  order  to  achieve  a  constant  heat  flux,  each  tube  was  tightly  wrapped 
with  six  silicone  rubber  flexible  electrical  heating  elements  rated  at  240  V,  200  W 
each.  Each  strip  was  2.5  cm  wide  and  50.6  cm  long.  The  heating  elements  were 
connected  in  parallel  to  a  240  Vac  power  supply  via  a  plug  panel  and  power 
switchboard.  Power  to  the  six  heaters  was  controlled  by  a  STACO,  240  Vac,  23.5  kVa 
rheostat  (see  Figure  3.2). 

To  minimize  heat  loss  from  the  heater  to  ambient,  the  tube  was  wrapped 
in  a  20  mm  thick  neoprene  type  insulation.  This  type  of  insulation  proved  to  be 
extremely  effective  giving  a  negligible  heat  loss.  It  also  provided  insulation  against 
heat  influx  from  ambient  at  the  colder  coolant  inlet  temperatures.  Actual  heat  loss 
and  heat  influx  experiments  and  their  results  are  discussed  in  more  detail  in 
Appendices  B  and  C  All  other  tubing  and  PVC  piping  was  wrapped  in  similar  type 
insulation. 

The  whole  test  section  was  then  placed  in  a  restraining  stand  which 
minimized  outside  vibration  effects  and  ensured  that  each  tube  remained  horizontal 
and  in  the  same  relative  position  (to  the  remainder  of  the  apparatus)  for  each 
experiment. 
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E.       DATA  ACQUISITION  AND  CONTROL  SYSTEM 

A  Hewlett-Packard  9125  microcomputer  (MC)  was  used  to  control  a  Hewlett- 
Packard  3497A  Automatic  Data  Acquisition  Unit  (DAU)  (see  Figure  3.2).  The  DAU 
read  the  output  from  the  twelve  tube  wall  thermocouples  and  three  coolant 
thermocouples  (one  inlet  and  two  outlet).  Additional  channels  were  assigned  to  sense 
heater  input  voltage  and  current  at  the  plug  panel.  The  voltage  was  converted  to 
heater  power  using  the  heater  resistance  and  the  calibration  discussed  in  Appendix 
D. 

Thermocouple  readings  were  made  in  microvolts  which  were  converted  to 
temperatures  (along  with  the  heat  flux  calculations)  using  the  data  reduction 
program.  The  data  reduction  program,  DRPSING,  will  be  discussed  in  more  detail 
in  a  later  chapter.  All  temperatures,  power,  and  heat  flux  measurements  were 
monitored  using  the  MC  and  DAU.  DAU  channel  assignments  are  listed  in  Table 
3.1. 
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TABLE  3.1   DATA  ACQUISITION  UNIT  CHANNEL  ASSIGNMENTS 


CHANNEL  NO. 

MEASURED  QUANTITY 

1 
2 
3 
4 
5 
6 
7 
•     8 

Thermocouple  Pos  1 
Thermocouple  Pos  2 
Thermocouple  Pos  3 
Thermocouple  Pos  4 
Thermocouple  Pos  5 
Thermocouple  Pos  6 
Thermocouple  Pos  7 
Thermocouple  Pos  8 

9 
10 
11 
12 
13 

Thermocouple  Pos  9 
Thermocouple  Pos  10 
Thermocouple  Pos  11 
Thermocouple  Pos  12 
Tube  Inlet 

14,  15 

Tube  Outlet 

25 
29 

Input  Current 
Input  Voltage 
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IV.   EXPERIMENTAL  PROCEDURES 

A.      CALIBRATION  EXPERIMENTS 

Several  calibration  experiments  were  required  prior  to  performing  the  main 
experiments.  Figure  4.1  depicts  a  typical  tube  being  tested.  The  primary  heat  source 
is  from  the  heating  elements  (qUp«).  However,  to  complete  the  heat  balance  on  the 
system,  several  other  heat  sources  and  heat  sinks  needed  to  be  accounted  for.  The 
first  experiment  conducted  was  to  determine  the  amount  of  heat  being  lost  from  the 
heating  elements  to  the  atmosphere  vice  going  to  the  coolant  Heat  loss  experiments 


are  discussed  in  Appendix  B. 

The  second  experiment  was  to  account  for  two  additional  potential  heat 
sources:  heat  leakage  into  the  coolant  from  ambient  (qieak)  and  frictional  temperature 
rise  (qfnct).  Approximately  8  cm  at  either  end  of  the  tube,  along  with  the  flexible 
tubing  section  prior  to  the  thermocouple  wells,  were  not  covered  by  the  heating 
elements.  Since  coolant  inlet  temperatures  were  primarily  below  ambient 
temperature,  a  potential  for  heat  leakage  to  the  coolant  from  ambient  existed.  Also, 
because  of  these  cold  inlet  temperatures,  frictional  heating  could  exist,  especially 
when  using  the  various  inserts.  Heat  leakage  into  the  system  and  frictional  effects 
were  considered  together  as  one  heat  influx.  These  experiments  were  conducted  with 
the  heater  tapes  turned  off;  any  temperature  rise  was  then  monitored  using  the  inlet 
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and  outlet  thermocouples.  Appendix  C  discusses  the  experiments  conducted  to 
determine  this  part  of  the  heat  balance  equation.  As  mentioned  earlier,  the  test 
section  and  all  other  tubing  were  insulated  to  minimize  these  effects. 

HEAT  BALANCE 

q  *—       q 

P-  X  i  =  \  it 


INLET 


OUTLET 


•  • 


^t«*  q-.*  q^-*  <*«-"  *.«. 

rV*  tap*  *••*  inn  ■**• 

Figure  4.1   Heat  Balance. 

B.       PROCEDURES  FOR  DETERMINATION  OF  INSIDE  HEAT  TRANSFER 
RATE 

1.      General  Procedures 

Each  tube  utilized  in  this  study  was  tested  with  no  insert  and  the 
appropriate  twisted  tape  and  HEATEX  inserts  corresponding  to  the  tube  inside 
diameter.  For  each  insert  condition,  four  nominal  coolant  inlet  temperatures,  -10°C, 
0°C,  +  10°C  and  +20°C  were  chosen  to  increase  the  range  of  experimental  Reynolds 
and  Prandlt  numbers  that  could  be  attained.  To  check  that  heat  flux  variations  did 
not  affect  the  results,  data  were  taken  at  heat  fluxes  of  10  kW/m2  and  15  kW/m2  for 
each  inlet  temperature.  For  each  insert,  coolant  inlet  temperature  and  heat  flux 
setting,  coolant  velocity  was  varied  between  0.2  m/s  and  1.4  m/s.  Since  the  inside 
diameter  of  the  Cu/Ni  finned  tube  was  smaller  than  the  other  two  tubes,  higher 
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coolant  velocities  of  up  to  2.4  m/s  could  be  attained.  A  matrix  of  the  individual  data 
runs  and  data  file  names  is  shown  in  Table  4.1. 

With  the  tube  installed  and  all  calibration  experiments  completed,  the 
experiments  to  determine  the  inside  heat  transfer  coefficient  were  conducted.  Since 
flowmeter  settings  for  specific  coolant  velocities  were  based  upon  the  inlet  coolant 
temperature,  the  system  was  initially  operated  at  a  100%  flowmeter  setting  with  no 
power  to  the  heating  elements  for  15  minutes.  This  allowed  the  inlet  thermocouple 
well  to  stabilize.  Once  the  desired  inlet  coolant  temperature  was  reached  and 
maintained,  power  to  the  heating  elements  was  set  to  the  predetermined  heat  flux 
setting  for  that  data  run  (10  or  15  kW/m2).  The  flowmeter  was  then  adjusted  to 
achieve  the  desired  coolant  velocity  based  upon  flowmeter  calibration  results 
(discussed  in  Appendix  A).  Flow  was  maintained  for  approximately  10  minutes  to 
ensure  steady  state  equilibrium  conditions  prior  to  taking  data. 

TABLE  4.1   DATA  FILE  NAMES 


TEMP(C)/ 

Cu  Siooth  Tube 

Cu/Ni  Finned  Tube 

Cu/Nl  Korodense  Tube 

No 

Insert 

Twisted 
Tape 

HEATEX 

NO 

Insert 

Twisted 

Tape 

HEAT EX 

No 

Insert 

Twisted 
Tape 

HEAT EX 

20/10 

SUNI04 

• 

3UHX04 

CFNI04 

(  +  15) 

CFT01 

(+15) 

CFHX01 

(+15) 

KDNI01 

KDTT01 

K0HX01 

20/15- 

SUN I 04 

• 

SMHX04 

CFNI04 

r»is) 

CFTT01 
(+15) 

CFHX01 

(  +  15) 

KDNI01 

KDTT01 

KDHX01 

10/10 

SMNI03 

SMTT01 

SMHX01 

CFNI01 

CFTT02 

CFHX02 

K0NI02 

K0TT02 

K0HX02 

10/15 

SMNI03 

SMTT01 

SMHX01 

. 

CFTT02 

CFHX02 

KDNI02 

K0TT02 

KDHX02 

0/10 

SMNI02 

SMTT02 

3MHX02 

CFNI03 

CFTT03 

CFHX03 

K0NI03 

K0TT03 

K0HX03 

0/15 

SMNI02 

SMTT02 

SMHX02 

CFNI03 

CFTT03 

CFHX03 

K0NI03 

K0TT03 

K0HX03 

-10/10 

SMNI01 

SMTT03 

SMHX03 

CFNI02 

CFTT04 

CFHX04 

K0NI04 

K0TT04 

K0HX04 

•10/15 

SMNI01 

SMTT03 

3MHX03 

CFNI02 

CFTT04 

CFHX04 

K0NI04 

KDTT04 

K0HX04 
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Velocities  were  increased  from  0.2  m/s  to  1.4  m/s  in  steps  of  0.2  m/s  and 
then  decreased  from  1.3  m/s  to  0.3  m/s  in  similar  steps.  This  method  of  increasing 
on  even  values  and  decreasing  on  odd  values  was  done  not  only  to  determine  if  any 
hysteresis  effects  existed  (i.e.,  a  difference  in  the  increasing  and  decreasing  data),  but 
also  to  ultimately  increase  the  number  of  data  points  acquired.  Also,  since  data  were 
immediately  plotted  (as  Nun  versus  Re),  the  method  allowed  for  an  easy 
determination  as  to  whether  a  data  point  had  reached  a  steady  state  condition. 

Inlet  and  outlet  coolant  temperatures,  tube  wall  temperatures,  heat  flux 
and  heating  element  power  settings  were  all  continuously  monitored  using  the 
computer  program  DRPSING  in  conjunction  with  the  Data  Acquisition  System 
discussed  earlier  (DRPSING  is  listed  in  Appendix  G).  Individual  data  runs  for  each 
tube  and  insert  condition  were  conducted  in  a  similar  manner. 

2.      Specific  Experiments 

a.      Copper  Smooth  Tube 

The  first  tube  tested  was  the  copper  smooth  tube.  This  tube  was 
chosen  to  provide  a  baseline  set  of  data  needed  for  validation  of  system  performance 
and  for  comparison  with  data  acquired  using  other  tubes.  It  should  also  be  noted  that 
two  complete  independent  sets  of  data  covering  all  insert  conditions  and  inlet 
temperatures  were  taken  over  a  period  of  two  weeks.  Data  from  both  sets  were 
compared  and  found  to  be  in  very  good  agreement  (±2%);  this  validated  the 
repeatability  of  data  obtained  from  the  system  under  the  same  conditions. 
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The  initial  intention  for  all  tubes  was  to  start  with  a  given  insert  at 
the  highest  coolant  inlet  temperature,  take  data  at  each  of  the  heat  flux  settings  and 
then  proceed  to  the  next  coolant  temperature.  Once  one  insert  was  completed,  the 
system  could  then  be  shut  down  and  the  insert  changed.  This  procedure  would  then 
be  repeated  for  all  three  insert  conditions.  However,  the  coolant  supply  tank  and 
coolant  system  were  used  concurrently  for  other  experiments  which  required  the 
coolant  temperature  to  be  -20°C.  Since  it  took  approximately  36-48  hours  for  the 
coolant  temperature  to  increase  from  -10°C  to  +20°C  (the  sump  was  not  fitted  with 
a  heater),  only  the  copper  smooth  tube  was  done  in  this  way.  The  method  for  the 
other  two  tubes  is  discussed  below. 

b.      Cu/Ni  Externally  Finned  Tube 

The  second  tube  tested  was  the  Cu/Ni  finned  tube.  For  a  given  inlet 
temperature  and  tube  insert,  the  same  procedure  as  with  the  copper  smooth  tube  was 
used.  However,  rather  than  decreasing  coolant  temperature  with  the  same  insert,  the 
coolant  temperature  was  maintained  and  the  insert  changed  until  all  three  of  the 
insert  conditions  were  completed  for  a  given  inlet  temperature.  Once  all  the  data  for 
a  given  inlet  temperature  had  been  taken,  the  supply  tank  temperature  was 
decreased  to  the  next  inlet  temperature  condition.  It  took  only  about  one  hour  for 
the  secondary  coolant  system  to  decrease  the  temperature  by  10°C.  Since  it  was  a 
relatively  short  task  to  change  inserts,  the  36  to  48  hour  waiting  period  required  to 
raise  coolant  temperatures  back  to  +20°C  between  insert  changes  was  eliminated. 


43 


It  should  be  noted  that  for  the  Cu/Ni  finned  tube,  the  highest  nominal 
inlet  temperature  attained  was  only  +  15°C  vice  +20°C.  Colder  than  normal  ambient 
temperatures  existed  preventing  supply  tank  temperatures  from  reaching  +20°C. 
Also,  due  to  the  smaller  inside  diameter,  coolant  velocities  could  be  varied  between 
0.3  m/s  and  2.4  m/s. 

c.      Cu/Ni  Korodense  Tube 

The  same  procedure  used  with  the  finned  tube  was  also  used  with  the 
Korodense  tube.  In  addition,  another  test  was  conducted  with  this  tube.  It  was  noted 
that  larger  than  normal  wall  temperature  variations  existed  at  one  longitudinal 
location  when  using  the  two  inserts.  It  was  suspected  that  the  point  of  contact 
between  the  insert  and  tube  wall  was  directly  below  a  thermocouple  position  causing 
a  higher  wall  temperature  at  that  point  To  verify  this,  the  system  was  shut  down  by 
securing  the  coolant  pump  and  the  heating  elements  so  that  the  flowmeter  and 
power  settings  remained  exactly  the  same.  The  outlet  end  of  the  tube  was  opened 
and  the  insert  moved  downstream  by  approximately  2  cm.  The  system  was  then 
closed  and  the  coolant  flow  and  power  to  the  heating  elements  restored.  Wall 
temperature  variations  still  existed  but  were  significantly  less.  This  procedure  was 
repeated  several  times  with  similar  results. 
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V.   DATA  REDUCTION 

A.      AVERAGE  INSIDE  HEAT  TRANSFER  COEFFICIENT 

The  main  objective  of  the  data  reduction  scheme  was  to  determine  the  average 
inside  heat  transfer  coefficient  (hj).  This  in  turn  was  used  to  determine  a  mean 
Nusselt  number  (for  the  whole  heated  length  of  the  tube)  to  be  used  as  the  basis  for 
establishing  the  desired  correlations. 

To  accomplish  this  task,  the  Data  Acquisition  Unit  discussed  in  Chapter  III 
made  use  of  the  computer  program  DRPSING.  A  listing  of  the  program  is  given  in 
Appendix  E.  The  program  allowed  the  user  te  monitor  information  from  the 
instrumentation  attached  to  the  apparatus,  acquire  and  then  process  the  data  to 
determine  the  desired  output  In  doing  this,  the  user  would  select  the  tube  and  insert 
condition  as  well  as  the  nominal  coolant  inlet  temperature  for  that  particular  data 
run.  The  desired  coolant  velocity  was  then  selected.  Based  on  the  actual  measured 
coolant  inlet  temperature,  coolant  mixture  properties  and  flowmeter  calibration 
results  (see  Appendix  A),  the  appropriate  flowmeter  setting  (as  a  %)  was  calculated 
and  set  on  the  flowmeter. 

The  program  then  monitored  coolant  inlet  and  outlet  temperatures,  the  twelve 
wall  thermocouple  temperatures  and  the  actual  measured  heat  flux  and  power 
settings.  Once  all  readings  reached  a  steady  state  condition  (approximately  10  min.), 
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the  user  could  terminate  the  display  session  and  acquire  data  for  that  coolant  velocity 
setting. 

The  DAU  measures  emf  from  the  thermocouples  and  then  converts  these 
values  to  temperatures  using  the  temperature  calibration  equation  developed  by 
Mabrey  [Ref.  5].  Heater  power  was  determined  by  measuring  the  voltage  at  the  plug 
panel  and  calculating  the  power  based  on  heater  resistance  and  the  power  calibration 
contained  in  Appendix  D. 

The  twelve  wall  thermocouples  (placed  axially  and  radially)  were  used  to 
minimize  the  uncertainty  in  determining  the  average  inside  wall  temperature.  These 
twelve  wall  temperatures  were  averaged  and  corrected  (to  account  for  the  depth  of 
burial  using  a  simple  conduction  calculation)  to  determine  the  average  inside  wall 
temperature,  T^. 

The  physical  and  thermodynamic  properties  of  the  ethylene  glycol/water 
mixture  (density,  dynamic  viscosity,  thermal  conductivity,  heat  capacity)  were 
determined  using  equations  and  figures  from  Cragoe  [Ref.  22]  and  Gallant  [Ref. 
23]  and  correlated  by  Mabrey  [Ref.  5].  From  these  quantities,  the  coolant  mass 
flow  rate  was  calculated  by: 

m  =  VcPAt  (5.1) 

where  A;  is  the  inside  cross  sectional  area  of  the  tube  and  Vc  is  the  set  coolant 
velocity.  The  coolant  Reynolds  number,  ReD  was  calculated  by: 
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V  Dp 
ReD  =  -1^  (5-2) 

where  D;  is  the  inside  diameter  of  the  tube.  The  dimensionless  distance  (x*)  in  the 
flow  direction  is  calculated  as: 

xlD. 
x'  =  '  (5.3) 

ReDPr 

where  x  is  the  axial  distance  from  the  beginning  of  the  heated  length. 

To  determine  the  heat  transfer  rate  to  the  coolant,  several  factors  were 

considered.  The  basic  equation  to  determine  the  heat  transfer  rate  is  given  by: 

H  =  mc,A7  (5.4) 

where  AT  is  the  measured  coolant  temperature  rise  along  the  tube.  The  temperature 
rise  required  is  that  for  the  actual  heated  length  of  the  tube  that  is  subjected  to  the 
constant  heat  flux  condition,  ATHL.  To  determine  ATHL  from  AT,  corrections  have  to 
be  made  to  account  for  heat  loss,  heat  leakage  and  heat  transferred  to  the  unheated 
ends.  Figure  5.1  depicts  a  typical  tube  and  the  correction  terms  considered  in  the 
heat  balance. 

The  overall  AT  is  first  corrected  for  heat  loss  to  the  ambient,  (qio,,)  and  heat 
leakage  into  the  system,  (qleak).  The  heat  leakage  term  consisted  of  any  temperature 
rise  due  to  the  influx  of  heat  from  ambient  as  well  as  that  caused  by  frictional 
effects.  Appendix  D  discusses  the  experiments  conducted  to  determine  both  this  heat 
loss  and  heat  leakage  correction. 
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TEMPERATURE  RISE  CORRECTIONS 


t/c.in 


1 C 


IOS8 


\ 


ATHL 
AT. 


a  9  <*. 


T  T 

'  out  t/c.out 


Figure  5.1  Temperature  Rise  Corrections. 

To  take  account  of  the  heat  transferred  to  the  coolant  from  the  unheated  ends, 
the  heat  transfer  rate  to  the  coolant  over  the  entire  length  of  the  tube  (heated  and 
unheated)  is  first  calculated.  The  average  heat  transfer  coefficient  can  then  be 
determined  by: 


K  = 


IffL 


(5.5) 


AtLMTD 

where  A,  is  the  inside  surface  area  based  on  the  heated  length  (LH)  and  LMTD  is 
the  log  mean  temperature  difference  defined  as: 

At  this  point  in  the  calculation,  however,  the  heat  transfer  rate  to  the  coolant 
from  the  heated  length  alone  (cJhl)  is  not  known.  Since  the  heating  elements  did  not 
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T    -T. 
LMTD  =  _ 


In 


7.-7 

7,-7_, 


(5.6) 


cover  the  entire  length  of  the  tube,  account  must  be  taken  of  heat  added  to  the 
coolant  from  the  unheated  lengths  at  the  inlet  and  outlet  of  the  tube.  This  was 
accomplished  by  considering  these  ends  as  annular  fins.  An  adiabatic  fin  tip  and 
outer  wall  boundary  condition  was  applied  and  the  heat  transfer  from  the  fins,  qf  (as 
given  by  Incopera  and  Dewitt  [Ref.  10]),  determined  from: 


qf  -  y  HPkAc    9btanhmLf 


where  Lf  =  fin  length 
P  =  fin  perimeter  =  irDi 


Ac  =  fin  cross  sectional  area 


4 


0b  =  Twj-T^  (for  inlet  fin) 
=  T ^-T,,,  (for  outlet  fin) 


m  = 


kA 


10  J 
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The  fin  calculation  requires  a  known  average  inside  heat  transfer  coefficient  (h;). 
Therefore,  an  initial  estimate  of  h4  is  calculated  by: 

Ti.  =       ktubt  (5.8) 

'       AtLMTD 

where  q^  is  known.  The  amount  of  heat  that  passes  through  the  unheated  ends  is 

then  calculated  using  the  fin  correction  above.  The  heat  transfer  rate  to  the  coolant 

from  the  heated  length  (qmj  is  then  calculated  by: 

4«.  «  4-.-*/  <5-9) 

A  new  average  inside  heat  transfer  coefficient  (h{)  can  then  be  calculated  using 
equation  5.5  and  an  iteration  process  carried  out  until  convergence  to  a  final  average 
inside  heat  transfer  coefficient  is  reached. 

An  additional  correction  for  the  Cu/Ni  finned  tube  is  needed  to  account  for  the 
change  in  inside  diameter  of  the  smooth  ends  (needed  to  fit  through  the  tube  sheet 
in  the  condenser).  The  heat  transfer  coefficient  based  on  the  smaller  inside  diameter 
corresponding  to  the  finned  portion  of  the  tube  is  what  ultimately  is  desired.  As 
shown  in  Figure  5.2,  the  heated  length  extended  beyond  the  smaller  diameter  section. 
To  correct  for  this,  the  annular  fin  correction  discussed  earlier  was  applied  to 
determine  the  temperatures  at  the  inlet  and  outlet  of  the  heated  length.  The  coolant 
temperature  rise  along  the  heated  length  was  assumed  to  be  linear,  and  therefore  the 
temperature  rise  along  the  smooth  larger  inside  diameter  lengths  (L{  and  Lj) 'could 
be  determined  by  the  simple  ratio: 
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A7\  =  _i ^    am*    AT2  =      2     "L 


(5.10) 


where  Lhl  is  the  heated  length  of  the  tube  and  AThl,  the  temperature  rise  along  the 
heated  length.  These  corrections  were  then  added  or  subtracted  from  the  corrected 
finned  inlet  and  outlet  temperatures. 


Cu/Ni  EXTERNALLY  FINNED  TUBE 
TEMPERATURE  RISE  CORRECTIONS 
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Figure  5.2   Cu/Ni  Finned  Tube  Temperature  Rise  Corrections. 

With  the  average  inside  heat  transfer  coefficient  known,  the  mean  Nusselt 
number  can  be  determined  by: 


iV«a  =  -LJ. 
k 


(5.11) 


where  k  is  the  thermal  conductivity  of  the  ethylene  glycol/water  mixture. 
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The  experimentally  determined  mean  Nusselt  number  was  then  compared  to 
various  theoretical  and  experimental  correlations,  based  on  the  insert  condition* 
These  correlations  were  discussed  earlier  but  are  briefly  mentioned  here. 

For  the  no  insert  condition,  the  data  were  compared  to  two  theoretical 
correlations,  the  Hausen  correlation  [Ref.  10]  given  by: 

Nu       ^      O/KfiKD/W* 

U0M[(DIL)ReDPr]m 

and  the  Shah  and  London  correlation  [Ref.  9]  given  by: 

Num  =  1.93  (x-y1*  (5.13) 

The  reader  is  reminded  that  the  Hausen  correlation  is  for  a  uniform  wall 
temperature  boundary  condition  and  was  used  in  Mazzone's  condensation 
experiments  [Ref.  6].  The  Shah  and  London  correlation  is  for  a  constant  wall  heat 
flux  condition.  It  should  also  be  noted  that  equation  5.13  is  used  since  x"  for  the  test 
section  was  always  less  than  0.03. 

For  the  twisted  tape  insert  the  experimental  mean  Nusselt  number  was 
compared  to  the  correlation  developed  by  Hong  and  Bergles  [Ref.  15]: 

Num  =  Q3$3Pr035(ReJy)0-622  (5- 14) 

Equation  5.14  is  used  because  values  for  Re,/y  were  always  greater  than  10. 
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For  the  HEATEX  insert,  the  data  were  compared  to  the  two  correlations 
developed  by  Mazzone  [Ref.  6].  For  the  13.26  mm  inside  diameter  tube,  the  equation 
is  given  by: 

Num  =  0.226Re065Pr0M  (5.16) 

and  for  the  10.16  mm  inside  diameter  tube  by: 

Nu     =  0.063Re0J6Pr0A6  (5.16) 


B.      ADDITIONAL  TUBE  LENGTHS 

In  addition  to  the  above  calculations  for  the  entire  tube  length,  mean  Nusselt 


numbers  for  three  additional  tube  "lengths"  could  be  determined.  This  was 
accomplished  by  sectioning  the  heated  length  of  the  tube  into  four  cells  as  shown  in 
Figure  5.3.  Each  cell  contained  one  longitudinal  thermocouple  position.  Three 
additional  tube  lengths  were  then  determined  by  adding  cell  lengths  together.  The 
first  length,  Lv  consisted  of  cell  1  alone,  L  consisted  of  cells  1  and  2  and  L3 
consisted  of  cells  1,  2  and  3.  The  fourth  length  was  actually  the  entire  heated  length 
as  previously  discussed. 

The  coolant  temperature  rise  along  the  heated  length  of  the  tube  was  assumed 
to  be  linear.  The  correction  to  the  wall  heat  flux  to  account  for  variation  in  electrical 
resistance  with  temperatures  was  assumed  to  be  negligible.  Therefore,  the 
temperature  rise  along  each  cell  could  be  determined  by: 
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SECTIONED  TUBE 


Lm.- 


T/C  POS  1  T/C  POS  2  T/C  POS  3  T/C  POS  4 


INLET 


CELL  1 


CELL    2 


CELL  3 


CELL  * 


OUTLET 


Figure  5.3  Additional  Tube  Lengths. 


AT     -   47*« 


Ctf 


4 


(5.17) 


Average  inside  wall  temperatures  were  based  on  the  thermocouples  within  that  tube 
"length."  Therefore,  the  temperature  rise,  heat  flux  and  average  inside  wall 
temperatures  were  all  known  for  each  tube  length.  Applying  the  same  calculation 
method  as  before,  the  mean  Nusselt  number  for  these  three  additional  tube  "lengths" 
could  be  calculated. 

This  adaption  to  the  program  was  not  done  until  after  the  smooth  tube 
experiments  had  been  completed;  however,  this  extra  data  were  taken  for  the 
remaining  two  tubes. 
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C.      LOCAL  NUSSELT  NUMBERS 

The  program  was  also  adapted  to  determine  local  inside  heat  transfer 
coefficients  at  the  four  axial  thermocouple  positions.  Corrections  discussed  earlier 
to  determine  the  temperatures  at  the  inlet  and  outlet  of  the  heated  length  were  first 
accounted  for.  Exact  thermocouple  longitudinal  positions  were  entered  into  the 
program  as  L/D  ratios  so  that  a  local  x*  could  also  be  calculated.  Again,  the 
temperature  rise  along  the  heated  length  was  assumed  to  be  linear;  therefore,  the 
bulk  coolant  temperature  at  the  particular  thermocouple  position  could  be 
determined  by  the  simple  ratio  method  discussed  earlier  (equation  5.10). 

The  local  heat  transfer  coefficient,  h^,  was  then  determined  by: 

a" 
K  =  -=^r  (513) 

where  T^  is  the  average  inside  wall  temperature  based  on  the  three  thermocouples 
at  that  position  and  Tc  is  the  bulk  coolant  temperature  at  that  position.  The  local 
Nusselt  number  was  then  given  by: 

Nu    =  -2-L  (5- 19) 

k 

The  local  Nusselt  number  for  all  four  thermocouple  positions  was  then  determined. 

For  the  smooth  tube  experiments,  local  Nusselt  numbers  were  calculated  only 

for  position  1.  The  program  was  later  adapted  to  include  local  calculations  at 

positions  2,  3,  and  4  for  the  Cu/Ni  finned  tube  and  Korodense  tube. 


VI.   RESULTS  AND  DISCUSSION 

A.      GENERAL 

The  experimental  heat  transfer  results  are  initially  presented  in  terms  of  the 
mean  Nusselt  and  Reynolds  numbers  for  all  insert  conditions.  From  these  graphs,  it 
can  be  seen  how  the  Nusselt  number  is  inlet  temperature  dependent,  indicating  the 
need  for  a  Prandtl  number  (raised  to  an  exponent)  in  the  final  correlation.  However, 
although  not  specifically  shown,  the  data  shows  no  dependence  on  the  heat  flux 
setting.  In  addition,  the  Nu  vs.  Re  graphs  for  the  no  insert  condition  show  when  a 
transition  to  turbulent  flow  occurs. 

The  heat  transfer  results  are  then  compared  to  existing  correlations  based  upon 
the  insert  condition.  For  the  no  insert  condition,  the  Nusselt  number  is  plotted 
against  x*  (dimensionless  x  coordinate).  Note  that  for  a  given  tube,  L/D  is  constant 
and  x"  varies  due  to  Re  and  Pr  changes.  The  data  are  then  compared  to  the  Hausen 
[Ref.  10]  correlation  (equation  2.10)  and  the  Shah  and  London  [Ref.  9]  correlation 
(equation  2.14).  Recall  that  the  Hausen  correlation  is  for  a  uniform  wall  temperature 
boundary  condition  and  was  used  by  Mazzone  [Ref.  6]  in  the  refrigerant 
condensation  experiments.  Since  the  present  data  were  taken  under  a  constant  heat 
flux  boundary  condition,  comparison  with  the  Shah  and  London  correlation  was  also 
required.  The  twisted  tape  insert  data  are  compared  to  the  correlation  developed  by 
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Hong  and  Bergles  [Ref.  15]  (equation  2.20)  which  was  also  used  by  Mazzone.  The 
HEATEX  insert  data  are  compared  to  the  two  correlations  developed  from  the 
limited  data  received  from  CAL-GAVIN  by  Mazzone  (equations  2.27  and  2.28). 

Correlations  are  developed  for  each  tube  and  insert  condition  based  upon  a 
least  squares  linear  regression.  An  r-squared  value,  a  measure  of  the  accuracy  of  the 
least  squares  curve  fitting  to  the  data,  is  used  to  determine  the  uncertainty  of  the 
correlation.  The  correlations  and  the  corresponding  r-squared  values  are  summarized 
in  Table  6.1.  Note  that  the  0.2  m/s  and  0.3  m/s  data  were  eliminated  prior  to 
developing  all  correlations  due  to  the  high  degree  of  uncertainty  in  these  data  points. 
TABLE  6.1   SUMMARY  OF  CORRELATIONS 


SUMMARY  OF  CORRELATIONS 


NO  INSERT 

TWISTED 
TAPE 

HEATEX 

SMOOTH 

TUBE 

Nu  - 
4.137Xstar~-0.234 

r-sq  "  71.9% 

Nu  • 
.5l2Pr*.35(Rea/y)*.568 

r-sq  -  99.5% 

Nu  • 
.164Pr.4efle~.64i 

r-sq  -  98.9% 

EXT 

FINNED 
TUBE 

Nu  - 
1.84Xstar~-0.34i 

r-sq  -  93.0% 

Nu  " 
.148Pr~35(Res/yr.e87 

r-sq  -  96.7% 

Nu  - 
.l26Pr'46fle'657 

r-sq  ■  95.9 

KORODENSE 
TUBE 

Nu  - 
2.047Xstar -0.326 

r-sq  -  89.2% 

Nu  - 
.569Prr35(P.9S/yr.498 

r-sq  -  98.9% 

Nu  - 
.282Pr~.4eRe~.499 

r-sq  -  94.3% 
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B.      SMOOTH  TUBE 

1.      Smooth  Tube  -  No  Insert 

The  first  tube  tested  was  the  copper  smooth  tube  with  no  insert.  Heat 
transfer  data  is  presented  in  Figure  6.1  in  terms  of  the  Nusselt  number  as  a  function 
of  Reynolds  number  for  each  nominal  inlet  temperature.  The  first  thing  to  notice  is 
the  sudden  increase  in  Nu  at  a  Re  number  of  about  2300.  This  is  indicative  of  a 
change  from  laminar  to  transition  flow  where  the  heat  transfer  is  greatly  increased 
due  to  improved  mixing.  This  contradicts  Mazzone's  [Ref.  6]  assumption  of  a 
transition  to  turbulence  occurring  at  Re  >  4000.  The  data  in  the  transition  region 
was  primarily  for  the  +20°C  and  +10°C  nominal  inlet  temperatures.  Figure  6.1  also 
shows  distinct  groupings  of  data  at  given  inlet  temperatures,  indicating  a  Prandtl 
number  effect.  These  groupings  are  shown  more  clearly  in  Figure  6.2  which  represent 
data  under  laminar  flow  conditions  only  (Re  ^  2300)  and  Figure  6.3  for  data  in  the 
transition  region  only  (Re  >  2300). 

To  develop  a  useful  heat  transfer  correlation,  the  mean  Nusselt  number 
is  plotted  as  a  function  of  x"  in  Figure  6.4.  Recall  that  x"  can  be  determined  by: 

*•-_££_  (2.9) 

ReDPr 

For  a  given  tube,  L/D  is  a  fixed  quantity.  Therefore,  for  a  given  Prandtl  number 
(essentially  fixed  inlet  temperature),  x"  is  dependent  only  upon  the  Reynolds  number. 
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Figure  6.4  also  shows  a  comparison  of  the  experimental  data  with  the 
Hausen  correlation  (equation  2.10)  for  a  uniform  wall  temperature  (used  by 
Mazzone)  and  the  Shah  and  London  correlation  (equation  2.14)  for  a  constant  heat 
flux.  The  Hausen  correlation  appears  to  significantly  underestimate  the  heat  transfer 
coefficient  found  in  the  present  experimental  data,  i.e.,  the  present  experimental  data 
indicate  a  lower  inside  thermal  resistance.  This  is  consistent  with  Mazzone's  [Ref.  6] 
results  in  that  the  Hausen  correlation  was  predicting  larger  inside  resistances  than 
the  overall  measured  resistance.  However,  it  must  always  be  remembered  that  the 
Hausen  correlation  is  for  a  uniform  wall  temperature  boundary  condition  and  this 
would  indeed  lead  to  lower  heat-transfer  coefficients. 

The  data  also  fall  slightly  above  (but  parallel  to)  that  predicted  by  the 
Shah  and  London  correlation.  Due  to  the  constant  heat  flux  boundary  condition  used 
by  Shah  and  London,  one  may  expect  better  agreement  (as  seen).  However,  the  fact 
that  the  data  fall  slightly  above  can  be  attributed  to  secondary  flow  induced  by  the 
apparatus  inlet  conditions  (90°  bend  due  to  the  thermocouple  well).  Also,  as  before, 
a  departure  from  laminar  flow  can  be  seen  in  Figure  6.4.  Since  the  transition  data 
forms  two  distinct  lines  (+20°C  and  +10°C),  it  indicates  that  x*  is  not  a  good 
correlating-  parameter  for  transitional  flow. 

After  removal  of  the  transitional  flow  data,  a  correlation  was  developed 
based  upon  the  form  of  the  Shah  and  London  correlation.  The  resulting  present 
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correlation  for  the  mean  Nusselt  number  for  laminar  flow  with  no  insert  is  given  by: 

Num  =  4.137 (x')-0**  (6.1) 

Equation  6.1  and  the  Hausen  correlation,  equation  2.10,  are  both  shown 
in  Figure  6.5  with  the  present  laminar  flow  data.  The  general  form  of  Figure  6.5  is 
explained  by  the  fact  that  at  low  x"  (high  Re)  the  heat  transfer  is  improved  due  to 
slightly  better  mixing.  As  x*  increases  (coolant  velocity  decreases),  the  heat  transfer 
decreases. 

Figure  6.5  also  shows  a  high  degree  of  scatter  in  the  data.  Initially,  it  was 
thought  that  this  scatter  was  due  to  the  flow  not  being  hydrodynamically  fully 
developed.  However,  it  is  seen  to  be  more  pronounced  at  large  x*.  Recall  that  for  a 
given  nominal  inlet  temperature,  large  x"  represents  a  low  Reynolds  number.  For  a 
low  Reynolds  number,  the  hydrodynamic  entry  length  (given  by  X/D  =  0.05Re)  is 
the  shortest.  Also,  since  we  are  dealing  with  relatively  large  Prandtl  numbers  (30- 
140),  the  assumption  that  the  velocity  profile  is  fully  developed  throughout  the 
thermal  entry  region  should  be  fairly  valid.  Therefore,  it  is  felt  that  this  scatter  is  not 
due  to  the  flow  being  hydrodynamically  undeveloped. 

This  scatter  is  believed  to  be  due  to  the  effects  of  the  different  inlet 
temperatures,  indicating  that  x*  is  not  an  ideal  correlating  parameter.  That  is  to  say, 
the  Prandtl  number  exponent  needs  to  be  adjusted  to  account  for  the  different  inlet 
temperatures.  It  should  also  be  noted  that  due  to  the  high  degree  of  detail  in  Figure 
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6.5,  the  scatter  appears  more  pronounced.  However,  in  Figure  6.4  with  larger  scale 
increments,  the  scatter  seems  minimal. 

2.  Smooth  Tube  -  Twisted  Tape 

Results  for  the  smooth  tube  with  a  twisted  tape  insert  are  shown  in 
Figures  6.6  and  6.7.  A  definite  Prandtl  number  dependence  can  be  seen  in  Figure  6.6 
where  the  mean  Nusselt  number  is  presented  as  a  function  of  Reynolds  number. 
Following  the  guidance  of  Hong  and  Bergles  [Ref.  15]  the  dimensionless  parameters 
Re,  and  y  were  combined  as  one  parameter  Re^y  and  the  mean  Nusselt  number  and 
Prandtl  number  were  combined  as  Nu/Pr0J5.  The  final  correlation  shown  in  Figure 
6.7  is: 

Num  =  0.5\2Pr035(ReJy)osa  (6.2) 

This  correlation  is  in  close  agreement  with  the  Hong  and  Bergles  correlation 
(equation  2.20).  The  slight  differences  can  be  attributed  to  differences  in  the 
experimental  apparatus  (including  the  insert  itself)  and  the  possible  existence  of 
secondary  flows.  Figure  6.7  shows  little  inlet  temperature  dependence,  indicating  that 
the  Prandtl  number  exponent  of  0.35  is  probably  fairly  accurate.  Furthermore,  there 
is  little  experimental  data  scatter,  giving  increased  confidence  in  this  correlation. 

3.  Smooth  Tube  -  HEATEX  Insert 

HEATEX  insert  results  for  a  smooth  tube  are  shown  in  Figures  6.8  and 
6.9.  Again,  on  a  Nu  vs.  Re  basis,  it  can  be  seen  from  Figure  6.8  that  there  is  a 
definite  dependence   on   the   inlet  temperature.    Using   the   same   form   as   the 
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correlation  of  Oliver  and  Aldington  [Ref.  20],  the  resulting  correlation  for  the  mean 
Nusselt  number  is: 

Nu     =  0A64Pr0MRe0Ml  (63) 

This  is  compared  in  Figure  6.9  to  the  correlation  developed  and  used  by  Mazzone 
(equation  2.23).  It  is  clearly  evident  that  Mazzone's  correlation  underpredicts  the 
inside  thermal  resistance  and  therefore  gives  a  larger  inside  heat  transfer  coefficient. 
Again,  Figure  6.9  shows  little  dependence  on  inlet  temperature  indicating  that  the 
Prandtl  number  exponent  is  fairly  accurate.  Comparisons  of  the  enhancement  in  heat 
transfer  gained  by  the  use  of  HEATEX  and  twisted  tape  insert,  over  the  no  insert 
condition  are  given  section  F  of  this  chapter. 
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Figure  6.1   Smooth  tube  -  No  insert  (Num  vs.  Re)  results. 
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Figure  6.2   Smooth  tube  -  No  insert  (Num  vs.  Re)  laminar  flow  results. 
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Figure  6.3   Smooth  tube  -  No  insert  (Num  vs.  Re)  transitional  flow  results. 
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Figure  6.4   Smooth  tube  -  No  insert  (Num  vs.  x")  results. 
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Figure  6.5   Smooth  tube  -  No  insert  laminar  flow  correlation  results. 
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Figure  6.6   Smooth  tube  -  Twisted  tape  (Num  vs.  Re)  results. 
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Figure  6.7   Smooth  tube  -  Twisted  tape  correlation  results. 


69 


200 


150 


i 


100 


50 


SMOOTH  TUBE  -  HEATEX 

Nu  vs  Re 


♦ 


A    A 


A    A  ♦ 

♦ 


*  A 


A  ♦  ■ 

A 


A*      «♦     a 


* 


1000 


2000  3000 

Reynolds 


4000 


5000 


+20C     ♦  +10C     A   0C         ,    -IOC 


SMHXNURE 


Figure  6.8   Smooth  tube  -  HEATEX  (Nu^,  vs.  Re)  results. 
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Figure  6.9   Smooth  tube  -  HEATEX  correlation  results. 
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C.      CU/NI  EXTERNALLY  FINNED  TUBE 

1.      Cu/Ni  Finned  Tube  -  No  Insert 

The  second  tube  tested  was  the  Cu/Ni  externally  finned  tube.  Recall  that 
this  tube  has  a  smaller  inside  smooth  diameter  than  the  smooth  tube  (10.15  mm  vs. 
13.26  mm).  The  no  insert  results  for  the  mean  Nusselt  number  as  a  function  of 
Reynolds  number  are  shown  in  Figure  6.10  and  in  more  detail  in  Figures  6.11  and 
6.12  for  laminar  and  transitional  flows,  respectively.  It  can  be  seen  that  a  slightly 
higher  Reynolds  number  (Re  «  2700)  for  the  departure  from  laminar  flow  exists  for 
this  tube.  Again,  the  dependence  on  inlet  temperature  is  clearly  evident 

The  data  are  compared  to  the  Hausen  correlation  (equation  2.10)  and  the 
Shah  and  London  correlation  (equation  2.14)  in  Figure  6.13.  As  expected,  the 
Hausen  correlation  again  underpredicts  the  inside  heat  transfer  coefficient.  However, 
the  laminar  flow  data  is  in  closer  agreement  with  the  Shah  and  London  correlation 
than  they  were  for  the  smooth  tube.  This  can  be  attributed  to  the  different  diameter 
affecting  the  thermal  boundary  layer  development.  The  resulting  correlation  for  the 
laminar  flow  mean  Nusselt  number  is  shown  in  Figure  6.14  and  is  given  by: 

Num  =  1.84(.r-)-0J41  (6.4) 

The  same  conclusions  about  data  scatter  can  be  drawn  as  for  the  smooth  tube 
results,  although  in  general  the  scatter  seems  to  be  less. 
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2.  Cu/Ni  Finned  Tube  -  Twisted  Tape  Insert 

Figures  6.15  and  6.16  show  the  results  for  the  Cu/Ni  finned  tube  using  a 
twisted  tape  to  enhance  tubeside  heat  transfer.  The  inlet  temperature  dependence 
in  Figure  6.15  seems  to  be  reduced,  but  still  evident  The  form  of  the  correlation  was 
again  that  of  Hong  and  Bergles'  correlation.  It  was  developed  in  the  same  manner 
as  for  the  smooth  tube  and  is  given  by: 

Num  =  Q.MPr^tReJy)0-6*7  (6.5) 

The  data  and  resulting  correlation  are  compared  to  the  Hong  and  Bergles  correlation 
(equation  2.20)  in  Figure  6.16  where  the  latter  is  seen  to  overpredict  the  mean 
Nusselt  number.  This  difference  can  be  attributed  to  several  factors.  A  small 
diameter  tube  with  a  larger  diameter  smooth  entrance  length  could  cause  severe 
secondary  flows.  Also,  because  of  the  smaller  diameter  tube  a  different  size  twisted 
tape  (twist  ratio,  y=3)  was  used,  which  had  a  slightly  looser  fit  than  the  tape  used 
with  the  smooth  tube.  Additionally,  Prandtl  number  effects  need  to  be  investigated 
more  closely,  as  Figure  6.16  shows  a  definite  trend  with  inlet  temperature. 

3.  Cu/Ni  Finned  Tube  -  HEATEX  Insert 

Data  for  the  Cu/Ni  finned  tube  using  the  HEATEX  insert  are  presented 
in  Figures  6.17  and  6.18.  Figure  6.17  (Nu  vs.  Re)  again  shows  a  dependence  of  the 
data  on  the  coolant  inlet  temperature.  Since  the  HEATEX  element  used  for  this 
tube  was  a  smaller  diameter  (10.16  mm)  with  a  different  loop  density  than  that  used 
for  the  smooth  tube,  a  new  correlation  needed  to  be  developed.  Again,  using  the 
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form  of  the  correlation  proposed  by  Oliver  and  Aldington,  the  final  proposed  correlation 
is  given  by: 

Num  -  0A26Pr0MReQ^  (6.6) 

The  data  and  resulting  correlation  are  compared  to  the  correlation  developed 
by  Mazzone  (equation  2.24)  for  this  insert  in  Figure  6. 18.  It  can  be  seen  that  Mazzone's 
correlation  is  adequate  for  predicting  the  inside  heat  transfer  coefficients  at  low  Reynolds 
number  (Re  <  1400),  where  both  correlations  merge.  However,  at  higher  Reynolds 
numbers,  Mazzone's  correlation  overpredicts  the  present  experimental  data.  The 
differences  in  the  two  correlations  are  primarily  due  to  the  limited  amount  of  data 
available  to  Mazzone  when  developing  his  correlation.  Furthermore,  the  data  he  used 
were  taken  on  a  different  experimental  facility  (by  the  manufacturers).  As  seen  before, 
the  data  shown  in  Figure  6. 18  still  indicate  a  dependence  on  inlet  temperature  and  hence 
the  value  of  0.46  for  the  Prandtl  number  exponent  needs  to  be  reinvestigated  and 
optimized. 
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Figure  6.10   Cu/Ni  externally  finned  tube  -  No  insert  (Num  vs.  Re)  results. 
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Figure  6.11    Cu/Ni  externally  finned  tube  -  No  insert  (NuB  vs.  Re)  laminar  flow 
results. 
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Figure  6.12  Cu/Ni  externally  finned  tube  -  No  insert  (Num  vs.  Re)  transitional  flow 
results. 
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Figure  6.13   Cu/Ni  externally  finned  tube  -  No  insert  (Num  vs.  x")  results. 
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Figure  6.14  Cu/Ni  externally  finned  tube  -  No  insert  laminar  flow  correlation  results. 


79 


Cu/Ni  FINNED  TUBE  -  TWISTED  TAPE 

Nu  vs  Re 


110 

- 

a 

80 

a 

100 

a 
a 
■■ 

a  " 

90 

♦   ■ 
.1 

as  so 

- 

a  bV 

u 

3 

A   ■  ♦ 

- 

Z  70 

- 

AA  ^S* 

§ 

A  4  ♦ 

S  60 

- 

•  A  % 

A   ▼ 

_ 

'    A* 

50 

- 

*   * 

■» 

40 

„  '  jA  B 

*■  s 

,  a  a 

A  J 
^A     + 

30 

— 

'   A  ■» 

.  A 

20 

1                1                I.I                1 

1000  2000  3000 

Reynolds 


4000 


5000 


6000 


+15       ♦  +10C     A  0C         ,    -IOC 


CFTTNURE 


Figure  6.15   Cu/Ni  externally  finned  tube  -  Twisted  tape  (Num  vs.- Re)  results. 
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Figure  6.16   Cu/Ni  externally  finned  tube  -  Twisted  tape  correlation  results. 
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Cu/Ni  FINNED  TUBE  -  HEATEX 
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Figure  6.17   Cu/Ni  externally  finned  tube  -  HEATEX  (Num  vs.  Re)  results. 
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Cu/Ni  FINNED  TUBE  -  HEATEX 
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Figure  6.18   Cu/Ni  externally  finned  tube  -  HEATEX  correlation  results. 
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D.      KORODENSE  TUBE 

1.      Korodense  Tube  -  No  Insert 

The  final  tube  tested  was  the  Korodense  tube.  This  tube  was  first  tested 
with  no  insert  The  results  for  the  mean  Nusselt  number  as  a  function  of  Reynolds 
number  are  shown  in  Figure  6. 19  and  in  more  detail  for  laminar  and  transitional  flow 
in  Figures  6.20  and  6.21,  respectively.  As  was  seen  in  the  smooth  tube  case,  a 
transition  to  turbulent  flow  occurred  at  a  Reynolds  number  of  approximately  2300. 

Figure  6.22  compares  the  data  to  Hausen's  correlation  (equation  2.10)  and 
the  Shah  and  London  correlation  (equation  2.14).  As  before,  the  transitional  flow 
data  are  seen  as  two  distinct  lines  for  data  taken  at  +20°C  and  +  10°C  nominal  inlet 
temperatures.  Under  laminar  flow  conditions,  the  data  are  seen  to  fall  very  closely 
to  Shah  and  Londons'  correlation,  but  significantly  above  Hausen's  correlation.  This 
was  also  found  for  the  other  two  tubes.  This  result  was  expected  since  a  mildly 
corrugated  wall  is  not  good  in  enhancing  laminar  flow  and  the  tube  behaves 
essentially  like  a  smooth  tube.  Significant  enhancement  for  Korodense  tubes  is  not 
expected  until  ReD  >  104  [Ref.  21]. 

The  resulting  correlation  for  laminar  flow  data  is  shown  on  Figure  6.23 
and  given  by: 

Num  =  2.05(O-0J26  (6.7) 

The  Hausen  correlation  (equation  2.10)  used  by  Mazzone  again  under  predicts  the 
present  measured  inside  heat  transfer  coefficient.  It  was  expected  that  the  mean 
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Nusselt  number  would  be  the  same  (or  maybe  slightly  higher)  as  for  the  smooth  tube 
case  for  the  reasons  already  mentioned.  However,  the  Korodense  data  fall  slightly 
below  the  smooth  tube  data  and  closer  to  Shah  and  Londons'  smooth  tube 
correlation  (equation  2.14).  No  explanation  can  be  given  for  this  result.  Data  scatter 
explanations  given  for  the  other  tubes  also  apply  for  the  Korodense  tube. 

2.      Korodense  Tube  -  Twisted  Tape  Insert 

The  results  for  the  Korodense  tube  with  a  twisted  tape  insert  are  shown 
in  Figures  6.24  and  6.25.  The  same  twisted  tape  that  was  used  for  the  smooth  tube 
was  used  for  the  Korodense  tube.  Note  that  the  Korodense  tube  nominal  inside 
diameter  (13.35  mm)  is  slightly  larger  than  the  smooth  tube  inside  diameter  (13.26 
mm).  This  small  difference  in  inside  diameters  (together  with  corrugation  in  the 
Korodense  tube)  led  to  a  looser  fit  for  the  twisted  tape  insert. 

The  mean  Nusselt  number  as  a  function  of  Reynolds  number  is  shown  in 
Figure  6.24  which  again  clearly  indicates  the  inlet  temperature  dependence.  The  data 
were  again  correlated  using  the  form  of  Hong  and  Bergles'  correlation.  The  present 
data  were  compared  to  their  correlation  (equation  2.20)  in  Figure  6.25.  The 
correlation  developed  for  the  present  data  is: 

Num  =  0.569Pr035(ReJy)0A9(i  (6:8) 

This  correlations  falls  below  the  Hong  and  Bergles  correlation  (equation  6.20)  as  well 
as  for  the  smooth  tube.  This  result  is  thought  to  be  primarily  due  to  the  looser  fit  of 
the  twisted  tape  inside  the  Korodense  tube.  The  corrugations  also  leave  areas  of  no 
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contact  between  the  tape  and  the  inside  wall  of  the  tube.  Stagnant  flow  or  even  air 
bubbles  may  form  in  these  gaps,  especially  at  low  Reynolds  numbers,  causing  an 
interruption  of  the  swirled  flow.  The  thermal  boundary  layer,  therefore,  is  not  as 
efficiently  mixed  as  in  the  smooth  tube.  For  these  reasons,  lower  inside  heat  transfer 
coefficients  than  those  assumed  by  Mazzone  should  be  expected. 

3.      Korodense  Tube  -  HEATEX  Insert 

The  last  tube-insert  condition  tested  was  the  Korodense  tube  with  the 
HEATEX  insert  The  same  HEATEX  insert  used  for  the  smooth  tube  was  used 
here.  Figure  6.26  shows  the  mean  Nusselt  number  as  a  function  of  Reynolds  number 
and  again  shows  an  inlet  temperature  dependence.  The  resulting  correlation  for  the 
present  data  is  given  by: 

Nu     =  O^S^r0-46^0499  (6.9) 

The  data  and  this  correlation  are  compared  with  the  correlation  Mazzone  used 
(equation  2.23)  in  Figure  6.27.  It  is  evident  that  Mazzone's  chosen  correlation 
overpredicts  the  inside  heat  transfer  coefficient  significantly.  This  again  is  thought 
to  be  primarily  due  to  the  incomplete  mixing  of  the  thermal  boundary  layer  caused 
by  non-contact  regions  existing  between  the  tube  wall  and  HEATEX  petals. 
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Figure  6.19   Korodense  tube  -  No  insert  (Num  vs.  Re)  results. 
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Figure  6.20   Korodense  tube  -  No  insert  (Nu,,,  vs.  Re)  laminar  flow  results. 
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Figure  6.21   Korodense  tube  -  No  insert  (Num  vs.  Re)  transitional  flow  results. 
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Figure  6.22   Korodense  tube  -  No  insert  (Num  vs.  x")  results. 
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Figure  6.23   Korodense  tube  -  No  insert  laminar  flow  correlation  results. 
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KORODENSE  TUBE  -  TWISTED  TAPE 
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Figure  6.24   Korodense  tube  -  Twisted  tape  (Num  vs.  Re)  results. 
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Figure  6.25    Korodense  tube  -  Twisted  tape  correlation  results. 
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Figure  6.26   Korodense  tube  -  HEATEX  (Num  vs.  Re)  results. 
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Figure  6.27   Korodense  tube  -  HEATEX  correlation  results. 


95 


E.      INSERT  COMPARISON  BETWEEN  TUBES 

Figures  6.28,  6.29  and  6.30  each  compare  the  three  correlations  for  each  tube 
for  no  insert,  twisted  tape,  and  HEATEX  inserts,  respectively.  It  can  be  seen  that 
with  no  insert  (Figure  6.28),  the  Korodense  correlation  is  about  10%  below  the 
smooth  tube  data.  Although  one  may  expect  no  increase  in  the  heat  transfer  with  a 
roped  tube  under  laminar  flow  conditions,  one  certainly  would  not  expect  to  find 
such  a  decrease  with  a  similar  diameter  tube.  This  difference,  therefore,  can  only  be 
attributed  to  experimental  uncertainty,  probably  in  the  wall  temperature 
measurement  The  finned  tube  correlation  is  lower  than  both  the  Korodense  («5%) 
and  smooth  tube  («15%)  correlations;  this  tube  had  a  smaller  inside  diameter  and 
it  may  be  that  different  induced  flow  patterns  are  set  up  causing  a  reduction  in  heat 
transfer.  Experimental  uncertainty  must  also  not  be  ruled  out  as  to  the  cause  of  the 
difference. 

Figure  6.29  shows  the  same  three  tube  correlations,  but  for  the  twisted  tape 
insert  The  smooth  tube  again  has  the  highest  performance.  The  finned  tube  is  =20% 
below  at  a  high  Re  number  and  nearly  40%  lower  at  a  low  Re  number.  This  could 
be  due  to  a  different  twist  ratio  used  for  the  finned  tube  due  to  the  smaller  inside 
diameter.  The  Korodense  tube  is  similar  to  the  finned  tube,  giving  slightly  lower  heat 
transfer  at  high  Re  number  and  slightly  higher  heat  transfer  at  low  Re.  The  reason 
for  this  tube  being  significantly  lower  than  the  smooth  tube  is  though  to  be  due  to 
the  tape  not  contacting  the  inner  surface  of  the  tube  as  well  as  for  the  smooth  tube, 
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thereby  leaving  regions  of  stagnant  flow.  These  regions  would  tend  to  inhibit  mixing 
and  lead  to  a  lower  heat  transfer  performance. 

For  the  HEATEX  comparison  (Figure  6.30),  similar  conclusions  can  be  made, 
although  the  finned  tube  seems  to  have  better  relative  performance  indicating  that 
the  small  diameter  HEATEX  insert  used  in  the  small  diameter  tube  is  more  effective 
than  the  twisted  tape.  The  difference  between  the  smooth  and  Korodense  tubes  are 
of  similar  magnitude  to  that  shown  in  Figure  6.29  (about  40%  at  high  Re). 
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Figure  6.28   No  insert  comparison  under  laminar  flow  for  all  three  tubes. 
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TWISTED  TAPE  INSERT  COMPARISON 
UNDER  LAMINAR  FLOW 
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Figure  6.29  Twisted  tape  insert  comparison  under  laminar  flow  for  all  three  tubes. 
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Figure  6.30   HEATEX  insert  comparisons  under  laminar  flow  for  all  three  tubes. 
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F.       LAMINAR  FLOW  ENHANCEMENT  RATIO 

The  comparisons  made  in  Figures  6.28  through  6.30  are  made  more 
complicated  by  the  fact  that  the  uncertainty  associated  with  the  wall  temperature 
measurement  is  unknown  and  most  likely  different  for  each  tube.  For  this  reason,  it 
is  not  known  how  much  of  the  divergence  seen  in  these  figures  is  due  to  differences 
in  thermocouple  attachment  (giving  uncertainty  to  the  wall  temperature 
measurement)  and  how  much  is  an  actual  flow  induced  phenomenon.  A  comparison 
of  the  relative  effects  of  the  three  inserts  for  a  given  tube  is  more  accurate,  however, 
due  to  the  fact  that  any  uncertainty  for  a  given  tube  wall  thermocouple  attachment 
combination  would  be  the  same. 

To  determine  the  degree  of  tubeside  heat  transfer  enhancement  achieved  by 
the  use  of  the  twisted  tape  and  HEATEX  inserts  under  laminar  flow,  an 
enhancement  ratio  of  the  mean  Nusselt  number  for  the  insert  condition  (Nu^,)  to 
that  for  the  no  insert  condition  (Nunouuert)  was  calculated.  The  enhancement  ratio  is 
therefore  defined  as: 


=     NuJ~~>  (6.10) 

no  insert 


Figures  6.31  through  6.33  present  mean  Nusselt  number  as  a  function  of  Reynolds 
number  under  laminar  flow  conditions  for  each  tube  under  the  three  insert 
conditions.  For  the  purposes  of  this  study,  enhancement  ratios  were  calculated  at  the 
two  extreme  Reynolds  numbers  of  200  and  2300  with  a  constant  Prandtl  number  of 
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70  using  the  correlations  developed  in  the  previous  sections.  Table  6.2  summarizes  the 
enhancement  ratios  for  each  tube  and  insert  condition. 

It  is  evident  from  these  figures  that  the  greatest  enhancement  is  achieved  using  the 
HEATEX  insert.  An  enhancement  factor  of  over  seven  was  achieved  for  the  smooth  tube 
using  the  HEATEX  insert  at  Re  =  2300.  This  reduced  to  about  2.6  at  low  Re. 
Significantly  smaller  enhancements  were  achieved  for  the  twisted  tape  inset  for  the 
smooth  tube,  indicating  that  for  laminar  flow,  the  stripping  and  mixing  effect  of  the 
HEATEX  element  is  considerably  more  effective  than  the  simple  swirl  flow  effect 
provided  by  the  twisted  tape.  It  should  be  mentioned  here  that  HEATEX  was  developed 
specifically  with  laminar  flow  enhancement  in  mind.  In  turbulent  flow  (Re  >  10*),  the 
enhancements  obtained  are  much  smaller  (around  1.5  to  2)  and  both  twisted  tape  and 
HEATEX  elements  give  similar  enhancement  values. 


TABLE  6.2     SUMMARY  OF  LAMINAR  FLOW  ENHANCEMENT  RATIOS 
(defined  by  equation  6.10) 


• 

SMOOTH  TUBE 

FINNED  TUBE 

KORODENSE  TUBE 

TWISTED 
TAPE 

HEATEX 

TWISTED 
TAPE 

HEATEX 

TWISTED 
TAPE 

HEATEX 

Low  Re 
(Re  =  200) 

1.79 

2.66 

1.49 

2.86 

1.92 

2.58 

High  Re 
(Re  =  2300) 

4.04 

7.19 

3.46 

6.71 

2.93 

4.25 
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For  the  finned  and  Korodense  tubes,  the  HEATEX  element  also  provides  the 
highest  enhancements.  At  high  Re,  the  enhancements  provided  by  the  twisted  tape 
and  HEATEX  inserts  are  lower  than  those  obtained  with  the  smooth  tube  for  the 
reasons  mentioned  previously.  However,  at  low  Re,  the  enhancements  for  the  three 
tubes  are  about  the  same  for  a  given  insert,  indicating  the  "stagnant"  regions 
mentioned  earlier  have  less  of  an  effect  on  the  total  heat  transfer. 
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INSERT  CONDITION  COMPARISON 
LAMINAR  FLOW  -  SMOOTH  TUBE 
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Figure  6.31    Insert  condition  comparison  for  laminar  flow  in  a  smooth  tube. 
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INSERT  CONDITION  COMPARISON 
LAMINAR  FLOW  -  Cu/Ni  FINNED  TUBE 
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Figure  6.32    Insert  condition  comparison  for  laminar  flow  in  a  Cu/Ni  externally 
finned  tube. 
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INSERT  CONDITION  COMPARISON 
LAMINAR  FLOW  -  KORODENSE  TUBE 
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Figure  6.33   Insert  condition  comparison  for  laminar  flow  in  a  Korodense  tube. 
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G.      REPROCESSED  REFRIGERANT  CONDENSATION  DATA 

Figure  6.34  shows  the  condensation  data  of  Mazzone  [Ref.  6]  with  the  outside  heat 
transfer  coefficient  recalculated  using  the  present  correlations.  The  data  presented  are  for 
the  top  copper  smooth  tube.  Similar  recalculations  can  also  be  carried  out  for  the  other 
two  tubes  (Cu/Ni  finned  and  Korodense  tubes)  as  well  as  for  condensation  bundle  data. 
The  data  of  Mazzone  are  reprocessed  using  the  Modified  Wilson  Plot  Technique  (see 
Mazzone  [Ref.  6]  for  details).  This  technique  essentially  lets  the  leading  coefficient  of 
the  correlation  'float'  to  try  to  smooth  out  any  uncertainty  in  the  correlations  being  used. 
As  previously  stated,  the  chosen  correlation  for  the  no  insert  condition  (even  when 
reprocessing  using  the  Wilson  Plot  Technique)  gave  values  of  inside  thermal  resistances 
greater  than  the  measured  overall  thermal  resistance,  thereby  yielding  negative  values  of 
the  outside  heat  transfer  coefficient,  h°.  This  explains  why  Figure  6.34  has  no  data  for 
Mazzone' s  no  insert  condition. 

It  can  be  seen  that  by  using  the  present  developed  correlations  for  the  no  insert 
condition,  the  calculated  values  of  h0  are  significantly  improved  (i.e.  positive  values). 
In  addition,  for  the  other  two  inserts,  the  new  correlations  increase  h„  by  about  10%  for 
the  twisted  tape  insert  and  nearly  20%  for  the  HEATEX  insert.  The  newly  calculated 
values  of  h0  when  using  these  inserts  are  now  very  close  and  just  above  the  Nusselt 
prediction  as  may  be  expected  for  low  vapor  velocity  film  condensation.  Mazzone 
commented      that      his      data      were      below      Nusselt      by      about      10% 
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and  thought  that  this  was  a  consequence  of  the  length  of  his  condenser  tubes  violating  the 
uniform  wall  temperature  assumption  of  the  Nusselt  theory.  The  fact  that  these  data  are 
so  close  gives  confidence  in  the  new  correlations,  as  the  insert  itself  should  have  no 
influence  on  the  magnitude  of  h,.  By  significantly  reducing  the  inside  resistance,  the 
insert  simply  increase  the  accuracy  of  h,.  Of  course,  the  no  insert  condition  should  also 
give  similar  values  of  h0.  Figure  6.34  shows  that  they  are  too  high,  indicative  of  the  fact 
that  the  correlation  for  the  no  insert  condition,  although  much  improved,  is  still  not  ideal. 
If  one  recalls  the  scatter  obtained  in  the  present  data  for  the  no  insert  case,  together  with 
the  fact  that  the  inside  resistance  is  still  dominating,  then  the  results  shown  in  Figure 
6.34  are  not  too  surprising  and  still  very  encouraging.  It  is  recommended,  therefore,  that 
all  future  refrigerant  condensation  experiments  carried  out  on  the  bundle  facility  should 
be  conducted  with  the  inserts  using  the  newly  formulated  correlations. 

One  final  note  is  the  fact  that  the  present  correlations  were  formulated 
using  a  constant  wall  heat  flux,  whereas  for  the  condensation  experiments,  one 
should  ideally  use  correlations  developed  under  uniform  wall  temperature  conditions. 
One  may  expect  that  the  inside  heat  transfer  coefficient  would  increase  slightly  under 
constant  wall  heat  flux  conditions,  thereby  decreasing  the  outside  heat  transfer 
coefficient.  However,  from  Figure  6.34  it  is  clear  that  the  use  of  the  constant  wall 
heat  flux  correlations  seems  to  increase  h0.  Thus  it  appears  that  with  the  inserts,  the 
inside  resistance   is   so  reduced   that  the  boundary  condition  imposed   has   little 
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influence.    Future  investigators  should  look  into  this  boundary  value  problem  in 
more  detail. 
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Figure  6.34   Reprocessed  refrigerant  condensation  data  using  new  correlations. 
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VII.  CONCLUSIONS 


The  following  conclusions  are  based  upon  results  obtained  in  this  study. 


GENERAL 


1.  A  test  apparatus  was  successfully  designed  and  built  to  accurately  determine 
inside  heat  transfer  coefficients  for  the  same  tubes  and  insert  conditions  used 
in  the  condensation  experiments. 

2.  Correlations  to  accurately  determine  inside  heat  transfer  coefficients  under 
laminar  flow  conditions  both  with  and  without  inserts  were  successfully 
developed. 

3.  Reprocessing  of  Mazzone's  [Ref.  6]  refrigerant  condensation  data  using  these 
new  correlations  has  yielded  more  accurate  values  of  the  outside  heat  transfer 
coefficient  (Note:  the  name  of  the  program  used  for  the  reprocessing  is 
DRPCON7.) 


B.      SPECIFIC 


1.  A  departure  from  laminar  flow  into  a  transitional  flow  exists  at  Reynolds 
numbers  of  approximately  2300  for  the  smooth  and  Korodense  tubes  and  2700 
for  the  Cu/Ni  finned  tube.  This  contradicts  Mazzone's  [Ref.  6]  assumption  of 
laminar  flow'  up  to  a  Reynolds  number  of  approximately  4000. 

2.  Laminar  flow  results  for  the  no  insert  condition  were  in  close  agreement  with 
the  Shah  and  London  [Ref.  9]  correlation  for  a  constant  wall  heat  flux 
boundary  condition  for  all  three  tubes  tested. 

3.  The  Hausen  [Ref.  10]  correlation  used  by  Mazzone  [Ref.  6]  significantly 
overpredicted  the  inside  thermal  resistance  and  led  to  nonsensical  results  for 
h„ 


■o* 
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4.  The  smooth  tube/twisted  tape  data  were  in  close  agreement  with  the  Hong 
and  Bergles  [Ref.  13]  correlation.  Since  this  correlation  was  used  by  Mazzone 
[Ref.  6],  his  results  for  outside  heat  transfer  coefficient  should  be  fairly 
accurate  for  this  tube. 

5.  For  the  Cu/Ni  finned  and  Korodense  tubes  with  the  twisted  tape  insert,  the 
Hong  and  Bergles  correlation  predicted  higher  inside  heat  transfer  coefficients 
than  the  experimental  data  over  the  entire  Reynolds  number  range  tested. 
This  could  have  been  due  to  a  different  twist  ratio  (for  the  Cu/Ni  finned  tube) 
and  stagnant  regions  (for  the  Korodense  tube). 

6.  Experimental  data  for  all  three  tubes  using  the  HEATEX  insert  yielded 
•  results  that  were  well  below  the  correlations  developed  and  used  by  Mazzone 

[Ref.  6]  over  the  entire  Reynolds  number  range  tested. 

7.  HEATEX  inserts  provided  the  best  inside  heat  transfer  enhancement  for  all 
three  tube.  Enhancement  ratios  of  over  a  factor  of  seven  were  found. 
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VIII.   RECOMMENDATIONS 


Based  on  the  results  of  this  study,  the  following  recommendations  are  made: 


1.  Further  refinement  of  the  correlations  to  better  account  for  different  coolant 
inlet  temperatures  should  be  made  by  optimizing  the  Prandtl  number 
exponent 

2.  These  correlations  should  be  used  to  reprocess  all  refrigerant  condensation 
data  and  should  be  used  in  future  refrigerant  condensation  experiments 
conducted  at  NPS.  Note  that  since  these  correlations  were  developed  by 
omitting  coolant  velocities  of  0.2  and  0.3  m/s,  then  the  condensation  data 
should  do  the  same. 

3.  HEATEX  inserts  provide  the  best  enhancement  of  tubeside  heat  transfer. 
Due  to  the  data  reduction  technique  used,  HEATEX  inserts  should  therefore 
be  used  in  all  tubeside  refrigerant  condensation  experiments  to  increase  the 
accuracy  of  the  outside  heat  transfer  coefficient 

4.  Transitional  flow  correlations  can  also  be  developed  for  the  no  insert 
condition  since  condensation  experiments  typically  move  into  this  flow  region. 
However,  since  the  transition  region  is  not  generally  well  developed,  any 
developed  correlations  for  this  region  will  be  useful  only  for  this  particular 
apparatus. 
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APPENDIX  A.   FLOW  METER  CALIBRATION 

The  flowmeter  used  to  measure  flow  through  the  instrumented  tube  was 
initially  calibrated  at  0°C  and  24°C  by  Mazzone  [Ref.  6].  However  an  additional 
calibration  temperature  of  -15°C  was  needed  for  this  study.  The  flowmeter  was 
calibrated  using  the  same  procedure  as  Mazzone  and  is  summarized  as  follows. 

The  outlet  of  the  flowmeter  was  disconnected  from  the  test  section  and 
connected  to  an  empty  55  gallon  drum.  The  drum  was  placed  on  a  scale  with  a 
resolution  of  0.5  lbs.  The  coolant  pump  was  started  with  the  flow  through  a  bypass. 
The  flowmeter  was  initially  opened  to  maximum  flow  and  then  throttled  back  to 


10%.  The  weight  of  fluid  in  the  drum  was  recorded  and  the  timer  started.  The  weight 
was  then  recorded  at  the  end  of  the  calibration  time.  This  procedure  was  repeated 
in  10%  increments  up  to  100%  flow.  The  calibration  run  was  then  repeated  from  the 
maximum  flow  back  to  10%  flow  in  10%  increments.  The  mass  flow  rate  was  then 
calculated  in  kg/s  based  on  an  average  of  the  two  readings.  The  data  for  the 
flowmeter  calibrations  as  a  function  of  nominal  flow  rate  are  shown  in  Table  A.l. 
A  least  squares  linear  regression  was  performed  for  the  flowmeter  at  each  of 
the  three  temperatures  for  which  the  calibration  runs  were  carried  out  The  linear 
regression  equation  takes  the  following  form: 


113 


m  =  a+bN 


(A.1) 


where  a  is  the  y-intercept  in  kg/s,  b  is  the  slope  in  kg/s/%  and  N  is  the  flowmeter 
setting  in  %.  The  results  are  shown  in  Table  A.2. 

TABLE  A.1.   MASS  FLOW  RATE  (KG/S)  AT  -15,  0  AND  24°C 


FLOW  METER  A 

NOMINAL  FLOW 
RATE 

(%) 

MASS  FLOW  RATE  (kg/s) 

-15°C 

0°C 

24°C 

10 
20 
30 
40 
50 
60 
70 
80 
90 
100 

0.0114 
0.0295 
0.0483 
0.0679 
0.0866 
TJ.1078 
0.1258 
0.1453 
0.1613 
0.1824 

0.0149 
0.0363 
0.0544 
0.0794 
0.1003 
0.1179 
0.1388 
0.1563 
0.1846 
0.2019 

0.0172 
0.0403 
0.0621 
0.0851 
0.1120 
0.1339 
0.1556 
0.1769 
0.2028 
0.2151 

TABLE  A-2.   FLOWMETER  CALIBRATION  REGRESSION  RESULTS 


FLOW  METER  A 

-15°C 

o°c 

24°C 

Slope 

0.0019 

0.00206 

0.00225 

Y-Intercept 

-0.0098 

-0.00414 

-0.00329 

In  order  to  obtain  flowrates  at  any  intermediate  temperatures,  a  simple 
interpolation  procedure  was  incorporated  into  the  data  reduction  program 
"DRPSING."  The  results  of  the  calibrations  for  this  flowmeter  are  shown  in  Figure 
A.1. 
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FLOWMETER 
CALIBRATION  @  -15,  0,  24C 
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Figure  A.l   Flowmeter  Calibration. 
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APPENDIX  B.   HEAT  LOSS  TO  AMBIENT  EXPERIMENT 

A  series  of  heat  loss  experiments  were  conducted  to  account  for  the  amount 
of  heat  lost  to  ambient  To  accomplish  this,  the  coolant  was  drained  from  the  test 
tube  and  the  inlet  and  outlet  capped  off  to  prevent  any  air  flow  through  the  tube. 
Heater  power  was  set  at  a  desired  level  and  wall  temperatures,  monitored  by  the 
DAU.  Average  wall  temperatures  were  determined  every  minute  for  the  first  15 
minutes  and  then  every  four  minutes  until  an  equilibrium  temperature  was  reached 
(after  approximately  60  minutes).  At  that  point,  any  additional  heat  added  could  be 
considered  lost  to  ambient  This  procedure  was  repeated  at  five  nominal  power 


settings,  1W,  3W,  5W,  7W,  and  9W.  All  data  were  taken  at  an  ambient  temperature 
of  22.0°C. 

Figure  B.l  is  a  plot  of  average  wall  temperature  versus  time  for  the  1W  and 
7W  power  settings.  It  can  be  seen  that  an  equilibrium  average  wall  temperature  was 
reached  after  approximately  60  minutes  for  these  two  power  settings  (as  it  was  for 
the  other  three  power  settings.)  This  graph  shows  that  if  the  average  wall 
temperature  in  the  experiments  reached  46°C  or  74°  C,  then  the  heat  lost  to  the 
atmosphere  was  1W  or  7W,  respectively. 

Table  B.l  shows  the  data  for  the  equilibrium  average  wall  temperature  as  a 
function  of  the  actual  power  setting  (qUpe).  The  data  are  shown  graphically  in  Figure 
B.2.  A  least  squares  linear  regression  was  performed  and  the  resulting  equation  was: 
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(7  .-44.15)°C 


floss 


(B.1) 


4.39  wrc 

Therefore,  for  a  given  average  wall  temperature,  the  amount  of  heat  lost  to 
ambient  could  be  determined.  Equation  B.l  was  incorporated  into  the  data  reduction 
program  "DRPSING." 


TABLE  B.1   EQUILIBRIUM  AVERAGE  WALL  TEMPERATURE  RESULTS 


NOMINAL  POWER 

ACTUAL  POWER  (W) 

EQUILIBRIUM  Tw  (°C) 

1W 

0.9412 

46.47 

3W 

3.122 

59.23 

5W 

5.126 

67.90 

6.790 

74.34 

7W 

9W 

9.560 

84.85 
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Figure  B.l  Heat  loss  to  ambient;  Average  wall  temperature  versus  time  for  1W 
andTW. 
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HEAT  LOSS  TO  AMBIENT 
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Figure  B.2  Heat  loss  to  ambient;  Average  wall  temperature  versus  heat  input. 
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APPENDIX  C.    FRICTIONAL/HEAT  LEAKAGE  EXPERIMENT 


As  coolant  flows  through  the  thermocouple  wells/mixing  box  and  along  the 
tube,  a  small  temperature  rise  results  from  frictional  dissipation.  When  the  coolant 
inlet  temperature  drops  below  ambient,  an  additional  temperature  rise  can  be 
attributed  to  heat  leakage  into  the  coolant  from  ambient.  Since  there  was  no  means 
available  to  determine  these  terms  separately,  they  were  considered  as  one  heat 
leakage  term  (qieaj- 

To  determine  this  heat  leakage  term,  the  power  to  the  heating  elements  was 
secured  and  coolant  passed  through  the  tube  being  tested  at  a  given  flowmeter 
setting.  Any  temperature  rise  detected  was  therefore  due  to  frictional  effects  or  heat 
leakage  from  ambient  This  temperature  rise  was  determined  at  flowmeter  settings 
of  15%  and  100%  for  each  tube  and  insert  condition  at  all  four  nominal  inlet 
temperatures. 

Data  showed  that  flowmeter  setting,  tube  and  insert  condition  had  no 
substantial  effect  on  the  temperature  rise  and  that  only  the  coolant  inlet  temperature 
made  any  difference.  It  can  therefore  be  concluded  that  the  heat  leakage  term  was 
primarily  due  to  heat  addition  from  ambient.  The  results  are  shown  in  Table  C.l. 
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TABLE  CI  HEAT  LEAKAGE  RESULTS 


INLET  TEMPERATURE 

TEMPERATURE  RISE 

+20°C 

o.o°c 

+  10°C 

0.07°C 

0°C 

0.08°C 

-10°C 

0.18°C 

These  corrections  were  incorporated  into  the  data  reduction  program 
"DRPSING"  by  directly  subtracting  this  temperature  rise  from  the  measured  outlet 
temperature. 
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APPENDIX  D.   POWER  CALIBRATION 

Actual  heating  element  power  (and  therefore  heat  flux)  was  based  upon  the 
applied  voltage  and  total  heating  element  resistance.  To  measure  voltage  and 
current,  the  Data  Acquisition  Unit  made  use  of  a  voltage  and  current  sensing  unit 
However,  a  power  factor  needed  to  be  applied  to  these  sensed  values  to  determine 
actual  values  required  for  the  power/heat  flux  calculation. 

Voltage  and  current  meters  were  placed  in  the  circuit  downstream  of  the  power 
source  and  prior  to  the  first  heating  element  as  shown  in  Figure  D.l.  The  voltage 
rheostat  was  varied  over  its  entire  range  and  voltage  and  current  sensed  by  the 
meters  and  in  the  DAU  were  recorded  (Table  D.l).  For  each  setting,  actual  power 
was  calculated  by: 

P  =  VI  (D.l) 

and  total  heating  element  resistance  by: 

R  =  -  =  —  (D.2) 

I        V2 

In  addition,  voltage  (PFV)  and  current  (PF,)  power  factors  were  determined  by: 
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PF   =    Vacwal  and  PF,  =   Iacwal  (D.3) 

r  CH29  lCH25 

where  the  subscripts  "CH25"  and  "CH29"  refer  to  the  DAU  sensing  unit  channel 
assignments  and  "actual"  to  the  meter  readings.  All  the  above  readings  and 
calculations  were  done  for  five  rheostat  settings  as  shown  in  Table  D.l. 

An  average  heater  resistance  (R),  current  (PF,)  and  voltage  (PFV)  power  factors 
were  then  calculated.  Heater  power  could  then  be  calculated  by  one  of  the  following 
forms  of  equation: 

P  »  (PF.XPFJK^  1^  (D.4) 

or 

p  =  QVVanffR  <D-5) 

The  current  power  factor  varied  too  much  to  be  utilized  in  a  simple  calculation. 
Therefore,  it  was  decided  to  make  use  of  equation  (D.5).  The  final  equation  to 
determine  the  average  heating  element  power  (qupe  in  Watts)  for  a  given  rheostat 
setting  was: 

<?»,»  ■  [59.66^/49.6  (D.6) 

This  equation  was  incorporated  into  the  data  reduction  program  "DRPSING." 
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TABLE  D.l  POWER  CALIBRATION  RESULTS 


POWER 

(W) 

V 
(VOLTS) 

I 

(AMPS) 

R 

(OHMS) 

Vch» 
(VOLTS) 

'cms 
(AMPS) 

PFV 

PF, 

3131 

30.7 

0.61 

5032 

03135 

0.0106 

59.78 

57.50 

76.14 

61.4 
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49.51 

1.032 

0.0430 

59.49 

28.80 

30330 

122.5 

2.48 

4939 

2.054 

0.1050 

59.63 

23.60 

J   75036 

192.4 

3.90 

4930 

3.226 

0.1746 

59.64 

2230 

113.84 

234.0 

4.76 

4931 

3.915 

0.2170 

59.77 

21.90 

POWER  CALIBRATION 
ELECTRICAL  SCHEMATIC 


HEATERS 


Figure  D.  1   Electrical  schematic  of  power  calibration  experiment. 
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APPENDIX  E.   TUBE  AND  INSERT  SPECIFICATIONS 


TABLE  E.1   TUBE  CHARACTERISTICS/SPECIFICATIONS 


PHYSICAL 
CHARACTERISTICS 

Cu 
SMOOTH 

Cu/Ni 
FINNED 

Cu/Ni 
KORODENSE 

Inside  Diameter  (D,) 

13.26mm 

10.16mm 

13.35mm 

Outside  Diameter  (D0) 

15.88mm 

14.00mm 

15.85mm 

T/C  Diameter  (Dr/C)  t/c 

14.36mm 

12.62mm 

14.53mm 

Wall  Thickness 

1.31mm 

1.22mm 

1.25mm 

Length 

Overall  Length  (L) 

152.4cm 

153.7cm 

152.7cm 

Heated  Length  (LH) 

140.6cm 

143.0cm 

142.4cm 

Heated  Finned  Length 
(Lf) 

- 

122.3cm 

- 

Smooth  Section  Length 
(Ls) 

- 

15.7cm 

- 

T/C  Location  (1) 

Pos  1 

15.2cm 

17.8cm 

15.6cm 

Pos  2 

55.9cm 

57.2cm 

56.2cm 

Pos  3 

96.5cm 

96.5cm 

96.8cm 

Pos  4 

137.2cm 

135.9cm 

134.5cm 

T/C  X/D  for  LH  (2) 

Pos  1 

6.77 

2.05 

7.91 

Pos  2 

37.41 

40.80 

38.35 

Pos  3 

68.06 

79.55 

68.79 

Pos  4 

98.71 

118.27 

96.99 

Inside  Surface  Area  for  LH 

A, 

0.0586m2 

0.0340m2  (3) 

0.0600m2 

T/C  Attachment  Method 

Silver 
Solder 

Soft  Solder 

Tach  Weld/ 
Liquid  Metal 

(1)  measured  trom  tube  inlet; 
finned  heated  length. 

(2)  measurec 

i  trom  start  of  h 

sated  length;  (3)  tor 
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TABLE  E.2    INSERT  CHARACTERISTICS/SPECIFICATIONS 


TWISTED  TAPE 

NARROW 

WIDE 

thickness,  5 

width,  w 

180°  pitch,  H 

twist  ratio,  y 

material 

0.559mm 

10.16mm 

3.05cm 

3 
brass 

0.559mm 

13.26mm 

5.30cm 

4 
brass 

HEATEX 

SMALL 

LARGE 

diameter 
material 

10.16mm 
stainless  steel 

13.26mm 
stainless  steel 
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APPENDIX  F.    UNCERTAINTY  ANALYSIS  FOR  THE  INSIDE 
HEAT  TRANSFER  COEFFICIENT 


Data  file  SMNI03  (smooth  tube,  no  insert,  +  10°C/10kW/m2),was  chosen  to 
conduct  an  uncertainty  analysis  of  the  experimental  determination  of  the  inside  mean 
heat  transfer  coefficient,  hs.  In  order  to  assess  the  effects  of  coolant  velocity  on  the 
uncertainty  in  the  measurement  of  hi(  this  uncertainty  analysis  was  performed  at 
coolant  velocities  of  0.4  and  1.4  m/s  (for  the  same  data  set).  The  calculations  below 
are  shown  for  the  0.4  m/s  coolant  velocity.  Table  F.l  summarizes  the  results  of  the 
uncertainty  analysis  for  both  0.4  m/s  and  1.4  m/s  coolant  velocities. 

The  uncertainty  analysis  utilizes  the  procedure  suggested  by.  Kline  and 
McClintock  [Ref.  24].  This  procedure  states  that  if 


R=R  (xxs2sv xn) 


then  the  uncertainty  in  R,  namely,  SR  is  given  by  the  equation: 


SR  = 


k  \ 

+ 

(«*  ] 

Sx 

dx.     2 

I    2     ) 

2 
+ + 

(«*  1 

Sx 

dx      " 

2 

I  0.5 


where  x„  is  the  measured  variable  and  S\  is  the  uncertainty  in  the  measured  variable. 

The  uncertainty  in  the  coolant  temperature  measurements  is  due  to  uncertainty 

in  the  voltage  measured  by  the  thermocouples.  For  the  thermocouples  used,  this 

value  was  assumed  to  be  4  microvolts  or  0.1°C.  The  uncertainty  in  A7  is  found  from 
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the  uncertainty  in  the  coolant  inlet  and  outlet  temperatures,  7d  and  T^,  and  can  be 
calculated  as  follows: 


.   7(14) +7(15) 
2 


Taking  the  partial  derivatives  we  obtain: 


*Tco  dTco  1 

CO  CO  A 


37(14)        37(15)        2 
The  uncertainty  in  the  measurement  of  7^  is  thus: 


dTco    = 

CO 


1  QJ 


37 


\ 


[[  37(14) 


-57 


37 


I  57(15) 


.57 


/  . 


Based  on  the  thermocouple  uncertainty,  57  =  0.1°C,  we  obtain: 


57     =  0.071  °C 

CO 


The  coolant  inlet  temperature  is  calculated  from  a  single  thermocouple.  The 
uncertainty  in  this  measurement  is  calculated  as: 


37(13) 


=  1 
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and 


ST.  = 


1  0 J 


/  \2 

57 


^  37X13) 


/  J 


=  o.rc 


or, 


Since, 


ST    =  0.1°C 


l3Tj\ 


CO 


[dlT) 


=1, 


aArj 


— i 


then  the  uncertainty  in  AT  is  calculated  from: 


SAT  =  [(57J2  +  (S7a)2fJ 


or, 


SAT  =  0.1225° C 


For  a  measured  A7  of  3.7°C,  the  uncertainty  is  therefore: 


SAT 
A7 


=  0.0332  =  3.32% 
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Coolant  velocity  is  calculated  from: 

V   -     "* 

e     pA 

For  a  coolant  velocity  0.4  m/s,  rh  was  calculated  to  be  0.05954  kg/s.  The  uncertainty 
in  the  reading  of  the  flowmeter  is  defined  as  the  scale  interpolation.  This  value 
corresponds  to  one-half  the  value  of  the  smallest  marked  increment  In  this  case,  the 
scale  interpolation  factor,  8%if  is  0.5%,  which  corresponds  to  0.00101  kg/s  resolution. 
In  addition,  a  time-wise  jitter,  5^,  of  0.5%  was  noted  during  the  runs.  Therefore,  the 
uncertainty  in  rh  is  given  by: 


Srh 


m 


(Srh 


l_v 


m 


*A  + 


(Srh^2 


m 


I  . 


OJ 


For  a  the  measured  coolant  mass  flow  rate: 


Srh 


m 


=  0.0239  =  239% 


The  uncertainty  in  the  cross-sectional  area  of  the  tubes  was  estimated  from: 


&4 


n  2 


2x 


<SD^ 


D 


where  3D;  is  given  as  0.1  mm  based  on  tolerances  supplied  from  the  manufacturer. 
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Therefore,  for  an  inside  diameter  of  13.26  mm, 


SA 


=  0.0107  =  1.07% 


Assuming  that  the  uncertainty  in  the  evaluation  of  the  fluid  properties  is  negligible, 
the  uncertainty  in  the  coolant  velocity  is  now  given  by: 


5K 


or  for  a  coolant  velocity  of  0.4  m/s, 


1  0 J 


Sm\ 


m 


'*0 


5K 


-i  =  0.0262  =  2.62% 


The  uncertainty  in  the  calculation  of  coolant  Reynolds  number  is  given  by: 


8Ree 
~Re~ 


-i  -0J 


/5K\2  i 


IV     «  I 


or  substituting  appropriate  numbers  we  find 


SRe 


Re 


1  =  0.0282  =  2.82% 


The  uncertainty  in  the  calculation  of  the  heat  transfer  rate  to  the  coolant  is 
related  to  the  uncertainty  in  the  measurements  of  coolant  mass  flow  rate  and  coolant 
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temperature  rise  where: 


54 

4 


-10 j 


(8m\\/8AT\z 
\  A7 


v  m  1 


or 


il  =  0.0409  =  4.09% 


The  uncertainty  in  the  heat  flux  is  then: 


.UJ 

2 

+ 

1 

2 

T  0 J 


where  the  uncertainty  in  As  is  given  by: 


SA 


1 0 j 


/sn2 


v*/ 


Here,  the  uncertainty  in  the  length  of  the  tube,  8L,  is  assumed  to  be  5  mm  while  8Dt, 
is  assumed  to  be  0.1  mm  (again  from  given  manufactures  tolerances).  Substitution 
of  the  appropriate  numbers  yields: 


ft* 


=  0.0147  =  1.47% 
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Finally,  since  h;  is  given  by: 


h.  = 


LMTD 


the  uncertainty  in  the  LMTD  must  be  calculated.  The  LMTD  is  given  by: 


LMTD  = 


A7 


In 


T  -T  . 
T  -T 


In  order  to  calculate  the  uncertainty  in  the  measurement  of  the  LMTD,  the  following 
equations  were  used: 


dLMTD 
dbT 


1 


In 


7-7 

wt        a 

T  -7 

Wl  CO 


dLMTD 
ST. 


^[(T^-TJ-iT^TJ} 


In 


T  -7 
T  -T 

Wl  CO 


(T^TXT^TJ 


dLMTD 


AT 


In 


T  -7 

m        a 

7-7 

Wl  CO 


(7  -7  ) 

V      wi  coJ 
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dLMTD 
dT. 


AT 


In 


T,-T_ 


11  2 


T  -T 

Wt  CO 


(T..-TJ 


The  uncertainty  in  the  LMTD  is  therefore  given  by: 


8LMTD  = 


BLMTD^T\\ 


aAT 


dLMTD 


> 

2  / 

«■- 

+ 

J 

\ 

dLMTD 
dT 


ST 


dLMTD 
dT. 


y 

ST. 

a 

J 

Before  the  uncertainty  in  LMTD  can  be  calculated,  the  uncertainty  in  the 
measured  average  inside  wall  temperature  must  be  determined.  The  average  inside 
wall  was  determined  from  a  total  of  twelve  thermocouples,  three  at  each  of  the  four 
longitudinal  positions.  The  uncertainty  in  a  single  wall  temperature  measurement  is 
due  to  the  uncertainty  in  the  voltage  measurement  of  the  thermocouple  (ST  =  0. 1°C) 
and  the  uncertainty  due  to  the  influence  caused  by  the  various  attachment  techniques 
used.  Obviously,  the  variations  seen  in  the  wall  thermocouple  readings  were  due  to 
this  latter  uncertainty  and  some  account  has  to  be  made  of  this  as  this  is  by  far  the 
dominant  uncertainty.  If  the  average  temperature  at  a  given  longitudinal  position  is 
assumed  to  be  the  required  temperature,  then  the  standard  deviation  of  the  three 
thermocouple  readings  at  that  position  could  be  assumed  to  be  this  latter  uncertainty. 
Therefore,  the  uncertainty  at  a  given  longitudinal  position  can  be  represented  by: 
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ST. ,,  =  ^BD.ST 


where  i  =  1,2,3,4  for  each  of  the  four  thermocouple  longitudinal  locations  and 


— LI    is  the  standard  deviation  at  each  of  these  locations  mentioned  ab 
dT 


ove. 


Since  the  average  inside  wall  temperature  at  each  position  increases  along  the 
length  of  the  tube,  the  average  inside  wall  temperature  for  the  entire  tube  must  be 
thought  of  as  the  average  of  the  four  longitudinal  positions'  average  temperatures 
(Tpo,(,.))-  The  average  inside  wall  temperature  is  then  determined  by: 


*p  posl  pos2  posJ  par* 

*"  4 


Taking  the  partial  derivatives  gives: 


Therefore,  the  uncertainty  in  the  average  inside  wall  temperature  measurement  is: 


0.5 


57    = 


dT 


_\ 


dT 


.ST 


posl 


posl 


(  dT 


.ST 


(  dT 


pos2 


dT 


.ST 


(  dT 


posJ 


pos3 


dT 


.ST 


pos4 


pos4 
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Considering  position  1,  the  average  inside  wall  temperature  is: 

T      r(i0)+r(ii)+:r(i2) 

lpo*  = 3 - 

For  the  three  thermocouple  readings  of 

T(10)  =  25.99°C,  T(ll)  =  25.74'C  and  T(12)  -  26.01'Q 
the  average  inside  wall  temperature  for  position  1  is 

T^  =  25.91' C 
The  standard  deviation  at  this  location  is  given  by: 

SD-  ■  *Tpo,i  -  (LV.-T^/n)0*  =  0.123° C 


Therefore, 


-=£?!-  =  4.7  x  10"3  =  0.47% 
T 


Repeating  this  for  the  other  three  locations  gives: 


6Tpas2  "  0.497° C 


a7^  "  0.634°  C 


57^  =  0.145°C 
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Therefore,  the  uncertainty  in  the  measurement  of  average  inside  wall  temperatures 


is: 


S7V  =[[(0.25)(0.123)f  +[(0.25)(0.497)f  +[(0.25)(0.d34)f  +[(0.25)(0.145)fj     =  0.207 

After    substituting    the    appropriate    values    for    dLMTD/dAT,    dLMTD/dT^ 
dLMlD/dT^  and  3LMTD/3T-J  the  uncertainty  in  the  measurement  of  LMTD  is: 


8LMTD  =  1.122°C 


and 


SLMTD        1.122 
LA/TD    "   24.13 


=  0.0465  =  4.65% 


The  uncertainty  in  the  overall  heat  transfer  coefficient,  U&  is  then  given  by: 


T. 


8LMTD\ 
LMTD 


or 


Sh. 

IT 


=  [(0.0417)2+(0.0463)2]1 


Sh. 

—L  =  0.0625  =  6.25% 

Sh. 
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8h. 

—L  =  0.0625  =  6.25% 

», 

Table  F.l  summarizes  the  results  of  the  uncertainty  analysis  for  the  smooth  copper 

tube  with  no  insert  for  coolant  velocities  of  0.4  and  1.4  m/s,  respectively  at  a  nominal 
inlet  coolant  temperature  of  +  10°C. 

The  uncertainty  analysis  neglects  any  uncertainty  in  the  physical  properties  of 
the  coolant  For  the  most  part,  this  is  probably  a  reasonable  assumption.  The 
exception  to  this  may  be  the  viscosity  due  to  its  sensitivity  to  temperature.  Hence,  the 
uncertainty  reported  here  is  likely  to  be  conservative. 

It  is  curious  that  the  relative  magnitudes  of  the  uncertainty  for  the  low  and  high 
coolant  velocities  are  nearly  the  same.  Two  different  mechanisms  dominate  in  this 
calculated  uncertainty  at  low  and  high  coolant  velocities.  In  terms  of  the  mass  flow 
rate,  the  time-wise  jitter  and  scale  interpolation  terms  are  the  same  for  both  coolant 
velocities.  However,  the  actual  mass  flow  rate  is  less  at  the  lower  flow  rates.  Hence, 
the  uncertainty  in  the  mass  flow  rate  calculation  is  more  significant  at  the  lower 
coolant  velocities.  At  higher  coolant  flow  rates,  the  uncertainty  in  the  calculation  of 
LMTD  becomes  the  dominant  term.  This  is  because  the  coolant  temperature  rise  is 
smaller  at  high  coolant  velocities  but  the  uncertainty  in  the  temperature 
measurements  remains  the  same. 
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TABLE  F.l  SUMMARY  OF  RESULTS  FOR  SAMPLE  UNCERTAINTY  ANALYSIS 


Variable 

Vc  =  0.4  m/s 

Vc  =  1.4  m/s 

ST* 

0.98% 

0.87% 

5T 

CO 

T 

CO 

0.51% 

0.56% 

5Ar 

IT 

3.32% 

10.67% 

5m 
m 

2.39% 

0.69% 

8A 

1.07% 

1.07% 

V 

c 

2.62% 

1.27% 

8Rec 
Re- 

c 

2.32% 

1.47% 

Sq_ 

4.09% 

10.69% 
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57 
7 

0.57% 

0.71% 

9 

4.17% 

10.72% 

8LMTD 
LMTD 

4.65% 

16.32% 

6.25% 

19.53% 
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APPENDIX  G.    PROGRAM  LISTING 


'SELECT  OPTION"" 

■0   TAKIN6  OATA  OR  REPROCESING  PREVIOUS  OATA*" 

■1   PURSE  PILES'" 


1000  \    file    qrpsing 

1005  l   PURPOSE   This  progran  collects  and  processes  single  instrumented  tube  da 

U. 

101S  !   CREATEO;  APRIL  I.  1992  (S.  MEMORY  ) 

1030  SEEP 

103S  PRINTER  IS  1 

1040  PRINT  USING  '4X,' 

104S  PRINT  USIN6  *6X,' 

1065  PRINT  USINS  -SX/ 

1080  PRINTER  IS  701 

1085  INPUT  leal  I 

1090  IF  leal  1-0  THEN  CALL  Main 

1095  IF  Icall-1  THEN  CALL  Purge 

1125  ENO 

1130  SUB  Main 

1135  COM  /Cc/  C(7) 

1140  COM  /Fid/  Ift 

1145  COM  /Nus/  Iin.Tsat.Qdol.Hnus.Kf ,Rhof .Hfg.Muf .Oo.Itube 

1150  COM  /Will/  0oa(4)^0ia(4),Kna<4),Iact,0root<4),Ld(3) 

1155  COM  /Uil2/  0«lta.Isat,Nsets,Hod,Cia<4),Alpa«<3) 

1150  OIM  E«f<  !7),Tp(4),T(  17 > ,Ho< 4  )  ,Qdp< 4  )  ,Uo( 4  )  ,Pc< 4  )  ,Xstar< 4 ) ,Nux<  4 ) , Nuss<  4 ) 

1170  OIM  Xxxl(4),Xxx2<4>,Xxx3< 4>.Xxx4< 4 >,Xxx5<4  ),Xxx6<4> 

1180  OIM  Nu»< 4>,Nuag<4).Ldd<4  ),Nu( 4  )  ,Xs< 4  ) 

121S  OATA  1.0,5. 172. 0.226.0. 063, S. 172 

1220  REAO  Cia<«) 

1225  OATA  0.10086091.25727.94369,-767345.9295,78025595.31 

1230  OATA  -92*7486589, 6. 97688E11. -2. 66192E13. 3. 94078E14 

1235  REAO  C<«> 

1240  OATA  0.01S875, 0.014000, 0.015850, 0.0159, 0.0 

1245  OATA  0.0132S9. 0.010160, 0.013350, 0.01346. 0.0 

12S0  OATA  386.0.42.98.42.98.21.9,0.0 

1255  OATA  0.0005588.3 

1260  REAO  0oa< • >,0ia< • >,Kna< • >,0elta. Hod   '  Hod-H/Ot 

1261  IF  Itube-0  THEN 

1263  OATA  6.765,37.41.68.06,98.71 

1254  REAO  Ld<»)    i    LENSTH/OIAMETER  RATIO  FOR  4  17C  L0C3  (SMOOTH) 

1265  L-l.406    !  Heated  length  of  smooth  tube 

1266  ENO  IF 

1267  IF  Itube-1  THEN 

1268  OATA  2.05.40.80.79.55,118.29 

1269  REAO  Ld( • >    '    FIN  LENGTH/OIAMETER  RATIO  FOR  4 

1270  L-1.223      i  Heated  finned  length  of  Cu/Ni 

1271  ENO  IF 

1273  IF  Itube-2  THEN 

1274  OATA  7.908,38.35.63.792,96.987  . 

1275  REAO  Ld<»)    '    LENGTH/OIAMETER  RATIO  FOR  4  T/C  LOCS  I KORQOENSE  ) 

1276  L-l.422    i  Heated  length  of  Korodense  tube 

1277  ENO  IF 

1283  Jset-0 

1284  Okaccot-1 

1285  SEEP 

1286  INPUT  'ENTER  MONTH,  OATE  ANO  TIME  < MM: 00: HH. MM: SS >* .OtgS 

1287  OUTPUT  709;-TO-.OtgS 
1290  SEEP 
1300  I  hard- 1 

1310  INPUT  -UANT  A  HAROCOPY  PRINTOUT  (  I  -OEF-YES , 0-N ) ' , Ihard 

1315  BEEP 

1325  Iin-1 

1345  IF  Ihard-l  THEN  PRINTER  IS  701 
1355       SEEP 


T/C  LOCS  (CU/NI  FIN) 
fin  tube 
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1360        INPUT  "GIVE  A  NAME  FOR  THE  NEW  OATA  FILE"  ,Fi leS 

1365       CREATE  BOAT  FileS,20 

1370       ASSIGN  ©File  TO  FileS 

1375       BEEP 

1380        INPUT  "ENTER  TUBE  COOE  ( 0-SMOOTH  CU.l-FINNEO  CU/NI , 2-K0R0DENSE.3-TIT 

FIN  )",Itube 

1383  BEEP 

1384  INPUT  "0-NO  INSERT,  1-TWISTED  TAPE  (WIDE),  2-HEATEX  <D,  3=HEATEX  <S 
,  4-TWISTED  TAPE  < NARROW  )", Insert 

1385  IF  Insert-l  THEN  Hod-4 

1386  IF  Insert-4  THEN  Hod-3 

1388  BEEP 

1389  INPUT  "NOMINAL  SUMP  TEMP.  0—10  C,  1-0  C,  2-+10  C,  3-+20  C,  I  temp 

1390  BEEP 

1391  Egrat-S4 

1392  INPUT  'ENTER  EG  CONCENTRATION  ( WT  PERCENT > .OEFAULT  -  54X",Egrat 

1393  BEEP 

1395  •     ENTER  709;Otg$ 

1400       OUTPUT  9File;OtgS 

1405       OUTPUT  8FUe;  Itube, Egrat,0dl,0d2,0d3,Dd4,0dS 

1430       PRINT  USING  "10X,"*FILE  NAME=    " " , 12A" ; Fi leS 

1435       PRINT 

1595  Iout-1 

1600  BEEP 

1605  INPUT  "WANT  TO  CREATE  AN  OUTPUT  FILE?  ( 1-OEF-YES ,0-N0 >' , lout 

1610  IF  Iout-l  THEN 

1615  •   BEEP 

1620        INPUT  'ENTER  A  NAME  FOR  OUTPUT  FILE'.Fout* 

1625       CREATE  BOAT  FoutS.S 

1630  ASSIGN  SFout  TO  FoutS 
1635  ENO  IF 

1640  i 

1646  IF  Itube-0  THEN  Otc-. 01436      "  DIAMETER  TO  T/C  8URIAL  (SMOOTH) 

1647  IF  Itube-1  THEN  Otc-. 0126  !  OIAMETER  TO  T/C  BURIAL  (CU/NI  FIN) 
1649  IF  Itube-2  THEN  Otc-. 01453  '  OIAMETER  TO  T/C  BURIAL  (KOROOENSE) 
1850  Oo-Ooa< Itube) 

1651  Ool-Ooa(  Itube) 

1553  Oi-Oia< Itube) 

1655  Km-Kma< Itube  ) 

1660  Ax-PI»0i"2/4   i  Cross-sectional  area 

1661  IF  Itube-2  THEN  A*-. 00014    i  BASED  ON  WOLVERINE  OESIGN  DATA 

1665  Ao-PI»Oo«L 

1666  Ai-PI»Oi«L 

1667  IF  Itube-2  THEN  Ao-.0S»L     '  BASED  ON  WOLVERINE  OESIGN  DATA 

1668  IF  Itube-2  THEN  Ai-Ao/ 1.186  '  BASED  ON  WOLVERINE  OESIGN  OATA 
1670  Penm-PI»0i 

1571  Anarea-PI»(0o'2-0i"2  )/4. 

1672  Ldd(  1  )-L/( 4. »Oi  ) 

1673  Ldd( 2  )-L/(2. »Oi  ) 

1674  Ldd( 3  )«3. »L/< 4. »Oi  ) 

1675  Ldd(4)-L/0i 

1676  IF  Itube-0  THEN 

1677  Ulenl-,0627     I  UNHEATED  INLET  ENO  OF  COPPER  SMOOTH  TU8E  ONLY 
1679  Ulen2-.0S72     i  UNHEATED  OUTLET  ENO  OF  COPPER  SMOOTH  TUBE  ONLV 
1679  ENO  IF 

1690  IF  Itube-1  THEN 

1631  'Jlenl».045     '  UNHEATED  INLET  END  OF  CU/NI  FIN  TUBE  ONLY 

1632  Ulen2=».06      '  UNHEATED  OUTLET  END  OF  CU/NI  FIN  TUBE  ONLY 
1S33  Ocuo-. 015875   !  UNHEATED  OUTSIDE  DIAMETER  OF  CU/NI  FIN  TUBE 
1634  Ocui-. 013259   '  UNHEATED  INSIDE  OIAMETER  OF  CU/NI  FIN  TUBE 


LS85 
[685 
LS87 
1588 
L689 
LS90 
1691 

:692 

1693 

.694 

L700    ! 

1701 

1702 

L703 

'a)'", 

:704 

1705 

L706 

L707 

L708 

L709 

L710 

1711 

L712 

1713 

1714 

.715 

.716 

.717 

.718 

.719 

L720 

L721 

L722 

.723 

.724 

.725 

.725 

;727 

.728 

:729 

1730 

.731 

.733 

L734   R 

.735    ! 

L737 

.738 

L739 

.740 

i741 

1742 

1744 

L74S    ' 

.746 

L747 

.748    i 

L749 

.750    i 

.751 


Lcui-.  112 

Lcuo-.097 

PeriPi-PI'Ocui 

Anarea-PI*( Ocuo~2-Ocui* 

ENO  IF 

IF  Itube-2  THEN 

Ulenl-,050     i 

Ulen2-.055     ' 

ENO  IF 


HEATED  INLET  SMOOTH  ENO  OF  CU/NI  FIN  TUBE 
HEATED  OUTLET  SMOOTH  ENO  OF  CU/NI  FIN  TUBE 


)/4. 


UNHEATED  INLET  ENO  OF  CU/NI  KORODENSE  TUBE  ONLY 
UNHEATED  OUTLET  ENO  OF  CU/NI  KOROOENSE  TUBE  ONL/ 


Hr-50 


i  HEATER  RESISTANCE 


PRINTER  IS  1 
BEEP 

INPUT  "WHAT  IS  REQUIRED  NOMINAL  COOLANT  VELOCITY  IN  n/s    (MAX 
Uegst 

OUTPUT  709; "AR  AF13  AL13  UR5" 
OUTPUT  709; 'AS  SA" 
Vsum-0 

FOR  J-l  TO  5 
ENTER  709 ;E 
Usum"Vaufn+E 
NEXT  J 

Emf< l3)-Vsum/5 
T<  13)-FNTv3v(Emf(  13)) 
Tavgl-T< 13) 
IF  Tavgl<-15  OR  Tavgl>24  THEN 

PRINT  -COOLANT  TEMP..  NOT  IN  RANGE  -15  TO  24  C" 
GOTO  1703 
END  IF 

Rhoegl-FNRhoeg( Tavgl ,Egrat  ) 
Mdotl"Rhoegl»Ax»Uegst 
PcalS-(Mdotl  +  .0098  )/.00l9 
Pca0-<Mdotl+.00414  )/. 00206 
Pca24-(Mdotl+. 00329)/. 00225 
IF  Tavgl<0  THEN 

Pc(0)-PcalS+<  <Tavgl  +  15)»(Pca0-Pcal5)  )/15 
ELSE 

Pc< 0  )-Pca0+<  Tavgl»( Pca24-Pca0  )  >/24 
END  IF 

PRINT  "SET  FLOWMETER  REAOINGS  CORRESPONDING  TO  ■  " 
PRINT  USING  *SX,"*X  OF  METER  A  -  *" .OOO.D" ;Pc< 0  ) 
PRINT 

PRINT  "HIT  CONTINUE  UHEN  REAOY"  • 
PAUSE 
epeat :     ' 

PRINTER  IS  701 
Ido-1 

ON  KEY  0, 15  RECOVER  1734 
PRINT  USING  "4X, ""SELECT  OPTION  """ 
"0-TAKE  DATA  *"" 

"1-SET  Heater  Tape  Heat  FLu.<  (DEFAULT)  

"2-CHECK  FLOWMETER  7.  SETTING  " " " 
"NOTE:  KEY  0  «  ESCAPE""" 


1.4  pi 


PRINT  USING  "SX 
PRINT  USING  "SX 
PRINT  USING  "6X 
PRINT  USING  "4X 


Ido  ■  desired  option 

BEEP 

INPUT  Ido 

Set  default  value  for  input 

IF  Ido>2  THEN  Ido-2 

Take  data  opt  ion 

IF  Ido-0  THEN  1816 
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1752 

1754  i 

1755 

1756 

1758 

1759 

17S1  i 

1762 

1764 

1765 

1766 

1767 

1768 

1769 

1770 

1771 

1772 

1773 

1784  i 

1785 

1789 

1793 

1794 

1795 

1796 

1797 

1798 

1799 

1800 

1801 

1803 

1804 

1805 

(  12  ) 

1806 

1807 

) 

1808 

1809 

) 

1810 

1811 


IF  Ido-2  THEN  1703 

Loop  to  set  heat  flux  and  check  temperatures 

IF  Ido-1  THEN 

INPUT  "ENTER  DESIRED  HEAT  FLUX", Qdpuan 

PRINT  USING  "4X,""  Qdp  (wanted)     Qdp  (actual)    Tin     Oelf"' 

PRINT  USING  "4X,""    <U/m*2)         (U/n"2)       <C)      (C) 

Read  thermocouple  voltages 

OUTPUT  709; "AR  AF1  AL15  VR5" 

FOR  1-1  TO  IS 

OUTPUT  709; "AS  SA" 
Usun"0 

FOR  J-l  TO  20 
ENTER  709 ;E 
Usun"Usun+E 
NEXT  J 

Emf< I >-Usum/20 
T<  I)-FNTv3v<Emf< I > > 
NEXT  I 
Compute  average  temperature  of  coolant  outlet 
Tout-(T( 14  )+T( 15)  )/2 
Oelt-Tout-T< 13) 
OUTPUT  709; "AR  AF25  AL25  UR5" 
OUTPUT  709; "AS  SA" 
ENTER  709; Amp 
OUTPUT  709; "AR  AF29  AL29  URS" 
OUTPUT  709; "AS  SA" 
ENTER. 709;Uolt 
Qtape-(S0»Uolt  )*2/Hr 
Qdptape"Qtape/Ao 
PRINT 

PRINT  USING  "4X,S(MZ.30E,3X )". Qdpuan, Qdptape,T< 13 >,Delt,Qtape 
PRINT 
PRINT  USING  "4X, ""TEMPS  AT  POSITION  1 " " ,3( 000  .DO , 3X  )" . T< 10  ) , T( 1 1  ) , T 

PRINT 

PRINT  USING  "4X, ""TEMPS  AT  POSITION  2" " ,3( 000 .DO, 3X  )" ; T( 7  ) , T( 8  ) , T( 9 

PRINT 

PRINT  USING  "4X, ""TEMPS  AT  POSITION  3" " ,3( ODD .DO, 3X  )" ; T< 4  ) , T( 5  ) , T< S 

PRINT 

PRINT  USING  "4X, ""TEMPS  AT  POSITION  4" " ,3( 0DD.0O,3X  )" ; T< 1  ) , T( 2 ) , T( 3 


1812  WAIT  2 
1913  GOTO  1762 
1814  ENO  IF 

1315  i  TAKE  OATA  IF  Im-0  LOOP 

1816  Tun»0 

1813  Twi-0 

1819  PRINTER  IS  701 

1320  OUTPUT  709. "AR  AF1  AL15  VR5" 

1821  FOR  1-1  TO  15 

1822  OUTPUT  709; "AS  SA" 

1823  Vsum-0 

1824  FOR  J-l  TO  20 

1825  ENTER  709. E 

1826  Usum-Usum+E 

1827  NEXT  J 

1828  Emf (  I  )-Vsum/20 

1829  NEXT  I 
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1840  I 

1845  i 

1850  i 

1860 

1865 

1870 

1871 

1872 

1873 

1874 

1876 

1880 

1881 

1883 

1884 

1885 

1886  ■ 

1887 

1888 

1830 

1891 

1893 

1894 

1895 

1896 

1897 

1898 

1923 

192S 

1926 

1930 

1935 

1940 

1941 

1942 

1943 

1945 

1946 

1948 

1949 

1951 

1952 

1953 

1954 

1955 

1956 

1960 

1970 

L990 

1996 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 


OATA  ANALYSIS 

FOR  1-1  TO  15 

T< I )-FNTvav(Emf<  I  )  ) 
NEXT  I 
FOR  1-1  TO  12 

Twi-Twi+T(  I  ) 
NEXT  I 

Twiav"Tuu/12. 
Tout-(T(  14)+T(  l5))/2 
Oelt-Tout-T(  13) 


1  Heat  leakage  correction 
THEN 

THEN  Hlc-.ia 
THEN  Hlc-.08 
THEN  HU-.07 
Hlc-0. 


THEN 


Hlc-0. 
IF  Itube-0 
IF  Itemp-0 
IF  Itemp-1 
IF  Itemp-2 
IF  Itemp-3 
Oelt-Oelt-Hlc 
ENO  IF 
Tout-Tout-Hlc 

OUTPUT  709; "AS  SA" 

ENTER  709; Amp 
OUTPUT  709; *AR  AF29  AL29  VR5" 

OUTPUT  709; "AS  SA" 

ENTER  709; Wo  It 
Qtape-<60»Volt  >*2/Hr 
Qdptape-Qtape/Ao 
Jaet-Jaet+l 

IF  Okaccpt-0  THEN  Jaet-Jaet-1 
PRINT 

PRINT  USING  "10X,""Oata  set  n 
PRINT 

IF  Itube-0  THEN  PRINT  USING  " 
IF  Itube-2  THEN  PPINT  USING  " 
IF  Itube-1  THEN  PRINT  USING  " 
IF  Itemp-0  THEN  PRINT  USING  " 
IF  Itenp-1  THEN  PRINT  USING 
IF  Itemp-2  THEN  PRINT  USING 
IF  Itemp-3  THEN  PRINT  USING 
IF  Inaert-0  THEN.  PRINT  USING 
IF  Insert-1  THEN  PRINT 
IF  Inaert-4  THEN  PRINT 
IF  Inaert-Z  THEN  PRINT 
IF  Inaert-3  THEN  PRINT 
10X 
10X 


I 


HEAT  LEAKAGE  CORRECTION  AT  -10  C 
HEAT  LEAKAGE  CORRECTION  AT  3  C 
HEAT  LEAKAGE  CORRECTION  AT  +10  C 
NO  HEAT  LEAKAGE  CORRECTION  AT  +20  C 


*,DO" ;  Jaet 


USING  " 
USING  " 
USING  •' 
USING  " 
"Nominal  set 
"Nominal  set 


Insert-1 
Insert-4 
Inaert-Z 
Inaert-3 
PRINT  USING 
PRINT  USING 
PRINT 

Tavg-T(  L3  )+Oelf  .5 

Rhoeg»FNRhoeg<  Tavg.Egrat  ) 

Nuueg-FNNueg( Tavg.Egrat  ) 

Mueg-Nuueg»Rhoeg 

Cpeg-FNCpeg< Tavg.Egrat  ) 

Keg-FNKeg( Tavg.Egrat  ) 

Pr eg -Cp eg •Mueg/Keg 

Mdot-FNFmcaK  T(  13  ),Pc<0)  ) 

Veg-Mdot/f  Rhoeg»A>.  ) 

Reeg»veg*Oi/Nuueg 

Rea-4*Mdot/(PI»Mueg»<0i-4»0elta ) ) 

Qdot=f1dot»Cpeg«Oelt 


"  10X," "Smooth  tube""" 

*  10X , "  "Korodenae  tube"" 

"  10X, ""Copper/Nickel  fin 

"  10X ," "Nominal  sump  temp 

" 10X , " "Nominal  sump  temp 

"  10X ," "Nominal  sump  temp 

"  10X ," "Nominal  sump  temp 
"10X,""No  insert  used"" 
"10X,""Uide  twiated  tape  inaert  ua 
"  10X ," "Narrow  twiated  tape  insert 

"10X. " "Large  HEATEX  insert  used 

"10X, ""Small  HEATEX  insert  used 

heat  flux   -  "" ,rtZ.30E" .Qdptape 


-10  C" 

0  C" 

+  10  C" 

+  20  C" 


ed""" 
used" 


coolant  velocity 


""  ,0.0" .  Vegat 


145 


2048'  ITERATE  TO  FIND  INSIDE  WALL  TEMPERATURES 

2049  Tnl-Twiav 

2050  Tn2»Tuiav 

2051  Tnl»(Tnl+Tm2 >/2. 

2052  Cti-Qdot»LQG<Otc/Oi >/<2.»PI*L»Km ) 

2053  Tn2-Tuiav-Cti 

2056  IF  ABS(Tn2-Tnl  >>.001  THEN  GOTO  20S1 

2057  FOR  I-l  TO  12 

2058  T(I)-T(I)-Cti 

2059  Tun-Tun+T< I ) 
20S0  NEXT  I 

2061  Tun-Tun/ 12 

2062  Nuuegu»FNNueg<  Tun.Egrat  ) 

2063  Muegw-Nuuegu*Rhoeg 

2064  Twnl-(T( 12)+T( 11  )+T< 10)  )/3 

2065  Tun2-<T(9>+T(8>  +  T<7>  >/3 

2066  Tun3-<T(6)  +  T<5>  +  T(4)  )/3 

2067  Twm4-(T<3)+T<2)+T(l  ))/3 

2068  Tunl2-<T( 12)+T( 11  )+T( 10 >+T< 9 )+T< 8 >+T< 7  )  )/6 

2069  Tunl23*<T<  12  )+T(  11  )+T<  10  )+T(  9  )+T<  8  )+T(7  )+T(  6  )+T(S)+T(4  )  )/9 

2070  Tln-<  Tout-T< 13  )  >/LOS< <  Tun-T< 13 )  )/<  Tun-Tout  )  ) 

2071  Qloss-(Tum-44. 15J/4.39         !    HEAT   LOSS    TO    ATMOSPHERE 

2072  IF    Tun<45    THEN   Qloss-0 

2073  Qtape-Qtape-Qloss 

2074  Qdpi-Qdot/Ai 

2075  Alfai-Qdpi/Tln 

2076  i     FIN  CORRECTION  TO  ACCOUNT  FOR  CONOUCTION  TO  ENOS  OF  TUBE 

2077  Ffl-<Kn«Penn«Anarea>*.5 

2078  Alfail-Alfai 

2079  AlfaiZ-Alfail 

2080  Alanl-<  <  Al  f  ai2»Perm  )/(  Kn»Anarea  )  )*  .5 

2081  Ff3«Ff l«Alfai2"  .5«( Tun-T<  13)  ) 

2082  Ff4«Ffl»Alfai2".5»<  Tun-Tout  ) 

2083  FuncS-Ff3»FNTanh( Alanl'Ulenl ) 

2084  Func7=-Ff4»FNTanh(  Alanl»Ulen2  > 

2085  Func8-Alfai2»Ai»( Tun-Tavg  ) 

2086  Func*-Func6+Func7+Func8-Qdot 

2087  0func6«<  . S»FuncS/Al f ai2  >  +  2 . »Ff 3»Ulenl/( 1 . +FNCosh<2. »A 1  an L »U1 en  1  )  ) 

2088  0func7-<  .S«Func7/Alfai2  )+2.»Ff4«Ulen2/<  1 . +FNCosh< 2 . »ALanl »Ulen2  )  ) 

2089  0func8-Ai»( Twn-Tavg ) 

2090  0funcx-0func6+0func7+0func8 

2091  Alfail-Alfai2-Func</0funcx 

2092  IF  A8S(Alfail-Alfai2  )^.0S  THEN  GOTO  2079 

2093  •  AlfaiZ-Alfail 

2095  Cfeg-(Mueg/Muegu)" . 14 

2096  '  FINO  TEMPERATURE  INCREASE  DUE  TO  UNHEATED  LENGTHS 

2097  0eltin-Func6/<Mdot»Cpeg) 

2098  0eltout=>Func7/'Mdot»Cpeg ) 

2099  Thlm=»T(  13  KOelUn 

2100  Thlout-Tout-Oeltout 

2101  Delhl-<  Thlout-Thlm  ) 

2105  IF  Itube-1  THEN 

2106  Heatin-Lcui»Oelhl/(L+Lcui+Lcuo)  I  HEATED  SMOOTH  INLET  SECTION  FOR  CU 
/NI  FIN  TUBE 

2107  Heatout«Lcuo»Oelhl/<L+Lcui+Lcuo >  '  HEATED  SMOOTH  OUTLET  SECTION  FOR 
CU/NI  FIN  TUBE 

2109  Thlm-Thlin+Heatm 

2110  Thiout-Thlout-Heatout 

2111  Oelhl»< Thlout-Thlin  ) 
2114  ENO  IF 
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21  IS 

2116 
2117 
2118 
2119 
2120 
2121 
2122 
2123 
2124 
212S 
2127 
2128 
2129 
2130 
2131 
2132 
2133 
2134 
213S 
2136 
2137 
2138 
2139 
2140 
2141 
2143 
_2144_ 
2145 
2146 
2147 
2148 
2149 
2150 
2151 
2152 
21S3 
2154 
21S5 
2156 
2157 
2158 
2159 
2151 
2162 

2163  ! 

2164  i 

2165  ' 
2166 
2167 
2168 
2169 
2170 
2171 
2172 
2173 
2175  ! 
2176 
2177 
2178 


Qhl«Mdot»Cpeg»Oelhl 
Qdphl-Qhl/Ai 
Qdp*x-Func8/Ai 
Tcol-Thlin+<0elhl/4 ) 
Tco2-Thlm+<0elhl/2> 
Tco3«Thlin+<  3*0elhl/4 ) 
Ttcl»T(13)+<Ld(0)»Oi«Oelhl/L) 
Ttc2-T< 13>+<Ld< 1 >«0i»0elhl/L ) 
Ttc3-T(  13)  +  <Ld<2)»0i»0elhl/L> 
Ttc4-T<  l3)+<Ld(3>*0i«0elhl/L) 

Tlmall«<  Thlout-Thlin  )/L06( <  Twn-Thlin )/<  Tum-Thlout ) ) 
Tlnl-<Tcol-Thlin)/LOG( ( Twnl-Thl in )/(  Tuml-Tcol ) ) 
Tlr»2-<  Tco2-Tcol  )/L0G(  <  Twm2-Tcol  )/<  Tum2-Tco2  >  ) 
Tlw3-<  Tco3-Tco2  )/L0G< ( Tum3-Tco2  >/<  Tum3-Tco3  )  ) 
T1m4-( Thlout-Tco3 )/L0G< ( Twm4-Tco3  )/<  Twn4-Thlout ) ) 
Tlml2-<Tco2-Thlin )/L0G( < Twml2-Thl in )/( Twml2-Tco2 ) ) 
Tlml23-<  Tco3-Thlin )/L0G< ( Tunl23-Thl in )/( Tuml23-Tco3 
Alpall-Qdphl/Tlmall 
Alpl-Qdphl/Tlml 
Alp2-Qdphl/Tlm2 
Alp3-Qdphl/Tln3 
Alp4-Qdphl/Tlm4 
Allocl-Qdphl/(Twnl-Ttcl ) 
Alloc2»Qdphl/< Twm2-Ttc2  ) 
Alloc3-Qdphl/( Twn3-Ttc3 ) 
Alloc4-Qdphl/(Twm4-Ttc4  ) 
Alpl2-Qdphl/Tlml2 
Alpl23-Qdphl/Tlml23 
Kegl-FNKeg( ( Tcol+Thl in )/2 , Egrat ) 
Keg2*FNKeg<  < Tco2+Tcol )/2, Egrat ) 
Keg3-FNKeg<  (Tco3+Tco2  )/2, Egrat  ) 
Keg4-FNKeg< ( Thlout*TCo3  )/2, Egrat  ) 
Kegl2-FNKeg(  (  Tco2+Thlin  )/2, Egrat  ) 
Kegl23-FNKeg( ( Tco3+Thl in  )/2 , Egrat ) 
Numall-Alpall'Qi/Keg 
Nul-nlpl»Oi/Kegl 
Nu2«Alp2»0i/Keg2 
Nu3-Alp3»0i/Keg3 
Nu4-Alp4»0i/Keg4 
Nulocl-Allocl*Oi/Kegl 
Nuloc2-Alloc2»0i/Keg2 
Nuloc3-Alloc3»0i/Keg3 
Nuloc4-Alloc4»0i/Keg4 
Nul2-Aipl2*0i/Kegl2 
Nul23-Alpl23»0i/Kegl23 
LOCAL  NUSSELT  FROM  SHAH/L0N00N  PAGE  128,  EQN  240-242 
CALCULATE  XSTAiR  FIRST  FOR  EACH  THERMOCOUPLE  LOCATION 
OEFINE  COOLANT  PECLET  NUM8ER 
Peag»Preg»Reeg 

Xs<  1  )-Ldd< I  >/Peeg   !  FOR  MEAN  VALUES  OF  h 
X3(2  >-Ldd(2 )/Peeg 
Xa<3)-Ldd<3>/Peeg 
Xs(4  )-Ldd(4  >/Peeg 
FOR  I- 1  TO  4 
Xstar<  I  )=>Ld<  1-1  >/Peeg  !  FOR  LOCAL  VALUES  AT 

Nua( I  )-l .302 »X star <  I  )*(-l/3)-.S 
DEVELOPING  HYORO  AND  THERMAL  BOUNOARY  LAYER 

Xxl-(  1  +  U220/PI  >«Xstar<  I  ))*(-10/9)  )"<3/S) 

Xx2-(  l+(Preg/.0207  >"<2/3  >  >*  .5 

Xx3-PI/( llS.2»Xstar(  I  )  ) 


ALL  (1,2,3,4) 
SECTION  1 
SECTION  2 
SECTION  3 
SECTION  4 
'  SECTION  1,2 
)!  SECTION  1,2.: 


T/C  POSITIONS 
(S/L  P148,  EQN. 


ZS3  I 


14: 


2179  Xx4-< l+(Xx3/(Xx2«Xxl  )  r<  5/3  >>*<  3/10  > 

2180  XxS>< l+< (2Z0/PI  >«Xstar< I  )  r<-10/9> r<3/10) 

2181  Nus3<  I  )-(X.<4*S.364»XxS)-l 

2188  IF  Xstar(  I  K.00005  THEN 

2189  Nux( I  )»1.302»Xstar<  I  >*<-l/3)-l 

2190  END  IF 

2194  IF  Xstar( I  )>.001S  THEN 

2195  Nux< I  )-4. 364  +  8. 58*  <  l000»Xstar< I )  >"(-.S06 >»EXP< -41»Xstar(I)) 

2196  ENO  IF 
2200  NEXT  I 

2202  IF  Inaert-0  THEN 

2203  PRINT  USIN6  "3X,""  LOCAL  NUSSELT  VALUES""" 

2205  PRINT  USING  "3X,""    Exp.  Nux(TLfl)    Exp.  NuxOULK)    Pred.  Nux(  1  ) 
Pred.  Nux(2)    Xstar<tc>""" 

2206  PRINT 

2207  PRINT  USING  "3X,"*P0S  1  " "  ,5( MZ . 3DE.5X  )" ,Nul .Nulocl ,Nux< 1  )  ,Nuis( 1  ) , X 
star< 1 > 

2208  PRINT  USING  "3X,""P0S  2  " "  ,5( MZ.30E.SX  )" ;Nu2,Nuloc2 ,Nux( 2  )  ,Nuss< 2  )  ,X 
star(2> 

2209  PRINT  USING  "3X,"*P0S  3  " "  ,5< MZ .3DE.5X  )" ;Nu3,Nuloc3,Nux( 3  )  ,Nuss< 3  )  ,X 
star<3) 

2210  PRINT  USING  "3X,"*P0S  4  " "  ,5( MZ .30E.SX  )" ;Nu4,Nuloc4,Nux< 4 >  ,Nuss< 4  ) , X 
star(4) 

2211  PRINT 

2215  ELSE 

2216  PRINT  USING  "3X ," 'POSITION  t    •«    Exp.  Nux  < TLM )     Exp.  Nux  (BULK) 
Xstar  at  tc  pos. " " " 

2217  PRINT 

2218  PRINT  USING  "3X  ,"  "POSITION  1  (LOCAL)  -  " " ,  3<  MZ  .3DE,,5X  )"  ;  Nul  ,Nuloc  1 ,  X 
start  1  ) 

2219  PRINT  USING  *3X ," "POSITION  2  (LOCAL)  -  " " ,3( MZ .30E.5X  )* ; Nu2,Nuloc2 , X 
star(2  ) 

2220  PRINT  USING  "3X ," "POSITION  3  (LOCAL)  -  " " ,3< MZ .30E , 5X  )" ; Nu3 ,Nuloc3 ,X 
star(3) 

2221  PRINT  USING  "3X ," "POSITION  4  (LOCAL)  -  " " , 3( MZ .3DE.5X  ) " , Nu4  ,Nuloc4  ,  X 
star(4) 

2222  ENO  IF 

2223  PRINT 

2230  FOR  1-1  TO  4 

2Z3Z  IF  Insert-0  THEN 

2233  Bbl».0668»Reeg»Preg/Ldd< I  ) 

2234  8b2=«l  +  .04.(Reeg»Preg/Ldd<  I  )  r  .5665  . 

2235  Nueg< I >-3.66+<Bbl/Bb2 ) 

2236  ENO  IF 

2237  IF  Insert-1  OR  Insert-4  THEN 

2239  Nueg< I >-Cia<  Insert  )•<  1>5. 484E-3»Preg" ,7»( Res/Hod  >* 1 .25  )* .5 

2240  .Nuoth-.383»Preg-.3S»< Res/Hod  )" .622 

2241  ENO  IF 

2242  IF  Insert-2  THEN 

2243  Nueg( I  >-Cia< Insert  )»(Reeg  .65  )»Pr eg '.46 

2244  ENO  IF 

2245  IF    Insert-3    THEN 

2246  Nueg< I )-Cia< Insert  >»<Reeg" .76  >»Preg~ .46 

2247  END  IF 

2248  8EEP 

2249t  CALCULATED  VALUE  OF  INSIDE  HT .  TRANSFER  COEFF .  FROM  CORRELATIONS 

2250  Hi«Nueg< 4  )»Keg/0i 

2251'  MEAN  NUSSELT  NUMBER  FROM  SHAH/LONOON  PAGE  128,  EQN  245-246 

2254  IF  Xs<  I  X.03  THEN 

2255  Num( I  )-L.953»X3( I  )"( -1/3  ) 
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2256  ELSE 

2258  Num< I >-4.364  +  <  .0722/Xs( I >  ) 

2259  ENO  IF 
22S2        NEXT  I 

2263        PRINT  USING  "3X ," "POSITION  *   -       Exp.  nean  Nu     Pred.  mean  Nu 

Xstar  at  section""" 
22S4        PRINT 

22SG        PRINT  USING  "3X ," "POSITION  1  (MEAN)  -     " "  ,3(  MZ  .  30E.SX  ) " . Nul ,Num<  L  ) 
,Xs<  1  ) 

2257  PRINT  USIN6  "3X ," "POSITION  12  (MEAN)  -    " " ,3< MZ .30E.SX  )" . Nul2 ,Num(  2 
>,Xs<2> 

2258  PRINT  USING  *3X ," "POSITION  123  (MEAN)  -   " "  ,3< MZ .30E.5X  )" . Nul23 ,Nun< 
3>,Xs<3> 

2269  PRINT  USING  "3X ," "POSITION  1234  (MEAN)  -  " " ,3( MZ . 30E.SX  ) " . Numall ,Num 
<4),Xs(4) 

2270  PRINT 

2271  PRINT 

2272  PRINT  USIN6  "10X,""Mass  flow  rate  (kg/s>    -  " "  ,MZ .  30E" .Mdot 
2273'  PRINT  USING  ' 10X ,"" Inside  Tube  Oia.  <m>    -  " " ,MZ .30E" . Oia( I  tube > 

2274  IF  Insert-1  THEN  PRINT  USING  *10X,"180  OEG  OVER  Oia.  (HOD)  -  "", 
MZ.30E";Hod 

2275  PRINT  USING  "10X, ""Inlet  temperature  (C)    -  " "  ,  MZ  . 30E" : T(  13  > 

2276  IF  Insert-1  OR  Insert-4  THEN 

2277  PRINT  USING  'lOX.-'OELTA  Tape  Thickness    -  " " ,MZ . 3DE" ; Oelta 

2278  ENO  IF 

2279  PRINT  USIN6  "10X,"'OELT   temp   Dif.  <C>     -  " " ,MZ . 30E" ; Oelt 

2280  PRINT  USING  "10X,""H.  T.  to  coolant  (U)     -  " "  ,MZ . 30E" ; Qdot 

2281  PRINT  USING  ,18X."H.  T.  fron  heater  tape  (U)  -  " "  ,MZ .  30E"  .Qtape 

2282  PRINT  USING  "  10X ," "Prandt 1  number  -  " " ,MZ . 30E" .Preg    — 

2283  PRINT  USING  "  10X ," "Reynolds  number         -  " " ,MZ . 3DE" ; Reeg 

2284  IF  Insert-1  OR  Insert-4  THEN-PRINT  USING  " 10X ," "Reynolds  number  H 
IB   S    -  "".MZ-SOE-iRes 

2285  PRINT  USING  "l0X,""Hi  from  correl.  (U/m"2.K)  -  " " ,MZ .3DE" ;Hi 

2286  PRINT  USING  "10X,""Hi  from  heat  bal.  (U/m"2.K)  =  " " .MZ.3DE" ; Al fai 
2 

2287  IF  Insert-0  THEN  PRINT  USING  "10X,""Nu  for  no  insert  ( Hausen  ) 
",MZ.30E";Nueg<4) 

2288  IF  Insert-1  OR  Insert-4  THEN  PRINT  USING  "10X,""Nu  for  twisted  ta 
pe  (Bergles2  )" " ,MZ .30E" .Nuoth 

2289  IF  Insert-l  OR  Insert-4  THEN  PRINT  USING  "10X,""Nu  for  twisted  ta 
pe  (8ergles3)"",MZ.30E" .Nueg(4) 

2290  IF  Insert-2  THEN  PRINT  USING  "10X,""Nu  for  large  heate.<  (Mazzone) 
""  ,MZ.3DE" .Nueg<  4 ) 

2291  IF  Insert-3  THEN  PRINT  USING  "10X,'-,,Nu  for  snail  Heate*  (Ma==one) 
"" .MZ.30E" .Nueg(4) 

2292  IF  Insert-0  THEN  PRINT  USING  "10X,""Nu  from  Shah/London     =  "", 
MZ.3DE" ;Num( 4  ) 

2293  PRINT  USING  "  10X , " "X5 tar  for  heated  length  =  " " ,MZ . 30E" ; Xs(  4 ). 

2294  PRINT  USING  "  10X  ,"  "Experimental  Nusselt  no.-  " " ,MZ . 3DE" ;Nunal 1 

2295  PRINT 

2296  PRINT  USING  "3X , " "Uai 1  t/c  12,11,10  (pos.  L)  <C)  -  '" , 3< ODD .00 , 3X 
)"  .T< 12 ),T( 11  ),T( 10  ) 

2297  PRINT  USING  "3X,',"Uall  t/c  9,8,7     (pos.  2)  <C>  -  " "  .  3<  000 .00 ,  3X 
)'•  ;T<  9  ),T(  8  ),T(7) 

2298  PRINT  USING  "3X," "Uall  t/c  6,5,4     (pos.  3)  (C)  ■  " "  ,3( 000 .00 , 3X 
>"  .  T(S  i  ,T(5  ),T(  4  ) 

2299  PRINT  USING  "3X , " "Ual 1  t/c  3,2,1     (pos.  4)  (C)  -  ""  ,3(000. DD.3X 
)"  ; T( 3  ),T(  2 ),T<  1 ) 

2300  PRINT  USING  "3X,""Av.  inside  wall  temp  (C)  ■  " " , MZ . 20E" . Turn 

2301  PRINT 

149 


■MEAN   VALUES   AT 
'EXP.    LOG(Nu) 


'COR.    LOG(Nu)   <S/L> 


2302  Nu( I >»Nul 

2303  Nu(2)-Nul2 

2304  Nu<3>-Nul23 

2305  Nu<4  )«Numall 

2307  IF  Insert-0  THEN 

2308  FOR  1-1  TO  4 

2310  Xx/K  I  >«LOG(Nu<  I  )  > 

2311  Xxx2(I)-LGG<Xs<I>) 

2312  X*x3<  I  )-LOG<  Nueg<  I  >  ) 

2313  Xxx4< I  )-LOG(Nun< I ) > 

2314  NEXT  I 

2315  PRINT  USING  *3X,' 
2       POS  123      POS  1234*" 

2317  PRINT  USING  *SX,' 
xxl<2),Xxxl(3),Xxxl(4) 

2318  PRINT  USING  "6X , " "COR.  LOG(Nu)  (H) 
xx3<2 ),Xxx3(3),Xxx3(4) 

2313  PRINT  USING  "SX, 

xx4(2),Xxx4(3),Xxx4(4) 

2320  PRINT  USING  *6X , " "LOG( Xstar > 
xx2<2  ),Xxx2(3>,Xxx2(4) 

2321  ENO  IF 

2323  IF  Insert-2  OR  Inaert-3  THEN 

2324  FOR  1-1  TO  4 

2325  XxxK  I  )-LOG(Nu<  I  )/Preg".4S> 
232S  Xxx2(I)-L0G<Reeg> 

2327  Xxx3< I  >«LOG<  Nueg( I  )/Preg*  .45  ) 

2328  NEXT    I 

2329  PRINT  USING  "3X, 
2       POS  123       POS  1234" 

2331  PRINT  USING  *SX, 
Xxxl(2  ),X*xl<  3  ),Xxxl( 4) 

2332  PRINT  USING  "SX, 
Xxx3<2  ),Xx.<3(3  ),Xx/3(4  ) 

2334  PRINT  USING  "SX, 
Xxx2<2  ),Xxx2(3),Xxx2(4) 

2335  ENO  IF 
233S  IF  Insert-1  OR  Inaert-4  THEN 

2337  FOR  1-1  TO  4 

2338  X,.:<1(  I  >-LOG(Nu<  I  >/Preg~.35> 

2339  Xxx2( I )-LOG( Res/Hod) 

2340  X*x3< I >-LOG<Nuoth/Preg" .35  ) 

2341  XxxK  I  )-LOG<Nu(  I  )  ) 
2342'  X*x5< I >-LOG< Preg»<  Res/Hod  )"  1 .78  ) 

2343  XxxS< I  )-LOG<Nueg<  I  )  ) 

2344  NEXT    I 

2345  PRINT  USING  '3X,  *"I1EAN  VALUES  .AT 
Z        POS  123       POS  1234 

2346  PRINT  USING  "6X 
Xxxl(2  ),Xxxl(3  ),XxaK4> 

2347  PRINT  USING  "SX 
Xxx3(2  ),X/.x3<3  >,Xxx3<  4  ) 

2348  PRINT  USING  "SX 
Xxx2(2  >,Xxx2(3),Xxx2(4  ) 

2349  PRINT 

2350  PRINT  USING  "SX 
Xxx4(  2  )  ,Xxx4<  3  )  ,Xxx4<  4  ) 

2351  PRINT  USING  "SX 
Xxx6(2  >,XxxS<3  ),XxxS< 4  ) 

2352  PRINT  USING  "SX 


"MEAN  VALUES  AT 
"EXP.  L0G(Nu/Pr*.4S > 
"COR.  LQG<Nu/Pr".4G) 
"LQG(Re) 


"EXP.  LQG(Nu/Pr' .35  ) 
"COR.  LOG(Nu/Pr ' .35) 
"LOG<Re3/y> 

"EXP.    LOG(Nu) 
"COR.    LOG(Nu) 


POS    1  POS    1 

,4(MZ.3DE,4X  )";Xxxl(  1  ),X 
,4(MZ.30E,4X  )" ;Xxx3< 1  ),X 
,4<MZ.30E,4X  )" .Xxx4< 1  ),X 
,4<MZ.30E,4X  )" ;Xxx2( 1  ),X 


POS    1  POS    1 

"  ,4(MZ.3DE,4X  )" ;Xxxl< 1  ), 
"  ,4(MZ.30E,4X  )".X.<x3(  1  ), 
" ,4(MZ.30E,4X  P;Xxx2< 1  ), 


"LOG<Pr<Re3/y  )"1.78)    - 
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POS    1  POS    1 

,4<  MZ.30E.4X  )"  .  Xx.-.  1(1), 
,4(MZ.30E,4X  )"  ,Xa.<3(  1  ), 
,4(MZ.30E,4X  )" ;Xxx2(  1  >, 

,4<MZ.30E,4X  )"  .Xa.-.4(  1  7, 
,4(MZ.30E,4X  )"  ;Xx.<.6( I  ), 
,4(MZ.30E,4X  )" ;XxxS<  1  >, 


Xx*5<2  >,Xxx5(3>,Xxx5U) 

2353  ENO  IF 

2354  IF  Inaert-O  ANO  Reeg>2300  THEN 

2355  Nutu-.027»Reeg".8»Preg" .3333'Cfeg 

235S  PRINT  USING  "6X ," "TURBULENT  S/T  CORRELATION  -  " "  ,  4< MZ . 30E , 4X  ) " . Nu 

tu 

2357  Xxx7-L0G<Nutu) 

2358  PRINT  USING  "SX,— TURB.  COR.  LOG<Nu)    -  " " ,4<MZ .30E.4X  )" ;Xxx7 

2359  ENO  IF 

2360  Okaccpt-1 

2361  BEEP 

2362  INPUT  "OK  TO  ACCEPT  THIS  SET  <  1-OEFAULT-YES,  0-NO  )">"  .Okaccpt 

2363  IF  Okaccpt-1  THEN  OUTPUT  8File;Pc< • >,Emf < • ), Tp< • ) 

2364  IF  In-0  THEN 

2365  Okrpt-1 

2366  BEEP 

2367  INPUT  "WILL  THERE  8E  ANOTHER  RUN  (  1-YES-OEFAULT ,0-NO >" , Okrpi 

2368  IF  Okrpt-1  THEN  1701 

2369  ELSE 

2370  IF  Jaet<Nseta  THEN  1701 

2371  ENO  IF 

2372  ASSIGN  BFile  TO  • 

2373  IF  Iout-1  THEN  ASSIGN  SFile  TO  • 

2374  SUBENO 

2375  OEF  FNGrad(T) 

2376  Srad— 3:8778S7E-S-2»4.71428S7E-8»T 

2377  RETURN  Grad 

2378  FNENO 

2379  OEF  FNKcu<T> 

2380  '   OFHC  COPPER  250  TO  300  K 

2381  Tk-T+273.1S 

2382  K»434-.112»Tk 

2383  RETURN  K 

2384  FNENO 

2385  OEF  FNNueg(Tc,Egr ) 

2386  '   RANGE  OF  UALIOITY;   -20  TO  20  DEG  C 

2387  Tk-Tc+273.15 

2388  Nul-7.1195S07E-3-fk»(7.4863347E-5-Tk»(2.6294943E-7-Tk»3.0833329E-10) ) 
2390  Nu2-4.9237638E-3-Tk»<4.9213912E-5-Tk»< 1 .6437534£-7-Tk »1 . 3333331 E- 10 ) ) 
2395  Nu3-8.S586293E-3-Tk»( 8. 8837902E-5-Tk •( 3.0495032E-7-U  »3 .  4999996E-10 ) ) 
2400  A2-<Nu3-2»Nu2+Nul )/200 

2405  Al-(Nu2-Nul-940*A2  )/10 

2410  A0-Nul-42»A1-17S4»A2 

2415  Nu-A0+Egr*< Al+Egr»A2 > 

2420  RETURN  Nu 

2425  FNENO 

2430  OEF  FNCpegf  Tc . Egr  ) 

2435  '   RANGE  OF  VALIDITY   0  TO  20  OEG  C 

2440  Tk=Tc+273. 15 

2445  Cpl-l.S70l550E+3+Tk»6.3 

2450  Cp2-1.4748125E+3+Tk.6.25 

2455  Cp3-9.5800S00E+2+Tk»7.3 

2460  A2-<Cp3-2»Cp2+Cpl )/200 

2465  Al-<Cp2-Cpl-900»A2>/10 

2470  A0-Cpl-40»A1-1600»A2 

2475  Cp="A0+E3r«<  Al+Egr»A2  ) 

2480  RETURN  Cp 

2485  FNENO 

2490  OEF    FNRhoeg( T.Egr  ) 

2495  Rol-1 .0607093E+3-T»(3.7031283E-l+T.4.0837183E-3 ) 
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2500  Ro2-1.0748272E+3-T*<4.426619SE-l+T.4„0939706E-3) 

2S0S  Ro3-l.0885934E+3-T»<S.73S5G53E-H-T«6.l28U05E-3> 

2510  A2»<Ro3-2»Ro2+Rol >/200 

2S15  Al»(Ro2-Rol-900«A2 )/t0 

2520  A0-Rol-40»A1-1600»A2 

2525  Ro-A0+Egr«< Al+Egr*A2 > 

2530  RETURN  Ro 

2535  FNENO 

2540  OEF  FNPreg<T,Egr > 

2545  Pr- FNCpeg<T,Egr  )«FNNueg( T.Egr  >»FNRoeg< T,Egr  )/FNKeg< T.Egr > 

2550  RETURN  Pr 

2555  FNEND 

2560  OEF  FNKeg<Tc,Egr) 

2SSS  '   RANGE  OF  VALIDITY:   -20  TO  20  DEG  C 

2570  TV-Tc+273.1S 

2575  Kl-2.2824708E-l+Tk.<5.S989286E-4+Tk»3.S71428SE-7) 

2580  K2-2.58466lSE-l+Tk»(2.397857lE-4+Tk*7.U28571E-7) 

2585  K3-3.2l38932E-l-Tk»(3.0042857E-4-Tk»1.428S714F.-6)  ' 

2590  A2«<K3-2«K2+K1 )/200 

2595  A1-<K2-K1-900»A2)/10 

2600  A0=»K1-40»A1-1600»A2 

2605  K-A0+Egr»<Al+Egr»A2> 

2610  RETURN  K 

2615  FNENO 

2620  OEF  FNTanMX  ) 

2625  P-EXP<X) 

2630  Q-l/P 

2635  Tanh-(P-Q)/(P+Q) 

2540 RETURN  T«anh 

2645  FNENO 

2G50  OEF  FNTvav(U) 

2655  COM  /Cc/  C<7) 

2660  T-C<0) 

2665  FOR  I-l  TO  7 

2570       .  T-T+C(  I  )HTI 

2675  NEXT  I 

2630  RETURN  T 

2685  FNENO- 

2690  OEF  FN8eta(T) 

2595  Rop-FNRho(T+. 1 > 

2700  Rom»FNRho<  T-. 1 ) 

2705  Beta»-2/( Rop+Ron  )•<  Rop-Ron )/ .2 

2710  RETURN  Beta 

2715  FNENO 

2720  OEF  FNPsat(Tc  ) 

2725  i   0  TO  90  deg  F  CURVE  FIT  OF  Paat 

2730  Tf-1.8*Tc+32 

2735  Pa-5.94S52-S+Tf«<  .  1 5352082 +  Tf  •(  l  .4840963E-3+Tf  *9.615067lE-6  )  ) 

2740  Pg=Pa-14.7 

2745  IF  Pg/0  THEN       !  +-PSIG,  —  in  Hg 

2750       Paat-Pg 

2755  ELSE 

2750       Psat-Pg»29.92/14.7 

2765  ENO  IF 

2770  RETURN  Psat 

277.5  FNENO 

2905  SU8  Purge 

2910  BEEP 

2915  INPUT  "ENTER  FILE  NAME  TO  BE  DELETED"  ,Fi leS 

2920  PURGE  Files 
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2925 

2930 

2935 

2940 

2945 

2950 

2955 

29S0 

29BS 

2970 

2975 

2976 

2977 

2980 

2985 

2990 

2995 

3000 

3005 

3010 

3015 

3020 

3025 

3030  i 

3035  ' 

3040 

3045 

3046 

3050 

3051 

3052 

30S3 

3054 

3055 

3056 

3057 

3058 

3059 

3060 

3062 

3063 

3065 

3070 

3075 

3080 

3085 

3090  ' 

3095 

3100 

3105 

3110 

3115 

3120 

3125 

3130 

3135 

3U0 

3145 

3150 

3155 


GOTO  2910 

SU8EN0 

SUB  Wilson 

COM  /Uill/  0aa<4 ),Oia<4 > ,Kma< 4  ) ,  Iact ,0root<4  ),Ld<3  ) 

COM  /Uil2/  Oelta, Isat ,Nset s ,Hod,Cia< 4  ),Alpaa< 3) 

DIM  Emf( 17  ),T( 17),Xa< 100  >,Ya(  100  >,Pc<4  ),Tp(4 ) 

BEEP 

INPUT  "PLEASE  RE-ENTER  NAME  OF  FILE 

ASSISN  9File  TO  FileS 

INPUT  "0-SMOOTH  COPPER  TUBE,  1-CUNI 

BEEP 

0-NO  INSERT,  1-  TUISTED  TAPE 


INPUT 

BEEP 

INPUT  "GIVE  A  NAME  FOR  XY  FILE'.XyS 

CREATE  BOAT  Xy$,20 

ASSIGN  9Xy  TO  XyS 

L-1.2192 

Oo-Ooa< Itube ) 

Di-Oia< Itube  ) 

Kpi-Kma( Itube  ) 

Ax-PI*0i*2/4   i  Cross-sectional  area 

Ao-PI»Oo»L 

Rm-0o*L0G(0o/0i  )/<2»Km) 

Initial  values 

Tf-Tsat 

Alpa-.655 

IF  Insert-0  THEN 

Cx-i.a. 

ENO  IF 

IF  Insert-l  THEN 

Ci-5.172 

ENO  IF 

IF  Insert-2  THEN 

Ci-.22S 

ENO  IF 

IF  Insert-3  THEN 

CI-.063 

ENO  IF 

G-9.31 

Ibeg-0 

Iend-0 

IF  lact<!0  THEN 

Ibeg-Iact 

Iend*Iact 
ENO  IF 


FOR  I«Ibeg  TO  lend 

Sa»0 

Sy»0 
SaS=0 

5.".y-0 

Jset*0 

ASSIGN  9File  TO  FileS 

ENTER  9File.0tg$, I  tube .Egrat , Odl ,0d2 ,0d3 ,0d4 ,0dS 

ENTER  9File;Pc(*),Emf(*),Tp(*) 

FOR  J»0  TO  17 

T(  J  )=>FNTvsv(Enf  (  J  )  ) 
NEXT    J 
Tvap»<  T<0  >+T(  1  )  +  T(2  )  >/3 


, FileS 

FIN,  2-KOROOENSE,  3-FIN  TITAN", Itube 

INSERT,  2-HITRAN1,  3-HITRAN2" , Insert 


'CHANGE  TO  4,  IF  FIVE  TUBES  IN  BUNDLE 
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31E0 
3151 
31S5 
3170 
3175 
3180 
3185 
3190 
3191 
3195 
3200 
3201 
3202 
3204 
3205 
320S 
3210 
3215 
3220 
3225 
3230 
3235 
3240 
3245 
3250 
3255 
32S0 
32SS 
3270 
3275 
3280 
3285 
328S 
3287 
3298 
3289 
3290 
3291 
3293 
3294 
3295 
3296 
3297 
3298 
3299 
3300 
3302 
3303 
3335 
3310 
3315 
3320 
3325 
3330 
3335 
3340 
3345 
3350 
3355 
33S0 


Tliq-(T(3)+T<4)  )/2 

Tliql-T(5) 

IF  Iaat-0  THEN 

Tsat-Tliql 
ELSE 

Taat»Tvap 
END  IF 

Grad«FNGrad<T< 1+6 ) > 
Tout-<  T<2»I+10)+T<2»I+ll ) )/2 
Oelt«Tout-T< 1+6) 
Tavg-T< 1+6  )+Oelt •  .5 
IF  I-Iend  THEN 
Oelt-T( 10)-T<9) 
Tavg-<T(  10)-T(  9)  )/2 
ENO  IF 

Water/Ethy lene  Glycol  Mixture  Properties 

Rhoeg-FNRhoeg<  Tavg.Egrat ) 

Nueg»FNNueg<  Tavg.Egrat ) 

Mueg»Nueg«Rhoeg 

Cpeg-FNCpeg<  Tavg.Egrat ) 

Keg-FNKeg(  Tavg.Egrat ) 

Pr eg-Cpeg*Mueg/K eg 

Mdot-FNFncaK  I.T(  1+5  ),Pc(  I  )  ) 

Ueg-Mdot/(Rhoeg«Ax ) 

Reeg-yeg»Oi/Nueg 

Res-4*Mdot/(Mueg«<PI»0i-4»0elta ) ) 

Qdot-Mdot»Cpeg»Oelt 

Qdo-Qdot/Ao 

Lmtd«0elt/L0G< (Taat-T(  1+6  )  >/(Tsat-T(  1+6 )-Delt ) ) 

Uo-Qdp/Lmtd 

IF  Inaert-0  THEN 

Bbl".0668»(Oi/L  )»Reeg»Preg 

Bb2-l+.04»< <Qi/L>»Reeg»Preg>".6666 

0mega-3.66+<8bl/8b2 ) 

ENO    IF 

IF  Inaert-l  THEN 

Onega-< 1+5 . 484E-3»PregA .7»( Res/Hod >* 1 .25  >* .5 

ENO  IF 

IF  Insert-2  THEN 

Onega-(Reeg" .65  >*< Preg" .46 ) 

ENO  IF 

IF  Insert-3  THEN 

Owega-( Reeg* .76  )•( Preg" .46  ) 

ENO  IF 

R-114/R113  Prooerties 
Hfg»FNHfg( Tsat  ) 
Kf-FNK( Tf  ) 
Rhof-FNRho<  Tf ) 
Muf«FNMu( Tf ) 

F-(Kf3»Rhof-2»G»Hfg/(Muf  •Oo«Qdp  )  >"  .33333 

Ho»Alpa«F 

Tuo-Tsat-Qdo/Ho 

Tf-T5at/3+:»Tuo/3 

Y-< 1/Uo-Rn )»F 

X»0o*F/( Keg*0mega  ) 

PRINT    "OMEGA-" ; Omega ; "F«" .F. "X-" .X; "Y-" ; Y 
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3365  Xa(Jset)-X    i   INEFFICIENT  (MODIFY  LATER) 

3370  YaUsetw 

3375  Sx-Sx+X 

3380  Sy-Sy+Y 

338S  Sx3«Sxs+X»X 

3390  Sxy-Sxy+X»Y 

3395  Jset-Jset>l 

3400  IF  Jset<N3eta  THEN  3135 

3405  ASSIGN  QFile  TO  • 

3410  Slope-(Nset3»Sxy-Sx»Sy  )/<N3et3»Sxs-Sx"2 ) 

3415  Intcpt-(Sy-Slope»Sx  )/Nsets 

3420  Cic-1/Slope 

3425  Alpac-1/Intcpt 

3430  Cerr- ASS< (Ci-Cic )/Cic ) 

3435  Aerr- A8S< ( Alpac-Alpa >/Alpac  > 

3440  IF  Cerr>.001  OR  Aerr,. 001  THEN 

3445  Alpa«< Alpa+Alpac )*.S 

3450  Ci«<Cx+Cic)«.S 

3455  !         PRINT  "CIO"  ;Cic  ;  "ALPA-" ;  Alpa 

3460  SOTO  3100 

3465  ENO  IF 

3470  SEEP 

3475  BEEP 

3480  •  .     PRINTER  IS  1 

3485  PRINT  ' C IC- "; C ic ;' ALPA- "; Alpa 

3490  Cia< Insert )-Cic 

3495  Alpaa< I  )-Alpac 

3500  PRINTER  IS  701 

3505  FOR  J-0  TO  Naet3-1 

3510  OUTPUT  9Xy ,Xa( J  ),Ya( J  ) 

3515  NEXT  J 

3520  PRINTER  IS  1 

3525  NEXT  I 

3530  ASSIGN  ®Xy  TO  • 

3535  SUBENO 

8270  OEF  FNFmcaK  T,Pc  > 

8280  Mdotl5-.0019»Pc-.0098 

8290  Mdot0».00206«Pc-. 00414 

8300  Mdot24-.00225«Pc-. 00329 

8310  IF  T'.'B  THEN 

8311  Mdt-< (T+15)*<ndot0-Mdotl5)/l5  )+Mdotl5 
8320  ELSE 

8470  Mdt-<  T»<  Mdot24-Mdot0  )/24  )-t-Mdot0 

8471  ENO  IF 
8480  RETURN  Mdt 
8490  FNENO 

8500  OEF  FNCa3h<X ) 

8510  P-EXP( X  ) 

85Z0  Q»EXP( -X  ) 

8530  Co3h-.5»(P+Q  ) 

8540  RETURN  Cosh 

8550  FNENO 
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